Miller Robin素数测试与Pollcard Rho因数分解

原文地址

<1>预备算法

ll mul(ll a,ll b,ll p)
{  
    ll rn=0,i;
    for(i=1; i<=b; i<<=1,a=(a+a)%p)
    if(b&i) rn=(rn+a)%p;
    return rn;
}//计算模意义下两大数乘积
ll ksm(ll a,ll b,ll p)
{
    ll rn=1;
    for(;b;a=mul(a,a,p),b>>=1)
      if(b&1) rn=mul(rn,a,p);
    return rn;
}//计算模意义下两大数乘方
ll gcd(ll a, ll b)
{
    ll tmp;
    if(a<b) tmp=a,a=b,b=tmp;
    while(b) tmp=a%b,a=b,b=tmp;
    return a;
}//求最大公约数

<2>Miller Robin素数测试

能在$(0.25)^S$的错误率下判定质数，相应地需要付出$O(S*log n)$的时间复杂度。

原理基于两个定理：

1.若p是质数，$0<a<p$，那么$a^{(p-1)}≡1(mod\quad p)$。 （可惜对于某些合数也成立，为此我们需要借助定理2

2.对于$0<x<p$，$x^2≡1(mod\quad p)$有且只有两解：$x=1$和$x=p-1$。 （关于这个的判定，可以在定理1的判定中做快速幂的时候顺便判下

基本思路：对于一切2以下的数和其他偶数，先特判掉。然后把$(p-1)$拆成$k*2^j$，先令$x=(p-1)^k$，对于$2^j$，我们一边平方x，一边利用定理2判定。最后用定理1检查一次。

bool isprime(ll n)
{
    if(n==2) return true;
    if(n<2||!(n&1)) return false;
    ll a,x,y,u=n-1;
    int t=0;
    while((u&1)==0) t++,u>>=1;
    for(i=0;i<S;i++)
    {
        a=rand()%(n-1)+1;
        x=ksm(a,u,n);
        for(int j=1; j<=t; j++)
        {
            y=mul(x,x,n);
            if(y==1&&x!=1&&x!=n-1) return false;
            x=y;
        }
        if(x!=1) return false;
    }
    return true;
}

<3>Pollcard Rho因数分解

复杂度貌似是$n^{\frac 1 4}$，这应该是因数分解的最快算法了。

void rho(ll n)  
{  
    if(isprime(n)){sta[++top]=n;return;}  
    ll x,y,z,c,d;
    int i,j;  
    for(;;)
    {
        x=rand()*rand()%(n-1)+1;
        c=rand()*rand()%(n-1)+1;
        for(y=x,i=j=2;;i++)
        {
            x=(mul(x,x,n)+c)%n;
            z=x-y;
            if(z<0) z=-z;
            d=gcd(z,n);
            if(d>1&&d<n){rho(d);rho(n/d);return;}
            if(x==y) break;
            if(i==j) y=x,j<<=1;
        }
    }
}