poj2191 pollard-rho大数分解质因子+Miller_Rabin判断质数

/*
Time:2019.12.11
Author: Goven
type:pollard-rho大数分解质因子+Miller_Rabin判断质数
ref:
*/
#include<iostream>
#include<ctime>
#include<cstdlib>
#include<algorithm> 
using namespace std;
typedef long long ll;
 
int cnt = 0;
ll factor[1000];

ll Mult_mod (ll a, ll b, ll mod) {//大数乘法 
    a %= mod;
    b %= mod;
    ll res = 0;
    while (b) {
        if (b & 1) res = (res + a) % mod;
        a = (a << 1) % mod;
        b >>= 1;
    }
    return res;
}

ll pow_mod (ll a, ll b, ll mod) {
    a %= mod;
    ll res = 1;
    while (b) {
        if (b & 1) res = Mult_mod(res, a, mod);
        a = Mult_mod(a, a, mod);
        b >>= 1;
    }
    return res;
}

bool Check (ll a, ll x, ll t, ll mod) {//二次探测判断是否为合数 
    a = pow_mod(a, x, mod); 
    ll last;
    for (int i = 0; i < t; i++) {
        last = a;
        a = Mult_mod(a, a, mod);
        if (a == 1 && last != 1 && last != mod - 1) return true;
    }
    if (a != 1) return true;
    return false;
} 

bool Miller_Rabin (ll n, int k) {//Miller_Rabin判断素数 
    if (n < 2) return false;
    if (n == 2) return true;
    if (n & 1 == 0) return false;//step1:粗筛
    //step2:二次探测
    ll t = 0, x = n - 1;//att1: x这里的赋值别忘了减1 
    while (x & 1 == 0) x >>= 1, t++;  
    for (int i = 0; i < k; i++) {
        ll a = rand() % (n - 1) + 1;
        if (Check(a, x, t, n))
            return false; 
    }
    return true;
}

ll gcd (ll a, ll b) {
    if (a < 0) return gcd(-a, b);
    return b ? gcd(b, a % b) : a;
}

ll f (ll x, ll c, ll mod) {//生成随机数:f(x) = (x * x + c) % mod 
    return (Mult_mod(x, x, mod) + c) % mod; 
}

ll Pollard_rho (ll n, ll c) {//Pollard_rho找到n的某个因子 
    ll a = rand() % n;
    ll b = a;
    ll i = 1, k = 2;
    while (1) {
        i++;
        a = f(a, c, n);
        ll d = gcd(a - b, n);
        if (d != 1 && d != n) return d;
        if (a == b) return n;//到了随机数的循环节
        if (i == k) {
            b = a;
            k += k;
        } 
    }
}
void findFac (ll n) {//求解n的素因子
    if (Miller_Rabin(n, 5)) {
        factor[cnt++] = n;
        return;
    } 
    ll p = n;
    while (p >= n) p = Pollard_rho(p, (ll)rand() % (n - 1) + 1);
    findFac(p);
    findFac(n / p); 
}

bool isPrime (int x) {
    for (int i = 2; i <= x >> 1; i++) {
        if (x % i == 0) return false;
    }
    return true;
}
int main()
{
    srand(time(NULL));
    
    int k;
    cin >> k;
    ll n;
    
    for (int i = 2; i <= k; i++) {
        if (!isPrime(i)) continue;
        n = ((ll)1 << i ) - 1;//err1:1前面要加(ll) 
        cnt = 0;
        findFac(n); 
        if (cnt == 1) continue;
        sort(factor, factor + cnt);
        printf("%lld", factor[0]);
        for (int j = 1; j < cnt; j++) {
            printf(" * %lld", factor[j]);
        }
        printf(" = %lld = ( 2 ^ %d ) - 1\n", n, i);
    }
    
    return 0;
}

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 201,468评论 5 473
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,620评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 148,427评论 0 334
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,160评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,197评论 5 363
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,334评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,775评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,444评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,628评论 1 295
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,459评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,508评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,210评论 3 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,767评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,850评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,076评论 1 258
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,627评论 2 348
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,196评论 2 341

推荐阅读更多精彩内容