%I #15 May 21 2024 09:40:52
%S 4,16,0,0,32,0,0,0,16,32,0,0,0,0,0,0,32,32,0,0,0,0,0,0,0,96,0,0,0,0,0,
%T 0,0,0,0,0,32,0,0,0,64,0,0,0,0,0,0,0,0,48,0,0,64,0,0,0,0,0,0,0,0,0,0,
%U 0,64,0,0,0,0,0,0,0,32,64,0,0,0,0,0,0,32,32,0,0,0,0,0,0,0,64,0,0,0,0,0,0,0,32,0,0,96
%N a(n) = r(n)*r(n+1), where r(n) = A004018(n) is the number of ways of writing n as a sum of two squares.
%C a(n)=0 unless n == 0, 1 or 4 (mod 8).
%D H. Iwaniec. Spectral methods of automorphic forms, volume 53 of Graduate Studies in Mathematics. American Mathematical Society, Providence, RI, 2002.
%H Robert Israel, <a href="/A330315/b330315.txt">Table of n, a(n) for n = 0..10000</a>
%H Fernando Chamizo, <a href="https://doi.org/10.2969/jmsj/05110237">Correlated sums of r(n)</a>, J. Math. Soc. Japan, 51(1):237-252, 1999.
%H Fernando Chamizo, and Roberto J. Miatello, <a href="https://arxiv.org/abs/1812.10725">Sums of squares in real quadratic fields and Hilbert modular groups</a>, arXiv preprint arXiv:1812.10725 [math.NT], 2018.
%p N:= 200: # for a(0)..a(N)
%p g1:= 1 + 2*add(x^(i^2),i=1..floor(sqrt(N+1))):
%p g2:= expand(g1^2):
%p R:= [seq(coeff(g2,x,i),i=0..N+1)]:
%p seq(R[i]*R[i+1],i=1..N+1); # _Robert Israel_, Jun 12 2020
%t a[n_] := SquaresR[2, n] SquaresR[2, n + 1]; a /@ Range[0, 100] (* _Giovanni Resta_, Jun 12 2020 *)
%Y Cf. A004018, A330316, A330317, A330318.
%K nonn
%O 0,1
%A _N. J. A. Sloane_, Dec 11 2019