OFFSET
0,2
COMMENTS
Binomial transform of A005651.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..440
N. J. A. Sloane, Transforms
FORMULA
E.g.f.: exp(x + Sum_{k>=1} Sum_{j>=1} x^(j*k)/(k*(j!)^k)).
a(n) = Sum_{k=0..n} binomial(n,k)*A005651(k).
a(n) ~ exp(1) * A247551 * n!. - Vaclav Kotesovec, Jul 21 2019
MAPLE
seq(coeff(series(factorial(n)*exp(x)*mul((1-x^k/factorial(k))^(-1), k=1..n), x, n+1), x, n), n = 0 .. 22); # Muniru A Asiru, Oct 15 2018
MATHEMATICA
nmax = 22; CoefficientList[Series[Exp[x] Product[1/(1 - x^k/k!), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
nmax = 22; CoefficientList[Series[Exp[x + Sum[Sum[x^(j k)/(k (j!)^k), {j, 1, nmax}], {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Binomial[n, k] Total[Apply[Multinomial, IntegerPartitions[k], {1}]], {k, 0, n}], {n, 0, 22}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 15 2018
STATUS
approved