OFFSET
0,11
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: Product_{k>0} (1 + x^k)^A025426(k).
MAPLE
N:= 100: # for a(0)..a(N)
P:= 1:
for i from 1 to floor(sqrt(N)) do
for j from i while i^2 + j^2 <= N do
P:= P * (1 + x^(i^2 + j^2))
od od:
S:= series(P, x, N+1):
seq(coeff(S, x, k), k=0..N); # Robert Israel, Apr 21 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 09 2018
STATUS
approved