{
	seen = 0;

	a = vector(10 000);

	for (n=1, #a,
		a[n] = 2;
		while (bittest(seen, a[n]) || (n%2==0 && gcd(a[n],a[n/2])==1),
			a[n]++;
		);

		print (n " " a[n]);

		seen += 2^a[n];
	);
}

quit