OFFSET
1,1
FORMULA
a(n) = A009003(n)^2.
EXAMPLE
25 = 5^2 = 4^2 + 3^2, and so 25 is in the sequence.
100 = 10^2 = 8^2 + 6^2, and so 100 is in the sequence.
169 = 13^2 = 12^2 + 5^2, and so 169 is in the sequence.
MATHEMATICA
c = {}; Do[Do[k = a^2 + b^2; If[IntegerQ[Sqrt[k]], AppendTo[c, k]], {a, 1, b - 1}], {b, 200}]; Union[c] (* Artur Jasinski *)
Select[Range[100]^2, Length[PowersRepresentations[#, 2, 2]] > 1 &] (* Alonso del Arte, Feb 11 2014 *)
PROG
(PARI) select(n->for(k=1, sqrtint(n\2), if(issquare(n-k^2), return(n>k^2))); 0, vector(100, i, i^2)) \\ Charles R Greathouse IV, Jul 02 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 25 2007
STATUS
approved