POJ 1811 Prime Test（ Pollard-rho整数分解经典题 ）

---

链接：传送门

题意：输入 n ，判断 n 是否为素数，如果是合数输出 n 的最素因子

思路：Pollard-rho经典题

---

/*************************************************************************
    > File Name: Pollard_rho_Test.cpp
    > Author:    WArobot 
    > Blog:      http://www.cnblogs.com/WArobot/ 
    > Created Time: 2017年05月24日 星期三 10时55分04秒
 ************************************************************************/

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;

#define ll long long
#define TIME 30

const int MAX_N = 1000;
const int C = 201;
ll number = 0 , MIN ;
int List[MAX_N];

// --------------------Miller_Rabin素数测试-------------------------

ll RANDOM(ll n){ return (ll)((double)rand()/RAND_MAX * n + 0.5); }

ll quick_multi(ll a,ll b,ll Mod){   // 快速计算 a * b % Mod
    ll ret = 0;
    while(b){
        if(b&1) ret = (ret + a) % Mod;
        a = ( a << 1 ) % Mod;
        b >>= 1;
    }
    return ret % Mod;
}

ll quick_pow(ll a,ll x,ll Mod){     // 快速计算 a^x % Mod
    ll ret = 1;
    while(x){
        if(x&1) { ret = quick_multi( ret , a , Mod); x--; }
        x >>= 1;
        a = quick_multi( a , a , Mod );
    }
    return ret % Mod;
}

bool Witness(ll a,ll n){    // Miller_Rabin 二次探测
    ll tmp = n - 1;
    int j = 0;
    while( !( tmp & 1 )  ){ tmp >>= 1; j++; }
    ll x = quick_pow( a , tmp , n );
    if( x == 1 || x == n-1 )    return false;
    while( j-- ){
        x = x*x % n ;
        if( x == n-1 )  return false;
    }
    return true;
}

bool Miller_Rabin(ll n){    // 检测n是否为素数
    if( n < 2 )     return false;
    if( n == 2 )    return true;
    if( !(n&1) )    return false;
    for(int i = 1 ; i <= TIME ; i++){
        ll a = RANDOM( n - 2 ) + 1; // 得到随机底数
        if( Witness( a , n ) )      return false;
    }
    return true;
}

// --------------------Pollard_rho大整数分解-------------------------

ll gcd(ll a,ll b){  return b == 0 ? a : gcd( b , a%b ); }

ll Pollard_rho(ll n,ll c){  // 找到n的一个因子
    ll x , y , d , i = 1 , k = 2;
    x = RANDOM( n - 1 ) + 1;
    y = x;
    while(true){
        i++;
        x = ( quick_multi( x , x , n ) + c ) % n;
        d = gcd( y - x , n );
        if( 1 < d && d < n )    return d;
        if( y == x )            return n;
        if( i == k ){
            y = x;  k <<= 1;
        }
    }
}

void find(ll n,ll c){
    if( n == 1 )    return;
    if( Miller_Rabin(n) ){
        List[ number++ ] = n ;  
        MIN = min( MIN , n );
        return;
    }
    ll p = n ;
    while( p >= n ) p = Pollard_rho( p , c-- );
    find( p , c );
    find( n/p , c );
}

int main(){
    int T;  ll tar , tmp = pow(2,54)-1;
    scanf("%d",&T);
    while(T--){
        scanf("%lld",&tar);
        if( Miller_Rabin(tar) )     printf("Prime\n");
        else{
            MIN = tmp;
            find( tar , C );
            printf("%lld\n",MIN);
        }
    }
    return 0;
}

转载于:https://www.cnblogs.com/WArobot/p/6898653.html