diff options
| author | Paul Jungeblut <paul.jungeblut@gmail.com> | 2016-10-15 23:30:44 +0200 |
|---|---|---|
| committer | Paul Jungeblut <paul.jungeblut@gmail.com> | 2016-10-15 23:30:44 +0200 |
| commit | 463e389ea499dcc26314196c014a85ee9b124847 (patch) | |
| tree | 64fc77a603df88c4d898fba452c60006194c818a /math/discreteLogarithm.cpp | |
| parent | 53c8e56d8b0ee3b4374ab90630673b971ddf2710 (diff) | |
Adding primitive root and discrete logarithm.
Diffstat (limited to 'math/discreteLogarithm.cpp')
| -rw-r--r-- | math/discreteLogarithm.cpp | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/math/discreteLogarithm.cpp b/math/discreteLogarithm.cpp new file mode 100644 index 0000000..6d3f656 --- /dev/null +++ b/math/discreteLogarithm.cpp @@ -0,0 +1,13 @@ +// Bestimmt Lösung x für a^x=b mod m. +ll solve (ll a, ll b, ll m) { // Laufzeit: O(sqrt(m)*log(m)) + ll n = (ll)sqrt((double)m) + 1; + map<ll,ll> vals; + for (int i = n; i >= 1; i--) vals[powMod(a, i * n, m)] = i; + for (int i = 0; i <= n; i++) { + ll cur = (powMod(a, i, m) * b) % m; + if (vals.count(cur)) { + ll ans = vals[cur] * n - i; + if (ans < m) return ans; + }} + return -1; +} |
