OFFSET
0,3
FORMULA
a(n) = Sum_{k=0..n} binomial(n-1,k-1) * A007837(k) * n! / k!.
MATHEMATICA
nmax = 20; CoefficientList[Series[Product[(1 + x^k/(k! (1 - x)^k)), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
PROG
(PARI) seq(n)={Vec(serlaplace(prod(k=1, n, (1 + x^k/(k!*(1 - x)^k)) + O(x*x^n))))} \\ Andrew Howroyd, Feb 13 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 13 2020
STATUS
approved