seen = 0
unseen = 1

K = 10

other(p) = my (v=unseen); while (bittest(seen, v) || hammingweight(p+v)<K, v++); return (v)

for (n=1, 2^13, a = if(n==1, 1, other(p)); seen += 2^a; while (bittest(seen, unseen), unseen++); print (n " " a); p=a)

quit