A302863 a(n) = [x^(n^2)] (1 + theta_3(x))^n/(2^n*(1 - x)), where theta_3() is the Jacobi theta function.

1, 2, 6, 29, 165, 1203, 9763, 83877, 793049, 7903501, 83570177, 933697153, 10905583809, 133352809334, 1695473999478, 22354920990148, 305096197935075, 4296142551821184, 62336908825014452, 930284705538262688, 14255992611680074754, 224065160215526683317, 3607018540134004189466

Offset: 0

Author

Ilya Gutkovskiy, Apr 14 2018

Keywords

• nonn

Comments

a(n) = number of nonnegative solutions to (x_1)^2 + (x_2)^2 + ... + (x_n)^2 <= n^2.

Links

• Eric Weisstein's World of Mathematics, Jacobi Theta Functions
• Index entries for sequences related to sums of squares

Crossrefs

Main diagonal of A302998.

Cf. A000603, A000604, A010052, A055403, A055404, A055405, A055406, A055407, A055408, A055409, A298329, A302861, A302862.

Programs

Mathematica

Table[SeriesCoefficient[(1 + EllipticTheta[3, 0, x])^n/(2^n (1 - x)), {x, 0, n^2}], {n, 0, 22}]
Table[SeriesCoefficient[1/(1 - x) Sum[x^k^2, {k, 0, n}]^n, {x, 0, n^2}], {n, 0, 22}]