login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A194471 E.g.f. A(x) satisfies: A(x) = exp(x) + x*A(x)^2. 3
1, 2, 9, 79, 1065, 19401, 445933, 12389021, 403897553, 15120448273, 639345572181, 30138682861365, 1567316344601593, 89137628104427033, 5503952108613407933, 366697176991277153341, 26220726323043177903009, 2002962250253424509250081 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

The radius of convergence r of the e.g.f. satisfies: r = exp(-r)/4 = limit (n+1)*a(n)/a(n+1) = 0.203888354702240... with A(r) = 1/(2*r) = 2.452322501352287...

LINKS

Table of n, a(n) for n=0..17.

FORMULA

E.g.f.: A(x) = (1 - sqrt(1 - 4*x*exp(x))) / (2*x).

a(n) = 1 + n*Sum_{k=0..n-1} C(n-1,k)*a(k)*a(n-1-k) for n>=0.

a(n) ~ sqrt(2)*sqrt(1+LambertW(1/4))*n^(n-1)/(4*exp(n)*LambertW(1/4)^(n+1)). - Vaclav Kotesovec, Aug 19 2013

a(n) = n!*sum(k=0..n, (k+1)^(n-k-1)*binomial(2*k,k)/(n-k)!). - Vladimir Kruchinin, Sep 01 2014

EXAMPLE

E.g.f.: A(x) = 1 + 2*x + 9*x^2/2! + 79*x^3/3! + 1065*x^4/4! +...

Related expansion:

A(x)^2 = 1 + 4*x + 26*x^2/2! + 266*x^3/3! + 3880*x^4/4! + 74322*x^5/5! +...

Illustrate the recurrence:

a(2) = 1 + 2*(1*1*2 + 1*2*1) = 1 + 2*4 = 9;

a(3) = 1 + 3*(1*1*9 + 2*2*2 + 1*9*1) = 1 + 3*26 = 79;

a(4) = 1 + 4*(1*1*79 + 3*2*9 + 3*9*2 + 1*79*1) = 1 + 4*266 = 1065;

a(5) = 1 + 5*(1*1*1065 + 4*2*79 + 6*9*9 + 4*79*2 + 1*1065*1) = 1 + 5*3880 = 19401.

MATHEMATICA

f[0] = 1; f[n_] := f[n] = 1 + n*Sum[ Binomial[n - 1, k]*f[k]*f[n - 1 - k] , {k, 0, n - 1}]; Array[f, 18, 0] (* Robert G. Wilson v, Aug 25 2011 *)

PROG

(PARI) {a(n)=n!*polcoeff((1 - sqrt(1 - 4*x*exp(x +O(x^(n+2))))) / (2*x), n)}

(PARI) {a(n)=1+n*sum(k=0, n-1, binomial(n-1, k)*a(k)*a(n-1-k))}

(Maxima) a(n):=n!*sum((k+1)^(n-k-1)*binomial(2*k, k)/(n-k)!, k, 0, n); /* Vladimir Kruchinin, Sep 01 2014 */

CROSSREFS

Sequence in context: A111196 A229211 A056918 * A215629 A221460 A122720

Adjacent sequences:  A194468 A194469 A194470 * A194472 A194473 A194474

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Aug 24 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 2 09:21 EDT 2020. Contains 335398 sequences. (Running on oeis4.)