STATUS
editing
approved
editing
approved
nmax = 100; A063725 = Rest[CoefficientList[Series[(EllipticTheta[3, 0, x] - 1)^2/4, {x, 0, nmax}], x]]; s = 1; Do[s *= Sum[Binomial[A063725[[k]], j]*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]]; , {k, 2, nmax}]; Take[CoefficientList[s, x], nmax + 1] (* Vaclav Kotesovec, Nov 09 2018 *)
approved
editing
proposed
approved
editing
proposed
Seiichi Manyama, <a href="/A321428/b321428.txt">Table of n, a(n) for n = 0..10000</a>
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed