大数分解Pollard_rho模板

最新推荐文章于 2021-11-26 11:29:54 发布

alan_cty

最新推荐文章于 2021-11-26 11:29:54 发布

阅读量678

点赞数 1

CC 4.0 BY-SA版权

分类专栏：模板模板库文章标签： Pollardrho millerabin 模板

本文链接：https://blog.youkuaiyun.com/alan_cty/article/details/78462621

模板同时被 2 个专栏收录

21 篇文章

订阅专栏

模板库

21 篇文章

订阅专栏

本文介绍了一个基于Pollard_rho算法实现的质因数分解模板，并提供了包括快速幂、欧几里得算法在内的相关辅助函数。此外还实现了Miller-Rabin素性检测算法，用以判断数的素性。

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

发现我之前写的Pollard_rho模板是假的。。。
学(抄)了一发真的，但跑的贼慢

ll mult(ll x,ll y,ll Mo) {
    x%=Mo;y%=Mo;
    ll tmp=(ll)((long double)x*y/Mo+1e-8)*Mo;
    return (x*y-tmp+Mo)%Mo;
}

ll pwr(ll x,ll y,ll Mo) {
    ll z=1;x%=Mo;y%=(Mo-1);
    for(;y;y>>=1,x=mult(x,x,Mo))
        if (y&1) z=mult(z,x,Mo);
    return z;
}

ll gcd(ll x,ll y) {
    if (!x) return 1;
    if (x<0) return gcd(-x,y);
    return y?gcd(y,x%y):x;
}

bool check(ll a,ll n,ll x,ll t) {
    ll res=pwr(a,x,n),lst=res;
    fo(i,1,t) {
        res=mult(res,res,n);
        if (res==1&&lst!=1&&lst!=n-1) return 1;
        lst=res;
    }
    if (res!=1) return 1;
    else return 0;
} 

bool Miller_Rabin(ll n) {
    if (n<2) return 0;
    if (n==2) return 1;
    if (!(n&1)) return 0;
    ll x=n-1,t=0;
    while (!(x&1)) x>>=1,t++;
    fo(i,1,S) {
        ll a=rand()%(n-1)+1;
        if (check(a,n,x,t)) return 0;
    }
    return 1;
}

ll Pollard_rho(ll n,ll t) {
    int i=1,k=2;
    ll x0=rand()%n,y=x0;
    while (1) {
        i++;
        x0=(mult(x0,x0,n)+t)%n;
        ll d=gcd(y-x0,n);
        if(d!=1&&d!=n) return d;  
        if(y==x0) return n;  
        if(i==k) y=x0,k<<=1;  
    }  
}

void find(ll n) {
    if (n==1) return;
    if (Miller_Rabin(n)) {
        pri[++tot]=n;
        return;
    }
    ll p=n;
    while (p>=n) p=Pollard_rho(p,rand()%(n-1)+1);
    find(p);find(n/p);
}