Perfect numbers: Difference between revisions
→{{header|PARI/GP}}
Line 2,822:
<syntaxhighlight lang="parigp">isPerfect(n)=sigma(n,-1)==2</syntaxhighlight>
Show perfect numbers
<syntaxhighlight lang="parigp">
forprime(p=2, 2281,
if(isprime(2^p-1),
print(p"\t",(2^p-1)*2^(p-1))))
</syntaxhighlight>
===Faster with Lucas-Lehmer test===
<syntaxhighlight lang="parigp">p=2;n=3;n1=2;
|