OFFSET
0,4
COMMENTS
Also number of ways to write n as an ordered sum of n nonnegative numbers one less than a square.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..2000
FORMULA
a(n) = [x^(2*n)] ((theta_3(x) - 1) / 2)^n, where theta_3() is the Jacobi theta function.
a(n) = [x^n] (Sum_{k>=0} x^(k*(k + 2)))^n.
a(n) = A337165(2n,n). - Alois P. Heinz, Feb 04 2021
MAPLE
b:= proc(n, t) option remember;
`if`(n=0, `if`(t=0, 1, 0), `if`(t<1, 0, add((s->
`if`(s>n, 0, b(n-s, t-1)))(j^2), j=1..isqrt(n))))
end:
a:= n-> b(2*n, n):
seq(a(n), n=0..39); # Alois P. Heinz, Feb 04 2021
MATHEMATICA
Table[SeriesCoefficient[((EllipticTheta[3, 0, x] - 1)/2)^n, {x, 0, 2 n}], {n, 0, 37}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 31 2021
STATUS
approved