diff options
Diffstat (limited to 'math/rho.cpp')
| -rw-r--r-- | math/rho.cpp | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/math/rho.cpp b/math/rho.cpp index 865a438..df70251 100644 --- a/math/rho.cpp +++ b/math/rho.cpp @@ -3,7 +3,7 @@ ll rho(ll n) { // Findet Faktor < n, nicht unbedingt prim. if (n % 2 == 0) return 2; ll x = 0, y = 0, prd = 2; auto f = [n](lll x){return (x * x) % n + 1;}; - for (ll t = 30, i = n / 2 + 7; t % 40 || gcd(prd, n) == 1; t++) { + for (ll t = 30, i = n/2 + 7; t % 40 || gcd(prd, n) == 1; t++) { if (x == y) x = ++i, y = f(x); if (ll q = (lll)prd * abs(x-y) % n; q) prd = q; x = f(x); y = f(f(y)); @@ -15,6 +15,5 @@ void factor(ll n, map<ll, int>& facts) { if (n == 1) return; if (isPrime(n)) {facts[n]++, return;} ll f = rho(n); - factor(n / f, facts); - factor(f, facts); + factor(n / f, facts); factor(f, facts); } |