Pollard-Rho大整数拆分模板

随机拆分，简直机智。

关于过程可以看http://wenku.baidu.com/link?url=JPlP8watmyGVDdjgiLpcytC0lazh4Leg3s53WIx1_Pp_Y6DJTC8QkZZqmiDIxvgFePUzFJ1KF1G5xVVAoUZpxdw9GN-S46eVeiJ6Q-zXdei

看完后，觉得随机生成数然后和n计算gcd，可以将随机的次数根号一下。思想很叼。

对于里面说的birthday trick，在执行次数上我怎么看都只能减一半。只是把平均分布，变成了靠近0的的分布。

不过怎么说，这个好像是大家都公认比较靠谱的。 所以，我就勉强相信了。

/*****************************
大整数拆分模板（long long范围内）
调用Divide(n,222);
返回的结果在divsor中，因子最小值为dmi
注意：复杂度为n^(1/4),多次调用初始化dcnt,dmi
*****************************/
 
#define INF 1e18
 
long long divsor[];
int dcnt=;
long long dmi=INF;
 
//输入一个long long 范围内的素数，是素数返回true，否则返回false。定义检测次数TIMES,错误率为(1/4)^TIMES
#define TIMES 10
 
long long GetRandom(long long n)
{
    //cout<<RAND_MAX<<endl;
    long long num = (((unsigned long long)rand() + )*rand())%n;
    return num+;
}
 
long long Mod_Mul(long long a,long long b,long long mod)
{
    long long msum=;
    while(b)
    {
        if(b&) msum = (msum+a)%mod;
        b>>=;
        a = (a+a)%mod;
    }
    return msum;
}
 
long long Quk_Mul(long long a,long long b,long long mod)
{
    long long qsum=;
    while(b)
    {
        if(b&) qsum=Mod_Mul(qsum,a,mod);
        b>>=;
        a=Mod_Mul(a,a,mod);
    }
    return qsum;
}
 
bool Miller_Rabin(long long n)
{
    if(n==||n==||n==||n==||n==) return true;
    if(n==||n%==||n%==||n%==||n%==||n%==) return false;
    int div2=;
    long long tn=n-;
    while( !(tn%) )
    {
        div2++;
        tn/=;
    }
    for(int tt=;tt<TIMES;tt++)
    {
        long long x=GetRandom(n-); //随机得到[1,n-1]
        if(x==) continue;
        x=Quk_Mul(x,tn,n);
        long long pre=x;
        for(int j=;j<div2;j++)
        {
            x = Mod_Mul(x, x, n);
            if(x==&&pre!=&&pre!=n-) return false;
            pre=x;
        }
        if(x!=) return false;
    }
    return true;
}
 
long long gcd(long long a,long long b)
{
    if(b==) return a;
    return gcd(b,a%b);
}
 
long long pollard_rho(long long dn,long long dc)
{
    long long x,y,d,i=,k=;
    x = GetRandom(dn-);
    y = x;
    while()
    {
        i++;
        x = (Mod_Mul(x, x, dn) + dc)%dn;
        d = gcd( y-x , dn );
        if( < d && d < dn )
            return d;
        if( y==x ) return dn;
        if( i==k )
        {
            y=x;
            k <<= ;
        }
    }
}
 
void Divide(long long dn,int dk)
{
    if(dn==) return ;
    if( Miller_Rabin(dn) == true )
    {
        divsor[dcnt++]=dn;
        dmi = min(dmi,dn);
        return ;
    }
    long long dtmp=dn;
    while(dtmp>=dn) dtmp = pollard_rho(dtmp,dk--);//随机寻找dn的因子，dtmp
    Divide(dtmp, dk);
    Divide(dn/dtmp,dk);
}
 
/*
int main() {
    int T;
    cin>>T;
    while(T--)
    {
        long long n;
        cin>>n;
        if( Miller_Rabin(n) ) printf("Prime\n");
        else
        {
            dmi=INF;
            dcnt=0;
            Divide(n,222);
            cout<<dmi<<endl;
        }
    }
    return 0;
}
*/