POJ - 2447 RSA (PollarRho大整数的因数分解)

最新推荐文章于 2024-06-21 08:00:00 发布

hyesuixin

最新推荐文章于 2024-06-21 08:00:00 发布

阅读量927

点赞数

分类专栏：数论文章标签： ACM 算法数论

本文链接：https://blog.csdn.net/hyesuixin/article/details/77372652

版权

数论专栏收录该内容

29 篇文章 0 订阅

订阅专栏

本文介绍了如何运用PollarRho算法解决RSA系统中的大整数因数分解问题。在RSA加密中，给定公钥{E,N}和C，需要找到M。常规方法无法处理大数据的质因子分解，因此引入PollarRho算法。通过该算法找到P和Q，进而计算出T、D，从而求得M。文章提供了AC代码，但PollarRho算法的复杂原理建议参考其他博主的深入讲解。" 113011935,10546102,Vue2与Node.js接口结合WebSocket本地测试教程,"['前端开发', 'WebSocket', 'Vue框架', 'Node.js']

摘要由CSDN通过智能技术生成

题目描述：

点击打开链接

描述一下题目中所给的RSA系统。首先找到两个大整数P,Q，N=PQ,然后令T=(P-1)(Q-1)，令gcd(E,T)==1,再搞一个D，令(E*D)%T=1，很明显这个D是E%T的逆元。得到公匙{E,N}，私匙{D,N},然后进行加密，C=M^E%N,M=C^D%N，现在题目给你C和{E,N}，需要你求M。

思路还是比较明确的，通过N找到P,Q然后求出T，得到E，再求D最后就可以求出M，后面的逆元和gcd都不是麻烦事，麻烦的在于找P,Q，这题N的数据太大，常规的质因子分解无法解决，这里要引入PollarRho算法，这个算法的原理很复杂，可能但讲一篇博客都很难讲清楚，这里推荐另一位博主的博客：点击打开链接。

AC代码：

#include<iostream>
#include<cstdlib>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<string>
#include<stack>
#include<queue>
#include<algorithm>
using namespace std;
const int MOD=21252;


long long multi_mod(long long x,long long y,long long c)
{
    x = (x % c + c) % c;
    y = (y % c + c) % c;

    long long ans=0;
    while(y){
        if(y&1)ans=ans+x;
        if(ans>=c)ans-=c;
        x=x+x;
        if(x>=c)x-=c;
        y>>=1;
    }
    return ans;
}
long long gcd(long long x,long long y)
{
    return y==0? x:gcd(y,x%y);
}
long long qpow(long long x,long long y,long long MOD)
{
    long long res=1;
    while(y)
    {
        if (y&1) res=multi_mod(res,x,MOD);
        x=multi_mod(x,x,MOD);
        y=y>>1;
    }
    return res;
}
long long PollarRho(long long n,long long c)
{
    srand(time(NULL));
    long long x=rand()%(n-1)+1;
    long long y=x;
    int i=1;
    int k=2;
    while(true)
    {
        i++;
        x=(qpow(x,2,n)+c)%n;
        long long d=gcd(n+y-x,n);
        if (1<d&&d<n) return d;
        if (y==x) return n;
        if (i==k)
        {
            y=x;
            k=k*2;
        }
    }
}
long long extgcd(long long a,long long b,long long &x,long long &y)//扩展欧几里德求逆元。
{
    long long d=a;
    if (b!=0)
    {
        d=extgcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else
    {
        x=1;
        y=0;
    }
    return d;
}

long long c,e,n;

int main()
{
    while(cin>>c>>e>>n)
    {
        long long p=PollarRho(n,10007);
       // cout<<p<<endl;
        long long q=n/p;
        long long t=(p-1)*(q-1);
        long long x0,y0;
        long long d=extgcd(e,t,x0,y0);
        x0=(x0%t+t)%t;
        cout<<qpow(c,x0,n)<<endl;
    }
    return 0;
}