OFFSET
1,2
PROG
(Python)
import math
for i in range(4000):
s = (i*i) ^ ((i+1)*(i+1))
t = math.isqrt(s)
if s == t*t:
print(i, end=', ')
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Alex Ratushnyak, Mar 27 2013
STATUS
approved
