OFFSET
0,2
LINKS
Eric Weisstein's World of Mathematics, Bell Polynomial
FORMULA
a(n) = n! * [x^n] exp(n*(exp(x) + x - 1)).
a(n) = Sum_{k=0..n} binomial(n,k) * BellPolynomial_k(n) * n^(n-k).
a(n) ~ c * exp((r^2/(1-r) - 1)*n) * n^n / (1-r)^n, where r = A333761 = 0.59894186245845296434937... is the root of the equation LambertW(r) = 1-r and c = 0.897950293373062982395233981707095204244165706668136925178217032608352851... - Vaclav Kotesovec, Jun 09 2020
MATHEMATICA
Table[n! SeriesCoefficient[Exp[n (Exp[x] + x - 1)], {x, 0, n}], {n, 0, 17}]
Join[{1}, Table[Sum[Binomial[n, k] BellB[k, n] n^(n - k), {k, 0, n}], {n, 1, 17}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 19 2020
STATUS
approved