OFFSET
1,1
COMMENTS
A025426(a(n)) > 2. - Reinhard Zumkeller, Feb 26 2015
LINKS
Robert Israel and Reinhard Zumkeller, Table of n, a(n) for n = 1..10000 (first 2508 terms from Robert Israel)
MAPLE
N:= 100000: # generate all entries <=N
SSQ:= {}: SSQ2:= {}: SSQ3:= {}:
for a from 1 to floor(sqrt(N)) do
for b from a to floor(sqrt(N-a^2)) do
n:= a^2 + b^2;
if member(n, SSQ2) then SSQ3:= SSQ3 union {n}
elif member(n, SSQ) then SSQ2:= SSQ2 union {n}
else SSQ:= SSQ union {n}
end if
end do end do:
SSQ3; # Robert Israel, Jan 20 2013
MATHEMATICA
okQ[n_] := Length[Select[PowersRepresentations[n, 2, 2][[All, 1]], Positive] ] > 2; Select[Range[5000], okQ] (* Jean-François Alcover, Mar 04 2019 *)
PROG
(Haskell)
a025294 n = a025294_list !! (n-1)
a025294_list = filter ((> 2) . a025426) [1..]
-- Reinhard Zumkeller, Feb 26 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved