a = vector(16383)
for (n=1, #a, pa = n; fordiv(n, d, pa = bitxor( pa, a[d]) ) ; a[n] = pa ; print(n " " a[n]))
quit