{
	i = 0;
	l = 1;
	for (k=1, oo,
		if (omega(k) <= 1,
			l = lcm(l, k);
			fordiv (l/k, d,
				print (i++ " " k*d);
			);
			if (i >= 10 000,
				break(2)
			)
		)
	)
}

quit