Wieferich primes: Difference between revisions
→{{header|Nim}}
(add APL solution) |
|||
Line 490:
println(filter(p -> (big"2"^(p - 1) - 1) % p^2 == 0, primes(5000))) # [1093, 3511]
</lang>
=={{header|Mathematica}}/{{header|Wolfram Language}}==
<lang Mathematica>ClearAll[WieferichPrimeQ]
WieferichPrimeQ[n_Integer] := PrimeQ[n] && Divisible[2^(n - 1) - 1, n^2]
Select[Range[5000], WieferichPrimeQ]</lang>
{{out}}
<pre>{1093, 3511}</pre>
=={{header|Nim}}==
| |||