OFFSET
1,2
COMMENTS
All terms except for a(2)=2 are odd.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..1000
MATHEMATICA
s={1}; m=1; Do[f=0; Do[If[!SquareFreeQ[s[[i]]^2+p^2], f=1; Break[]], {i, m}]; If[f<1, AppendTo[s, p]; m++], {p, 2, 10^3}]; s
PROG
(PARI) v=List([1]); for(m=2, 1e3, for(j=1, #v, if(issquare(m^2+v[j]^2), next(2))); listput(v, m)); Vec(v) \\ Charles R Greathouse IV, Jul 30 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jul 28 2012
STATUS
approved