OFFSET
0,8
COMMENTS
x + x^4 + x^9 + x^16 + x^25 + ... is the expansion of (theta_3(0, x) - 1)/2, where theta_3 is the Jacobi theta function.
LINKS
Max Alekseyev, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Series Reversion
FORMULA
For n>1, a(n) = Sum_{j2,j3,...} (-1)^(j2+j3+...) * (n-1+j2+j3+...)! / (j2!*j3!*...) / n!, where the sum is taken over all nonnegative integers j2, j3, ... such that (2^2-1)*j2 + (3^2-1)*j3 + ... = n-1. - Max Alekseyev, Jul 06 2021
MATHEMATICA
InverseSeries[(EllipticTheta[3, 0, x] - 1)/2 + O[x]^30][[3]]
PROG
(PARI) Vec( serreverse( sum(i=1, 32, x^i^2) + O(x^33^2) ) ); \\ Max Alekseyev, Jul 06 2021
CROSSREFS
KEYWORD
easy,sign
AUTHOR
Vladimir Reshetnikov, Jul 09 2015
STATUS
approved