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}}==
4,107

edits