POJ1811 Prime Test （Pollard-Rho算法）

最新推荐文章于 2024-10-16 21:33:04 发布

Sunshine_cfbsl

最新推荐文章于 2024-10-16 21:33:04 发布

阅读量785

点赞数 1

CC 4.0 BY-SA版权

文章标签： poj 数论质因数分解素数判定

本文链接：https://blog.youkuaiyun.com/Sunshine_cfbsl/article/details/52517077

数论同时被 3 个专栏收录

7 篇文章

订阅专栏

POJ

4 篇文章

订阅专栏

素数判定

2 篇文章

订阅专栏

本文介绍了一种解决大规模素数分解问题的有效方法——Pollard-Rho算法，并提供了详细的算法实现代码。通过Miller-Rabin素性测试辅助，该算法能够高效地找出合数的非平凡因子。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

显然，这么大的数据不能用朴素算法过。
需要用到一种叫做Pollard-Rho的算法。
Pollard-Rho算法戳这里

看懂了就变成了裸题。。。
直接上代码：

#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;

const int MAXN = 65;
long long x[MAXN], ans;

long long multi(long long a, long long b, long long p) {
    long long ans = 0;
    while(b) {
        if(b&1LL) ans = (ans+a)%p;
        a = (a+a)%p;
        b >>= 1;
    }
    return ans;
}

long long qpow(long long a, long long b, long long p) {
    long long ans = 1;
    while(b) {
        if(b&1LL) ans = multi(ans, a, p);
        a = multi(a, a, p);
        b >>= 1;
    }
    return ans;
}

bool Miller_Rabin(long long n) {
    if(n == 2) return true;
    int s = 20, i, t = 0;
    long long u = n-1;
    while(!(u & 1)) {
        t++;
        u >>= 1;
    }
    while(s--) {
        long long a = rand()%(n-2)+2;
        x[0] = qpow(a, u, n);
        for(i = 1; i <= t; i++) {
            x[i] = multi(x[i-1], x[i-1], n);
            if(x[i] == 1 && x[i-1] != 1 && x[i-1] != n-1) return false;
        }
        if(x[t] != 1) return false;
    }
    return true;
}

long long gcd(long long a, long long b) {
    return b ? gcd(b, a%b) : a;
}

long long Pollard_Rho(long long n, int c) {
    long long i = 1, k = 2, x = rand()%(n-1)+1, y = x;
    while(true) {
        i++;
        x = (multi(x, x, n) + c)%n;
        long long p = gcd((y-x+n)%n, n);
        if(p != 1 && p != n) return p;
        if(y == x) return n;
        if(i == k) {
            y = x;
            k <<= 1;
        }
    }
}

void find(long long n, int c) {
    if(n == 1) return;
    if(Miller_Rabin(n)) {
        ans = min(ans, n);
        return;
    }
    long long p = n, k = c;
    while(p >= n) p = Pollard_Rho(p, c--);
    find(p, k);
    find(n/p, k);
}

int main() {
    long long n, i;
    int t;
    scanf("%d", &t);
    while(t--) {
        ans = 1LL<<62;
        scanf("%lld", &n);
        find(n, 120);
        if(ans == n) printf("Prime\n");
        else printf("%lld\n", ans);
    }
    return 0;
}