OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..5000
MATHEMATICA
nmax = 40; CoefficientList[Series[Product[(1 + k*x^k)^(k^2), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 40; s = 1 + x; Do[s*=Sum[Binomial[k^2, j]*k^j*x^(j*k), {j, 0, Floor[nmax/k] + 1}]; s = Select[Expand[s], Exponent[#, x] <= nmax &]; , {k, 2, nmax}]; CoefficientList[s, x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Apr 15 2017
STATUS
approved