OFFSET
0,2
FORMULA
G.f.: Product_{k>=1} (1+x^k)^(k^2+1).
a(0) = 1; a(n) = (-1/n) * Sum_{k=1..n} ( Sum_{d|k} (-1)^(k/d) * d * (d^2+1) ) * a(n-k).
PROG
(PARI) my(N=40, x='x+O('x^N)); Vec(prod(k=1, N, (1+x^k)^(k^2+1)))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Jun 10 2023
EXTENSIONS
Name suggested by Joerg Arndt, Jun 11 2023
STATUS
approved