|
|
A305304
|
|
Expansion of e.g.f. 1/(1 + LambertW(-x/(1 + x))).
|
|
1
|
|
|
1, 1, 2, 9, 52, 405, 3786, 42301, 542984, 7924041, 129110230, 2327399481, 45940938924, 986045445853, 22856850513602, 569163515043285, 15150885843083536, 429364157810169105, 12905794670246364078, 410108007771441394129, 13736898888997174964660, 483740530150449507164901, 17866185834825657429606682
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n-1,k-1)*k^k*n!/k!.
a(n) ~ n^n * (exp(1) - 1)^(n - 1/2) / exp(n - 1/2). - Vaclav Kotesovec, Aug 18 2018
|
|
MAPLE
|
a:=series(1/(1+LambertW(-x/(1+x))), x=0, 23): seq(n!*coeff(a, x, n), n=0..22); # Paolo P. Lava, Mar 26 2019
|
|
MATHEMATICA
|
nmax = 22; CoefficientList[Series[1/(1 + LambertW[-x/(1 + x)]), {x, 0, nmax}], x] Range[0, nmax]!
Join[{1}, Table[Sum[(-1)^(n - k) Binomial[n - 1, k - 1] k^k n!/k!, {k, n}], {n, 22}]]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|