OFFSET
0,2
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
FORMULA
a(n) ~ exp(3 * Pi^(1/3) * ((sqrt(2)-1) * Zeta(3/2))^(2/3) * n^(1/3) / 2) * sqrt(2/3) * ((sqrt(2)-1) * Zeta(3/2) / Pi)^(1/3) / (4*n^(5/6)). - Vaclav Kotesovec, Dec 09 2016
MATHEMATICA
nmax = 100; CoefficientList[Series[Product[(1 + x^(k^2))^2, {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 20; poly = ConstantArray[0, nmax^2 + 1]; poly[[1]] = 1; poly[[2]] = 2; poly[[3]] = 1; Do[Do[Do[poly[[j + 1]] += poly[[j - k^2 + 1]], {j, nmax^2, k^2, -1}]; , {p, 1, 2}], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Dec 09 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Dec 08 2016
STATUS
approved