Wieferich primes: Difference between revisions
Content added Content deleted
(add fermat) |
(add parigp) |
||
Line 534: | Line 534: | ||
<pre>Wieferich primes less than 5000: |
<pre>Wieferich primes less than 5000: |
||
1093 |
1093 |
||
3511</pre> |
|||
=={{header|PARI/GP}}== |
|||
<lang parigp>iswief(p)=if(isprime(p)&&(2^(p-1)-1)%p^2==0,1,0) |
|||
for(N=1,5000,if(iswief(N),print(N)))</lang> |
|||
{{out}}<pre>1093 |
|||
3511</pre> |
3511</pre> |
||