OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..411
Eric Weisstein's World of Mathematics, Lambert W-Function.
FORMULA
E.g.f.: 4*LambertW(-x/2*exp(x/2))^2 / (x^2*exp(x)).
E.g.f.: exp( L(x) ) where L(x) = -2*LambertW(-x*exp(x/2)/2) is the e.g.f. of A100526.
a(n) ~ sqrt(1+LambertW(exp(-1)))*n^(n-1)/(2^(n-1)*exp(n-2)*LambertW(exp(-1))^n). - Vaclav Kotesovec, Jun 23 2016
From Seiichi Manyama, Feb 11 2023: (Start)
E.g.f. satisfies A(x) = exp( x * ( exp(x) * A(x) )^(1/2) ).
a(n) = (1/2^(n-1)) * Sum_{k=0..n} k^(n-k) * (k+2)^(k-1) * binomial(n,k). (End)
EXAMPLE
E.g.f.: A(x) = 1 + x + 3*x^2/2! + 13*x^3/3! + 77*x^4/4! + 581*x^5/5! + 5347*x^6/6! + 58213*x^7/7! + 732937*x^8/8! + 10487737*x^9/9! + 168217811*x^10/10! + 2990748509*x^11/11! + 58397418037*x^12/12! +...
such that
A(x) = 1 + x*exp(x/2)*A(x)^(1/2) + x^2/2!*exp(x)*A(x) + x^3/3!*exp(3*x/2)*A(x)^(3/2) + x^4/4!*exp(2*x)*A(x)^2 + x^5/5!*exp(5*x/2)*A(x)^(5/2) + x^6/6!*exp(3*x)*A(x)^3 +...
The logarithm of A(x) begins:
log(A(x)) = x + 2*x^2/2! + 6*x^3/3! + 28*x^4/4! + 180*x^5/5! + 1476*x^6/6! + 14728*x^7/7! + 173216*x^8/8! + 2346480*x^9/9! + 35981200*x^10/10! + 616111056*x^11/11! + 11652662880*x^12/12! +...+ A100526(n)*x^n/n! +...
which equals -2*LambertW(-x*exp(x/2)/2).
MATHEMATICA
CoefficientList[Series[4*LambertW[-x/2*E^(x/2)]^2 / (x^2*E^x), {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, Jun 23 2016 *)
PROG
(PARI) {a(n) = my(A=1+x); for(i=1, n, A = sum(m=0, n, x^m/m!*exp(m/2*x +x*O(x^n))*A^(m/2)) ); n!*polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
(PARI) a(n) = sum(k=0, n, k^(n-k)*(k+2)^(k-1)*binomial(n, k))/2^(n-1); \\ Seiichi Manyama, Feb 11 2023
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Paul D. Hanna, Jun 14 2016
STATUS
approved