OFFSET
0,2
FORMULA
a(0) = 1; a(n) = n*a(n-1) + n^(2*n).
PROG
(PARI) a(n) = n!*sum(k=0, n, k^(2*k)/k!);
(PARI) a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=i*v[i]+i^(2*i)); v;
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 23 2022
STATUS
approved