(Translated by https://www.hiragana.jp/)
A341429 - OEIS
login
A341429
Number of positive solutions to (x_1)^2 + (x_2)^2 + ... + (x_10)^2 <= n^2.
7
56, 1108, 12098, 84624, 439371, 1785368, 6078048, 18139393, 48586117, 118929400, 270250734, 578320470, 1169522013, 2261784392, 4193751331, 7509793133, 13008356489, 21921125415, 35951569269, 57666975238, 90464266824, 139295784464, 210514511189, 313228848537
OFFSET
4,1
FORMULA
a(n) is the coefficient of x^(n^2) in expansion of (theta_3(x) - 1)^10 / (1024 * (1 - x)).
MAPLE
b:= proc(n, k) option remember; `if`(k=0, 1, `if`(n=0, 0,
add((s->`if`(s>n, 0, b(n-s, k-1)))(j^2), j=1..isqrt(n))))
end:
a:= n-> b(n^2, 10):
seq(a(n), n=4..27); # Alois P. Heinz, Feb 11 2021
MATHEMATICA
Table[SeriesCoefficient[(EllipticTheta[3, 0, x] - 1)^10/(1024 (1 - x)), {x, 0, n^2}], {n, 4, 27}]
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 11 2021
STATUS
approved