Wieferich primes: Difference between revisions
Added Perl
(Added C++ solution) |
SqrtNegInf (talk | contribs) (Added Perl) |
||
Line 100:
println(filter(p -> (big"2"^(p - 1) - 1) % p^2 == 0, primes(5000))) # [1093, 3511]
</lang>
=={{header|Perl}}==
{{libheader|ntheory}}
<lang perl>use feature 'say';
use bignum;
use ntheory 'is_prime';
say 'Weiferich primes less than 5000: ' . join ', ', grep { is_prime($_) and not ( (2**($_-1) -1) % $_**2 ) } 1..5000;</lang>
{{out}}
<pre>Weiferich primes less than 5000: 1093, 3511</pre>
=={{header|Phix}}==
| |||