|
%I #14 Apr 21 2024 11:39:43
%S 1,0,1,0,0,1,0,1,1,0,2,0,1,2,0,3,0,2,3,1,4,1,2,4,1,6,3,4,6,2,7,5,6,8,
%T 5,9,7,9,10,9,12,10,13,14,13,18,13,19,17,18,25,19,28,24,25,33,26,36,
%U 35,33,46,35,47,48,44,61,48,62,65,60,78,68,79,87,79,101,93
%N Expansion of Product_{1 <= i <= j} (1 + x^(i^2 + j^2)).
%H Robert Israel, <a href="/A321435/b321435.txt">Table of n, a(n) for n = 0..10000</a>
%F G.f.: Product_{k>0} (1 + x^k)^A025426(k).
%p N:= 100: # for a(0)..a(N)
%p P:= 1:
%p for i from 1 to floor(sqrt(N)) do
%p for j from i while i^2 + j^2 <= N do
%p P:= P * (1 + x^(i^2 + j^2))
%p od od:
%p S:= series(P,x,N+1):
%p seq(coeff(S,x,k),k=0..N); # _Robert Israel_, Apr 21 2024
%Y Cf. A025426, A321423, A321428, A321436, A321437.
%K nonn
%O 0,11
%A _Seiichi Manyama_, Nov 09 2018
|
