OFFSET
1,3
COMMENTS
MATHEMATICA
Select[Range[0, 15000], OddQ[Sqrt[8*BitXor[#^2, (#(#+1))/2]+1]]&] (* Harvey P. Dale, Jul 22 2024 *)
PROG
(Python)
import math
for n in range(100000000):
a = (n*n) ^ (n*(n+1)/2)
r = int(math.sqrt(a*2))
if r*(r+1)==a*2: print n,
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Alex Ratushnyak, Jun 08 2013
STATUS
approved