Wieferich primes: Difference between revisions

m
→‎{{header|Perl}}: better efficiency
(Added Sidef)
m (→‎{{header|Perl}}: better efficiency)
Line 386:
{{libheader|ntheory}}
<lang perl>use feature 'say';
use ntheory 'qw(is_prime' powmod);
use bignum;
use ntheory 'is_prime';
 
say 'Wieferich primes less than 5000: ' . join ', ', grep { is_prime($_) and not powmod(2, (2**($_-1) -1) %, $_**2$_) )== 1 } 1..5000;</lang>
{{out}}
<pre>Wieferich primes less than 5000: 1093, 3511</pre>
2,747

edits