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!)
A215003 E.g.f. satisfies: A(x) = x + A(x)^2*exp(A(x))/2. 2
1, 1, 6, 51, 610, 9360, 175371, 3881269, 99083916, 2866173795, 92650139560, 3309877247886, 129495296417565, 5506601964331741, 252882022593123210, 12473008516763238375, 657621285983402967766, 36908182475893682882532, 2196924777991866630353571 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

LINKS

Table of n, a(n) for n=1..19.

FORMULA

E.g.f.: Series_Reversion(x - x^2*exp(x)/2).

E.g.f.: x + Sum_{n>=1} d^(n-1)/dx^(n-1) exp(n*x)/2^n*x^(2*n) / n!.

E.g.f.: x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) exp(n*x)/2^n*x^(2*n-1) / n! ).

O.g.f.: Sum_{n>=0} (2*n)!/n!/2^n * x^(n+1) / (1 - n*x)^(2*n+1).

a(n) = Sum_{k=0..n-1} k^(n-k-1)/(n-k-1)! * (n+k-1)!/k!/2^k.

a(n) ~ r*sqrt((1+r)/(2+r*(4+r))) * n^(n-1) / (exp(n) * (r*(1+r)/(2+r))^n), where r = 0.4912251835444738797155054345... is the root of the equation r*(2+r)*exp(r) = 2. - Vaclav Kotesovec, Dec 28 2013

EXAMPLE

E.g.f.: A(x) = x + x^2/2! + 6*x^3/3! + 51*x^4/4! + 610*x^5/5! + 9360*x^6/6! +...

where A(x - x^2*exp(x)/2) = x and A(x) = x + A(x)^2*exp(A(x))/2.

Related expansions:

A(x)^2 = 2*x^2/2! + 6*x^3/3! + 54*x^4/4! + 630*x^5/5! + 9570*x^6/6! +...

exp(A(x)) = 1 + x + 2*x^2/2! + 10*x^3/3! + 85*x^4/4! + 1011*x^5/5! + 15466*x^6/6! + 289213*x^7/7! +...

A(x) = x + exp(x)*x^2/2 + d/dx exp(2*x)*x^4/2!/4 + d^2/dx^2 exp(3*x)*x^6/3!/8 + d^3/dx^3 exp(4*x)*x^8/4!/16 +...

log(A(x)/x) = exp(x)*x/2 + d/dx exp(2*x)*x^3/2!/4 + d^2/dx^2 exp(3*x)*x^5/3!/8 + d^3/dx^3 exp(4*x)*x^7/4!/16 +...

Ordinary Generating Function:

O.g.f.: x + x^2 + 6*x^3 + 51*x^4 + 610*x^5 + 9360*x^6 + 175371*x^7 +...

O.g.f.: x + 2*x^2/(1-x)^3/2 + 6*2!*x^3/(1-2*x)^5/4 + 20*3!*x^4/(1-3*x)^7/8 + 70*4!*x^5/(1-4*x)^9/16 + 252*5!*x^6/(1-5*x)^11/32 +...+ (2*n)!/n!*x^(n+1)/(1-n*x)^(2*n+1)/2^n +...

MATHEMATICA

Rest[CoefficientList[InverseSeries[Series[x - x^2*E^x/2, {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Dec 28 2013 *)

PROG

(PARI) {a(n)=sum(k=0, n-1, k^(n-k-1)/(n-k-1)! * (n+k-1)!/k!/2^k )}

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

for(n=1, 25, print1(a(n), ", "))

(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}

{a(n)=local(A=x); A=x+sum(m=1, n, Dx(m-1, exp(m*x+x*O(x^n))/2^m*x^(2*m)/m!)); n!*polcoeff(A, n)}

(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}

{a(n)=local(A=x+x^2+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, exp(m*x+x*O(x^n))/2^m*x^(2*m-1)/m!)+x*O(x^n))); n!*polcoeff(A, n)}

(PARI) /* O.g.f.: */

{a(n)=polcoeff(sum(m=0, n, (2*m)!/m!/2^m*x^(m+1)/(1-m*x+x*O(x^n))^(2*m+1)), n)}

CROSSREFS

Cf. A213643.

Sequence in context: A246189 A293128 A304185 * A134525 A125865 A210810

Adjacent sequences:  A215000 A215001 A215002 * A215004 A215005 A215006

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Jul 30 2012

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 January 21 05:52 EST 2021. Contains 340333 sequences. (Running on oeis4.)