RosettaCodeData/Task/Arithmetic-Rational/PARI-GP/arithmetic-rational.parigp

6 lines
66 B
Plaintext

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