summaryrefslogtreecommitdiff
path: root/math/millerRabin.cpp
diff options
context:
space:
mode:
authorGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
committerGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
commit98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch)
tree5113d5cc24d1ad5f93810b6442ce584a36950dc8 /math/millerRabin.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'math/millerRabin.cpp')
-rw-r--r--math/millerRabin.cpp19
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;
-}