OFFSET
1,3
PROG
(PARI) an=0; for(i=1, 10^5, an=bitxor(an, i^2); issquare(an)&print1(i", "))
(Python)
import math
x = 0
for i in range(1<<20):
x ^= i*i
t = int(math.sqrt(x))
if x == t*t:
print(str(i), end=', ')
# Alex Ratushnyak, Mar 27 2013
CROSSREFS
KEYWORD
easy,base,nonn
AUTHOR
M. F. Hasler, Oct 20 2008
EXTENSIONS
More terms from Alex Ratushnyak, Mar 27 2013
STATUS
approved