POJ 1811 Prime Test 素性测试 分解素因子

题意:

给你一个数n(n <= 2^54),判断n是不是素数,如果是输出Prime,否则输出n最小的素因子


解题思路:

自然数素性测试可以看看Matrix67的  素数与素性测试

素因子分解利用的是Pollard rho因数分解,可以参考 Pollard rho因数分解

存个代码~


/* **********************************************
Author      : JayYe
Created Time: 2013-9-25 16:02:25
File Name   : JayYe.cpp
*********************************************** */

#include <stdio.h>
#include <string.h>
#include <time.h>
#include <algorithm>
using namespace std;
#define Time 12 // Miller测试次数
typedef __int64 ll;

const ll INF = (1LL << 62) + ((1LL<<62)-1);
const int maxC = 240;

ll big_mul(ll a, ll b, ll m) {
    ll ret = 0;
    a %= m;
    while(b) {
        if(b & 1) {
            ret += a;
            if(ret >= m)    ret -= m;
        }
        a *= 2;
        if(a >= m)  a -= m;
        b /= 2;
    }
    return ret;
}

ll pow_mod(ll x, ll n, ll m) {
    ll ret = 1;
    x %= m;
    while(n) {
        if(n & 1)   ret = big_mul(ret, x, m);
        x = big_mul(x, x, m);
        n /= 2;
    }
    return ret;
}

// 以a为基对n进行Miller次测试并进行二次探测,返回true则是合数
bool Wintess(ll a, ll n) {
    ll m = n-1;
    int top = 0;
    // n-1 = m*(2^top)
    while(m % 2 == 0) {
        m /= 2;
        top++;
    }
    ll x = pow_mod(a, m, n), y;
    for(int i = 0;i < top; i++) {
        y = big_mul(x, x, n);
        if(y == 1 && (x != 1 && x != n-1))
            return true;
        x = y;
    }
    if(y > 1)   return true;
    return false;
}

// 对n进行ts次 Miller素性测试
bool Miller_Rabin(int ts, ll n) {
    if(n == 2)  return true;
    if(n == 1 || n % 2 == 0)  return false;
    srand(time(NULL));
    for(int i = 0;i < ts; i++) {
        ll a = rand() % (n-1) + 1;
        if(Wintess(a, n))   return false;
    }
    return true;
}

ll ans;

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

// 对n进行因式分解,找出n的一个因子,该因子不一定是最小的
ll Pollard(ll n, int c) {
    srand(time(NULL));
    ll i = 1, k = 2, x = rand()%n, y = x;
    while(true) {
        i++;
        x = (big_mul(x, x, n) + c) % n;
        ll d = gcd(y - x, n);
        if(d > 1 && d < n)  return d;
        if(y == x)  return n; // 如果该数已经出现过,直接返回
        if(i == k) {
            y = x; k <<= 1;
        }
    }
}

// 找出所有素因子
void solve(ll n, int c) {
    if(n == 1)  return ;
    // 判断是否为素数
    if(Miller_Rabin(Time, n)) {
        if(ans > n) ans = n;
        return ;
    }
    ll m = n;
    while(m == n) {   // 找出n的一个因子
        m = Pollard(n, c--);
    }
    solve(m, c);
    solve(n/m, c);
}

int main() {
    int t;
    ll n;
    scanf("%d", &t);
    while(t--) {
        scanf("%I64d", &n);
        if(Miller_Rabin(Time, n))
            puts("Prime");
        else {
            ans = INF;
            solve(n, maxC);
            printf("%I64d\n", ans);
        }
    }
    return 0;
}

 


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值