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>