OFFSET
1,1
COMMENTS
This sequence is complete - there are no other terms.
LINKS
EXAMPLE
666 = 15^2 + 21^2, therefore 666 is a term.
PROG
(Magma) lst:=[]; for e in [1..3] do a:=Truncate((10^e-1)/9); for x in [1..9] do n:=a*x; f:=Factorization(n); if IsSquare(n) and not IsZero([x: x in [1..#f] | f[x][1] mod 4 eq 1]) or not IsSquare(n) and IsZero([x: x in [1..#f] | IsOdd(f[x][2]) and f[x][1] mod 4 eq 3]) then Append(~lst, n); end if; end for; end for; lst;
CROSSREFS
KEYWORD
nonn,base,fini,full
AUTHOR
Arkadiusz Wesolowski, Feb 15 2017
STATUS
approved
