login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A138014 E.g.f. A(x) satisfies exp(A(x)) = x + exp(A(x)^2) where A(0) = 0. 6
1, 1, 2, 16, 174, 1988, 27124, 453136, 8791980, 191869392, 4668291000, 125662750464, 3706032771336, 118759029538368, 4109063510399088, 152696171895135744, 6065376023980289424, 256455323932682550528, 11499944141042532006432, 545124523779848580648960 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..100

FORMULA

E.g.f.: A(x) = Series_Reversion( exp(x) - exp(x^2) ).

a(n) = -sum(k=1..n-1, (k^n/n!+(k^(n/2)*(-1)^k*((-1)^n+1))/(2*(n/2)!)+sum(j=1..k-1, (-1)^(k-j)*binomial(k,j)*sum(m=1..n, (j^(2*m-n)*(k-j)^(n-m)*binomial(m,n-m))/m!)))*a(k)), n>1, a(1)=1. - Vladimir Kruchinin, Jun 25 2011

a(n) ~ n^(n-1) * sqrt((2*s-1)/(2-2*s+4*s^2)) / (exp(1+s)*(1-1/(2*s)))^n, where s = 0.6310764773894916166238... is the root of the equation exp(s) = 2*s*exp(s^2). - Vaclav Kotesovec, Jan 06 2014

MATHEMATICA

Rest[CoefficientList[InverseSeries[Series[E^x-E^(x^2), {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Jan 06 2014 *)

PROG

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

(Maxima)

  a(n):=if n=1 then 1 else -sum((k^n/n!+(k^(n/2)*(-1)^k*((-1)^n+1))/(2*(n/2)!)+sum((-1)^(k-j)*binomial(k, j)*sum((j^(2*m-n)*(k-j)^(n-m)*binomial(m, n-m))/m!, m, 1, n), j, 1, k-1))*a(k), k, 1, n-1); /* Vladimir Kruchinin, Jun 25 2011 */

CROSSREFS

Sequence in context: A215702 A155659 A108999 * A206988 A217360 A052606

Adjacent sequences:  A138011 A138012 A138013 * A138015 A138016 A138017

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Feb 27 2008

EXTENSIONS

a(19)-a(20) from Vincenzo Librandi, Feb 20 2018

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 26 07:15 EST 2022. Contains 350573 sequences. (Running on oeis4.)