OFFSET
0,4
LINKS
FORMULA
a(n) = [x^(n*(n+1)/2)] (Sum_{k>=1} x^(k^2))^n.
EXAMPLE
a(4) = 6 because fourth triangular number is 10 and we have [4, 4, 1, 1], [4, 1, 4, 1], [4, 1, 1, 4], [1, 4, 4, 1], [1, 4, 1, 4] and [1, 1, 4, 4].
MAPLE
b:= proc(n, t) option remember; local i; if n=0 then
`if`(t=0, 1, 0) elif t<1 then 0 else 0;
for i while i^2<=n do %+b(n-i^2, t-1) od; % fi
end:
a:= n-> b(n*(n+1)/2, n):
seq(a(n), n=0..25); # Alois P. Heinz, Feb 05 2018
MATHEMATICA
Table[SeriesCoefficient[(-1 + EllipticTheta[3, 0, x])^n/2^n, {x, 0, n (n + 1)/2}], {n, 0, 25}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 01 2018
STATUS
approved