OFFSET
1,1
LINKS
MAPLE
N:= 10^5: # generate all entries <=N
V:= Vector(N, datatype=integer[4]):
for a from 1 to floor(sqrt(N)) do
for b from a do
n:= a^2 + b^2;
if n > N then break fi;
V[n]:= V[n]+1
od od:
select(t -> V[t] >= 5, [$1..N]); # Robert Israel, Jun 01 2025
MATHEMATICA
nn = 30000; t = Table[0, {nn}]; lim = Floor[Sqrt[nn - 1]]; Do[num = i^2 + j^2; If[num <= nn, t[[num]]++], {i, lim}, {j, i}]; Flatten[Position[t, _?(# >= 5 &)]] (* T. D. Noe, Apr 07 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
