OFFSET
1,1
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000 (first 98 terms from Robert Price)
EXAMPLE
27625 is in the sequence via 20^2 + 165^2 = 27^2 + 164^2 = 45^2 + 160^2 = 60^2 + 155^2 = 83^2 + 144^2 = 88^2 + 141^2 = 101^2 + 132^2 = 115^2 + 120^2. - David A. Corneth, Jun 01 2025
MATHEMATICA
nn = 114985; 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, _?(# >= 8 &)]] (* T. D. Noe, Apr 07 2011 *)
PROG
(PARI) upto(n) = {my(v = vector(n)); for(i = 1, sqrtint(n), i2 = i^2; for(j = i, sqrtint(n - i^2), v[i2 + j^2]++)); select(x->x >= 8, v, 1)} \\ David A. Corneth, Jun 01 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
