Wieferich primes: Difference between revisions
m
no edit summary
(Added Wren) |
mNo edit summary |
||
Line 32:
V{ 1093 3511 }
</pre>
=={{header|Julia}}==
<lang julia>using Primes
is_weiferich(p) = (big"2"^(p - 1) - 1) % p^2 == 0
function weiferich_to(N)
n = 0
while (n = nextprime(n + 1)) < N
is_weiferich(n) && print(n, " ")
end
end
weiferich_to(5000) # prints 1093 3511
</lang>
=={{header|Raku}}==
|