a(n) = my (v=1, i=0, x=1); \ for (d=2, oo, \ if (n==0, return (v)); \ if (gcd(x,d)==1 && omega(d)==1, \ if (bittest(n, i), \ v *= d; \ n -= 2^i; \ x *= d; \ ); \ i++; \ ); \ ); for (n=0, 10 000, print (n " " a(n))) quit