poj1181 大数分解

 //Accepted    164 KB    422 ms
 //类似poj2429 大数分解
 #include <cstdio>
 #include <cstring>
 #include <ctime>
 #include <time.h>
 #include <iostream>
 #include <algorithm>
 using namespace std;
 ;
 __int64 gcd(__int64 a,__int64 b)
 {
     ) return a;
     return gcd(b,a%b);
 }
 __int64 mult_mod(__int64 a,__int64 b,__int64 p)
 {
     __int64 res=,temp=a%p;
     while (b)
     {
         )
         {
             res+=temp;
             if (res>=p) res-=p;
         }
         temp<<=;
         if (temp>=p) temp-=p;
         b>>=;
     }
     return res;
 }
 __int64 exp_mod(__int64 a,__int64 b,__int64 p)
 {
     __int64 res=,exp=a%p;
     )
     {
         )
         res=mult_mod(res,exp,p);
         exp=mult_mod(exp,exp,p);
         b>>=;
     }
     return res;
 }
 bool miller_rabin(__int64 n,__int64 times)
 {
     ) return true;
      || !(n&)) return false;
     __int64 a,u=n-,x,y;
     ;
     ==)
     {
         t++;
         u/=;
     }
     srand(time());
     ;i<times;i++)
     {
         a=rand()%(n-)+;
         x=exp_mod(a,u,n);
         ;j<t;j++)
         {
             y=mult_mod(x,x,n);
              && x!= && x!=n-)
             return false;
             x=y;
         }
         ) return false;
     }
     return true;
 }
 __int64 pollar_rho(__int64 n,int c)
 {
     __int64 x,y,d,i=,k=;
     srand(time());
     x=rand()%(n-)+;
     y=x;
     while (true)
     {
         i++;
         x=(mult_mod(x,x,n)+c)%n;
         d=gcd(y-x,n);
          && d<n) return d;
         if (y==x) return n;
         if (i==k)
         {
             y=x;
             k<<=;
         }
     }
 }
 __int64 min_ans;
 void findMinFactor(__int64 n,int c)
 {
     ) return ;
     ))
     {
         if (n<min_ans) min_ans=n;
         return ;
     }
     __int64 p=n;
     while (p>=n)
     p=pollar_rho(p,c--);
     findMinFactor(p,c);
     findMinFactor(n/p,c);
 }
 int main()
 {
     int T;
     scanf("%d",&T);
     while (T--)
     {
         __int64 n;
         scanf("%I64d",&n);
         )==true)
         {
             printf("Prime\n");
         }
         else
         {
             min_ans=inf;
             findMinFactor(n,);
             printf("%I64d\n",min_ans);
         }
     }
     ;
 }