OFFSET
0,3
COMMENTS
Logarithm of e.g.f. yields the e.g.f. of A304866.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..385
FORMULA
E.g.f. A(x) = Sum_{n>=0} a(n) * x^n / n! satisfies:
(1) Sum_{n>=0} (n*x - log(A(x)))^n / n! = 1.
(2) Sum_{n>=0} (n*x - p*log(A(x)))^n / n! = (1 + LambertW(-x))^(p-1).
(3) Sum_{n>=0} ((n + m)*x - log(A(x)))^n / n! = ( LambertW(-x)/(-x) )^m.
(4) Sum_{n>=0} ((n + m)*x - p*log(A(x)))^n / n! = ( LambertW(-x)/(-x) )^m * (1 + LambertW(-x))^(p-1).
a(n) ~ 2^(1/2 - exp(-1)/2) * sqrt(Pi) * n^(n + exp(-1)/2 - 1/2) / Gamma(exp(-1)/2) * (1 + sqrt(2)*log(n)*Gamma(exp(-1)/2 + 1) / (sqrt(n)*Gamma(exp(-1)/2 - 1/2))). - Vaclav Kotesovec, Sep 01 2020
EXAMPLE
E.g.f.: A(x) = 1 + x + 2*x^2/2! + 9*x^3/3! + 70*x^4/4! + 760*x^5/5! + 10451*x^6/6! + 173159*x^7/7! + 3350796*x^8/8! + 74100408*x^9/9! + 1842574557*x^10/10! + ...
RELATED SERIES.
log(A(x)) = x + x^2/2! + 5*x^3/3! + 40*x^4/4! + 434*x^5/5! + 5921*x^6/6! + 97152*x^7/7! + 1861224*x^8/8! + 40757712*x^9/9! + ... + A304866(n)*x^n/n! + ...
such that Sum_{n>=0} (n*x - log(A(x)))^n / n! = 1.
MATHEMATICA
nmax = 20; CoefficientList[Series[(1 + LambertW[-x])^(x/LambertW[-x]), {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Sep 01 2020 *)
PROG
(PARI) {a(n) = my(W = serreverse(-x*exp(x +x*O(x^n)))); n!*polcoeff( (1 + W)^(x/W), n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 18 2018
STATUS
approved