login
A341402
Number of nonnegative solutions to (x_1)^2 + (x_2)^2 + ... + (x_7)^2 <= n.
6
1, 8, 29, 64, 106, 169, 281, 422, 548, 702, 961, 1276, 1556, 1864, 2326, 2893, 3390, 3852, 4545, 5455, 6253, 7002, 8080, 9361, 10453, 11496, 12903, 14618, 16194, 17643, 19589, 22011, 24027, 25714, 28143, 31188, 33792, 36137, 39203, 42920, 46294, 49108, 52580, 57165, 61365
OFFSET
0,2
COMMENTS
Partial sums of A045849.
FORMULA
G.f.: (1 + theta_3(x))^7 / (128 * (1 - x)).
a(n^2) = A055406(n).
MAPLE
b:= proc(n, k) option remember; `if`(n=0, 1, `if`(n<0 or k<1, 0,
b(n, k-1)+add(b(n-j^2, k-1), j=1..isqrt(n))))
end:
a:= proc(n) option remember; b(n, 7)+`if`(n>0, a(n-1), 0) end:
seq(a(n), n=0..44); # Alois P. Heinz, Feb 10 2021
MATHEMATICA
nmax = 44; CoefficientList[Series[(1 + EllipticTheta[3, 0, x])^7/(128 (1 - x)), {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 10 2021
STATUS
approved