OFFSET
1,1
COMMENTS
"Primitive" means that x and y are coprime in the representations x^2+y^2.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 10^4: # to get all terms <= N
V:= Vector(N):
for a from 1 to floor(sqrt(N)) do
for b from 1 to min(a, floor(sqrt(N-a^2))) do
if igcd(a, b) > 1 then next fi;
r:= a^2 + b^2;
V[r]:= V[r]+1;
od od:
select(n -> V[n] > 2, [$1..N]); # Robert Israel, Feb 07 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
R. J. Mathar, Feb 01 2017
STATUS
approved