login
A025300
Numbers that are the sum of 2 nonzero squares in 9 or more ways.
5
71825, 93925, 122525, 138125, 143650, 156325, 160225, 173225, 187850, 204425, 209525, 223925, 226525, 235625, 244205, 245050, 257725, 267325, 273325, 276250, 287300, 292825, 296225, 300625, 308125, 308425, 312650, 320450, 333125, 337025
OFFSET
1,1
MAPLE
N:= 4*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] >= 9, [$1..N]); # Robert Israel, Jun 01 2025
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
STATUS
approved