0,4

Robert Israel, Table of n, a(n) for n = 0..2000

Index entries for Gaussian integers and primes

N:= 100: A:= Array(0..N):

for a from 1 to N do

for b from 0 to a-1 do

r:= ceil(sqrt(a^2+b^2));

if r > N then break fi;

A[r]:= A[r]+1;

od od:

seq(A[i], i=0..N); # Robert Israel, Dec 18 2018

Cf. A036700.

Sequence in context: A244230 A277814 A006162 * A053756 A210436 A003004

Adjacent sequences: A036698 A036699 A036700 * A036702 A036703 A036704

nonn

Clark Kimberling

approved