s = 0 unseen = 1 seen(v) = bit test(s, v) see(v) = s = bit or(s, 2^v); while (seen(unseen), unseen++) other(p) = { see(p); while (bigomega(unseen)!=2, unseen++; ); for (v=unseen, oo, if (!seen(v) && bigomega(v)==2 && bigomega(p+v)==2, return (v); ); ); } for (n=1, 10 000, print (n " " v=if (n==1, 4, other(v)))) quit