poj 1811 Prim test

基本上一个裸的Miller_Rabin大素数判定和一个裸的Pollard_rho素数分解算法，当模板用吧！

 #include<cstdio>
 #include<algorithm>
 #include<stdlib.h>
 #define LL long long
 using namespace std;
 const int s=;
 int tol;
 LL factor[];

 LL mult_mod(LL a,LL b,LL c)//计算a*b%c;
 {
     LL ret=;
     a%=c;
     b%=c;
     while(b>)
     {
         if(b&) ret=(ret+a)%c;
         a<<=;
         if(a>=c) a%=c;
         b>>=;
     }
     return ret;
 }

 LL pow_mod(LL a,LL b,LL mod)//计算a^b%b
 {
     if(b==) return a%mod;
     a%=mod;
     LL tmp=a;
     LL ret=;
     while(b>)
     {
         if(b&) ret=mult_mod(ret,tmp,mod);
         tmp=mult_mod(tmp,tmp,mod);
         b>>=;
     }
     return ret;
 }

 //以a为基,n-1=x*2^t      a^(n-1)=1(mod n)  验证n是不是合数
 //一定是合数返回true,不一定返回false
 bool check(LL a,LL n,LL x,LL t)
 {
     LL ret=pow_mod(a,x,n);
     LL last=ret;
     for(int i=; i<=t; i++)
     {
         ret=mult_mod(ret,ret,n);
         if(ret==&&last!=&&last!=n-) return ;
         last=ret;
     }
     if(ret!=) return ;
     return false;
 }

 // Miller_Rabin()算法素数判定
 //是素数返回true.(可能是伪素数，但概率极小)
 //合数返回false;
 bool Miller_Rabin(LL a)
 {
     if(a<) return ;
     if(a==) return ;
     if((a&==)) return ;
     LL x=a-;
     LL t=;
     while((x&)==)
     {
         x>>=;
         t++;
     }
     for(int i=; i<s; i++)
     {
         long long b=rand()%(a-)+;
         if(check(b,a,x,t))
             return ;
     }
     return ;
 }

 LL gcd(LL a,LL b)
 {
     if(a==)return ;
     if(a<) return gcd(-a,b);
     while(b)
     {
         LL t=a%b;
         a=b;
         b=t;
     }
     return a;
 }

 LL Pollard_rho(LL x,LL c)
 {
     LL i=,k=;
     LL x0=rand()%x;
     LL y=x0;
     while()
     {
         i++;
         x0=(mult_mod(x0,x0,x)+c)%x;
         LL d=gcd(y-x0,x);
         if(d!=&&d!=x) return d;
         if(y==x0) return x;
         if(i==k)
         {
             y=x0;
             k+=k;
         }
     }
 }

 void findfac(LL n)
 {
     if(Miller_Rabin(n))//素数
     {
         factor[tol++]=n;
         return;
     }
     LL p=n;
     while(p>=n)p=Pollard_rho(p,rand()%(n-)+);
     findfac(p);
     findfac(n/p);
 }

 int main()
 {
     LL a;
     int n;
     scanf("%d",&n);
     while(n--)
     {
         scanf("%lld",&a);
         if(Miller_Rabin(a))
         {
             puts("Prime");
             continue;
         }
         tol=;
         findfac(a);
         LL ans=factor[];
         for(int i=;i<tol;i++)
           if(factor[i]<ans)
              ans=factor[i];
         printf("%lld\n",ans);
     }
     return ;
 }