diff options
| author | Gloria Mundi <gloria@gloria-mundi.eu> | 2024-11-16 01:24:14 +0100 |
|---|---|---|
| committer | Gloria Mundi <gloria@gloria-mundi.eu> | 2024-11-16 01:24:14 +0100 |
| commit | 98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch) | |
| tree | 5113d5cc24d1ad5f93810b6442ce584a36950dc8 /math/millerRabin.cpp | |
| parent | ad3856a6b766087df0036de0b556f4700a6498c9 (diff) | |
| parent | 8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff) | |
mzuenni tests
Diffstat (limited to 'math/millerRabin.cpp')
| -rw-r--r-- | math/millerRabin.cpp | 19 |
1 files changed, 0 insertions, 19 deletions
diff --git a/math/millerRabin.cpp b/math/millerRabin.cpp deleted file mode 100644 index cb27d29..0000000 --- a/math/millerRabin.cpp +++ /dev/null @@ -1,19 +0,0 @@ -constexpr ll bases32[] = {2, 7, 61}; -constexpr ll bases64[] = {2, 325, 9375, 28178, 450775, - 9780504, 1795265022}; -bool isPrime(ll n) { - if (n < 2 || n % 2 == 0) return n == 2; - ll d = n - 1, j = 0; - while (d % 2 == 0) d /= 2, j++; - for (ll a : bases64) { - if (a % n == 0) continue; - ll v = powMod(a, d, n); //with mulmod or int128 - if (v == 1 || v == n - 1) continue; - for (ll i = 1; i <= j; i++) { - v = ((lll)v * v) % n; - if (v == n - 1 || v <= 1) break; - } - if (v != n - 1) return false; - } - return true; -} |
