s = 0 unseen = 1 seen(v) = bittest(s, v) see(v) = s = bitor(s, 2^v); while (seen(unseen), unseen++) cc = [0] other(p) = { my (w=1, b, c); for (e=0, oo, if (!bittest(p, e), if (#cc==w, cc=concat(cc, vector(#cc)); ); b = 2^e; for (k=1, w, c = cc[w+k] = b+cc[k]; if (!seen(c) && gcd(p,c)==1, return (c); ); ); w*=2; ); ); } for (n=1, 10 000, see(v=other(n)); print (n " " v)) quit