题意:

给你一个数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 << 61;
const int maxC = 240; ll big_mul(ll a, ll b, ll n) {
ll ret = 0;
a %= n;
while(b) {
if(b & 1) {
ret += a;
if(ret >= n) ret -= n;
}
a *= 2;
if(a >= n) a -= n;
b /= 2;
}
return ret;
} ll pow_mod(ll x, ll n, ll m) {
ll ret = 1;
x %= n;
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;
}

POJ 1811 Prime Test 素性测试 分解素因子的更多相关文章

  1. Miller_rabin算法+Pollard_rho算法 POJ 1811 Prime Test

    POJ 1811 Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 32534   Accepted: 8 ...

  2. POJ 1811 Prime Test (Rabin-Miller强伪素数测试 和Pollard-rho 因数分解)

    题目链接 Description Given a big integer number, you are required to find out whether it's a prime numbe ...

  3. Miller&&Pollard POJ 1811 Prime Test

    题目传送门 题意:素性测试和大整数分解, N (2 <= N < 254). 分析:没啥好讲的,套个模板,POJ上C++提交 收获:写完这题得到模板 代码: /************** ...

  4. 数论 - Miller_Rabin素数测试 + pollard_rho算法分解质因数 ---- poj 1811 : Prime Test

    Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 29046   Accepted: 7342 Case ...

  5. poj 1811 Prime Test 大数素数测试+大数因子分解

    Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 27129   Accepted: 6713 Case ...

  6. POJ 1811 Prime Test( Pollard-rho整数分解经典题 )

    链接:传送门 题意:输入 n ,判断 n 是否为素数,如果是合数输出 n 的最素因子 思路:Pollard-rho经典题 /************************************** ...

  7. POJ 1811 Prime Test (Pollard rho 大整数分解)

    题意:给出一个N,若N为素数,输出Prime.若为合数,输出最小的素因子.思路:Pollard rho大整数分解,模板题 #include <iostream> #include < ...

  8. POJ 1811 Prime Test(Miller-Rabin & Pollard-rho素数测试)

    Description Given a big integer number, you are required to find out whether it's a prime number. In ...

  9. POJ 1811 Prime Test

    题意:对于一个大整数,判断是否质数,如果不是质数输出最小质因子. 解法:判断质数使用Miller-Rabin测试,分解质因子使用Pollard-Rho,Miller-Rabin测试用的红书模板,将测试 ...

随机推荐

  1. sendto() 向广播地址发包返回errno 13, Permission denied错误

    http://blog.csdn.net/guanghua2_0beta/article/details/52483916 sendto() 向广播地址发包返回errno 13, Permission ...

  2. Swift的基础,操作符,字符串和集合类型

    这篇文章主要讲解苹果Swift官方指南的第二章前四节的要点内容,如果想看完整的英文文档可以去苹果开发者页面下载. Basic 声明常量let 声明变量var 注释依旧使用"//" ...

  3. LoadRunner性能测试中Controller场景创建需注意的几点

    在LR工具做性能测试中,最关键的一步是Controller场景的设计,因为场景的设计与测试用例的设计相关联,而测试用例的执行,直接影响最终的测试结果是怎么的,因此,我们每设计一种场景,就有可能是一个测 ...

  4. centos添加sudoers

    CentOS添加sudoers 1.chmod 740 /etc/sudoers 2. ## Allow root to run any commands anywhere     root    A ...

  5. Foundation--NSArray+NSMutableArray

    C语言中数组只能存放一种类型 OC语言中数组可以存放多种类型,但是只能是对象类型,不可以是基本类型 我们实际中最好在一个数组中只存放一种类型的元素 nil表示数组的结尾,所以我们不能填写值为nil的对 ...

  6. 【转载自友盟消息推送iOS文档】在appDelegate中注册推送

    1.2   基本功能集成指南 提示 请先在友盟的消息推送管理后台中创建App,获得AppKey和AppSecret 导入SDK 下载 UMessage_Sdk_All_x.x.x.zip并解压缩 导入 ...

  7. 控制器的view的加载优先级

    拿到控制器后,控制器的view是在什么时候按照什么优先级创建加载的? 1.控制器内部的view是延迟加载 1> 用到时再加载(loadView) 2> 加载完毕后会调用控制器的viewDi ...

  8. IOS中设置状态栏的状态

    IOS上 关于状态栏的相关设置(UIStatusBar) 知识普及 ios上状态栏 就是指的最上面的20像素高的部分 状态栏分前后两部分,要分清这两个概念,后面会用到: 前景部分:就是指的显示电池.时 ...

  9. js 去掉字符串最后一个字符

    var str = str.substring(0,str.length-1); alert(str); 注:length 需小写.

  10. weblogic启动时报错 java.lang.OutOfMemoryError: PermGen space

    PermGen space的全称是Permanent Generation space,是指内存的永久保存区域.这一部分用于存放Class和Meta的信息,Class在被 Load的时候被放入Perm ...