hdu 3864 素数分解

题意：求n是否只有4个因子，如果是就输出除1外的所有因子。

模板题，就不排版了

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<map>
 #include<ctime>
 using namespace std;
 #define MOD 1000000007
 const int INF=0x3f3f3f3f;
 const double eps=1e-;
 #define cl(a) memset(a,0,sizeof(a))
 #define ts printf("*****\n");
 const int MAXN=;
 int n,m,tt;
 const int S = ; //随机算法判定次数，一般8~10就够了
 // 计算ret = (a*b)%c a,b,c < 2^63
 long long mult_mod(long long a,long long b,long long c)
 {
 a %= c;
 b %= c;
 long long ret = ;
 long long tmp = a;
 while(b)
 {
 if(b & )
 {
 ret += tmp;
 if(ret > c)ret -= c;//直接取模慢很多
 }
 tmp <<= ;
 if(tmp > c)tmp -= c;
 b >>= ;
 }
 return ret;
 }
 // 计算 ret = (a^n)%mod
 long long pow_mod(long long a,long long n,long long mod)
 {
 long long ret = ;
 long long temp = a%mod;
 while(n)
 {
 if(n & )ret = mult_mod(ret,temp,mod);
 temp = mult_mod(temp,temp,mod);
 n >>= ;
 }
 return ret;
 }
 // 通过 a^(n-1)=1(mod n)来判断n是不是素数
 // n-1 = x*2^t 中间使用二次判断
 // 是合数返回true, 不一定是合数返回false
 bool check(long long a,long long n,long long x,long long t)
 {
 long long ret = pow_mod(a,x,n);
 long long last = ret;
 for(int i = ;i <= t;i++)
 {
 ret = mult_mod(ret,ret,n);
 if(ret ==  && last !=  && last != n-)return true;//合数
 last = ret;
 }
 if(ret != )return true;
 else return false;
 }
 //**************************************************
 // Miller_Rabin算法
 // 是素数返回true,(可能是伪素数)
 // 不是素数返回false
 //**************************************************
 bool Miller_Rabin(long long n)
 {
 if( n < )return false;
 if( n == )return true;
 if( (n&) == )return false;//偶数
 long long x = n - ;
 long long t = ;
 while( (x&)== ){x >>= ; t++;}
 srand(time(NULL));/* *************** */
 for(int i = ;i < S;i++)
 {
 long long a = rand()%(n-) + ;
 if( check(a,n,x,t) )
 return false;
 }
 return true;
 }
 //**********************************************
 // pollard_rho 算法进行质因素分解
 //
 //
 //*********************************************
 long long factor[];//质因素分解结果（刚返回时时无序的）
 int tol;//质因素的个数，编号0~tol-1
 long long gcd(long long a,long long b)
 {
 long long t;
 while(b)
 {
 t = a;
 a = b;
 b = t%b;
 }
 if(a >= )return a;
 else return -a;
 }
 //找出一个因子
 long long pollard_rho(long long x,long long c)
 {
 long long i = , k = ;
 srand(time(NULL));
 long long x0 = rand()%(x-) + ;
 long long y = x0;
 while()
 {
 i ++;
 x0 = (mult_mod(x0,x0,x) + c)%x;
 long long d = gcd(y - x0,x);
 if( d !=  && d != x)return d;
 if(y == x0)return x;
 if(i == k){y = x0; k += k;}
 }
 }
 //对 n进行素因子分解，存入factor. k设置为107左右即可
 void findfac(long long n,int k)
 {
 if(n == )return;
 if(Miller_Rabin(n))
 {
 factor[tol++] = n;
 return;
 }
 long long p = n;
 int c = k;
 while( p >= n)
 p = pollard_rho(p,c--);//值变化，防止死循环k
 findfac(p,k);
 findfac(n/p,k);
 }
 //POJ 1811
 //给出一个N(2 <= N < 2^54),如果是素数，输出"Prime",否则输出最小的素因子
 int main()
 {
     int T;
     long long n;
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     while(scanf("%I64d",&n)==)
     {
         if(n==)
         {
             printf("is not a D_num\n");
             continue;
         }
         tol=;
         findfac(n,);
         if(tol!= && tol!=)
         {
             printf("is not a D_num\n");
             continue;
         }
         sort(factor,factor+tol);
         if(tol==)
         {
             if(factor[]!=factor[])
             {
                 printf("%I64d %I64d %I64d\n",factor[],factor[],factor[]*factor[]);
                 continue;
             }
             else
             {
                 printf("is not a D_num\n");
                 continue;
             }
         }
         if(tol==)
         {
             if(factor[]==factor[]&&factor[]==factor[])
             {
                 printf("%I64d %I64d %I64d\n",factor[],factor[]*factor[],factor[]*factor[]*factor[]);
                 continue;
             }
             else
             {
                 printf("is not a D_num\n");
                 continue;
             }
         }
     }
     return ;
 }