OFFSET
0,2
FORMULA
a(0) = 1; a(n) = (n-1)! * Sum_{k=1..n} A033286(k) * a(n-k)/(n-k)!.
MAPLE
a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*
ithprime(j)*j*binomial(n, j)*j!, j=1..n)/n)
end:
seq(a(n), n=0..20); # Alois P. Heinz, Apr 28 2022
MATHEMATICA
a[0] = 1; a[n_] := a[n] = (n-1)! Sum[k Prime[k] a[n-k]/(n-k)!, {k, 1, n}];
Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Dec 28 2022 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(sum(k=1, N, prime(k)*x^k))))
(PARI) a(n) = if(n==0, 1, (n-1)!*sum(k=1, n, k*prime(k)*a(n-k)/(n-k)!));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 28 2022
STATUS
approved