OFFSET
0,12
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..700
FORMULA
E.g.f.: Product_{k>=5} (1 + x^k/k!).
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1,
`if`(i>n, 0, b(n, i+1)+binomial(n, i)*b(n-i, i+1)))
end:
a:= n-> b(n, 5):
seq(a(n), n=0..31); # Alois P. Heinz, Apr 28 2021
MATHEMATICA
nmax = 31; CoefficientList[Series[Product[(1 + x^k/k!), {k, 5, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
a[0] = 1; a[n_] := a[n] = -(n - 1)! Sum[DivisorSum[k, # (-#!)^(-k/#) &, # > 4 &] a[n - k]/(n - k)!, {k, 1, n}]; Table[a[n], {n, 0, 31}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 28 2021
STATUS
approved