A025300 Numbers that are the sum of 2 nonzero squares in 9 or more ways.

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

Author

David W. Wilson

Keywords

• nonn

Links

• Robert Israel, Table of n, a(n) for n = 1..10000
• Index entries for sequences related to sums of squares

Crossrefs

Cf. A025292, A025299, A025301, A025337.

Programs

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 *)