OFFSET
0,3
FORMULA
Expansion of e.g.f. Sum_{k>=0} x^k / (k! * (1 - k*x^k/k!)).
a(n) = n! * Sum_{d|n} 1/(d * (d-1)!^(n/d)) for n > 0.
a(p) = 1 + p! for prime p.
MATHEMATICA
a[n_]:= n! * DivisorSum[n, 1/(# * (# - 1)!^(n/#)) &]; a[0] = 1; Array[a, 23, 0] (* Amiram Eldar, Aug 22 2022 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(sum(k=0, N, x^k/(k!-k*x^k))))
(PARI) a(n) = if(n==0, 1, n!*sumdiv(n, d, 1/(d*(d-1)!^(n/d))));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 22 2022
STATUS
approved