for(n=2,1<<19,
  s=0;
  fordiv(n,d,s+=1/d);
  if(s==2,print(n))
)