Wieferich primes: Difference between revisions

Content added Content deleted

Inline

@@ Line 386: / 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 ( (2**($_-1) -1) % $_**2 ) } 1..5000;</lang>
+say 'Wieferich primes less than 5000: ' . join ', ', grep { is_prime($_) and powmod(2, $_-1, $_*$_) == 1 } 1..5000;</lang>
 {{out}}
 <pre>Wieferich primes less than 5000: 1093, 3511</pre>