OFFSET
0,3
FORMULA
a(n) = Sum_{k=0..n} binomial(n-1,k-1)*A000272(k)*n!/k!.
a(n) ~ (1 + exp(-1))^(n + 3/2) * n^(n-2). - Vaclav Kotesovec, Jan 27 2019
MAPLE
seq(n!*coeff(series(1-LambertW(-x/(1-x))*(2+LambertW(-x/(1-x)))/2, x=0, 21), x, n), n=0..20); # Paolo P. Lava, Jan 29 2019
MATHEMATICA
nmax = 20; CoefficientList[Series[1 - LambertW[-x/(1 - x)] (2 + LambertW[-x/(1 - x)])/2, {x, 0, nmax}], x] Range[0, nmax]!
Join[{1}, Table[Sum[Binomial[n - 1, k - 1] k^(k - 2) n!/k!, {k, n}], {n, 20}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 27 2019
STATUS
approved