OFFSET
1,1
MATHEMATICA
nn = 337025; 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, _?(# >= 9 &)]] (* T. D. Noe, Apr 07 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved