MillerRabin与PollardRho

最新推荐文章于 2023-06-09 00:03:03 发布
最新推荐文章于 2023-06-09 00:03:03 发布
阅读量406 收藏
点赞数
分类专栏： # 系列 文章标签： 素数检测
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/71112434
版权

BZOJ 4802 欧拉函数

Description

已知N，求phi(N)
Input

正整数N。N<=10^18
Output

输出phi(N)
Sample Input

8
Sample Output

4

直接套板子

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<vector>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define vll vector<ll> 
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10LL+ch-'0';ch=getchar();}
    return x*f;
}
inline ll gcd(ll a,ll b){if (!b) return a;return gcd(b,a%b);} 
inline ll abs2(ll x){if (x>=0) return x;return -x;} 
inline ll add(ll a,ll b,ll p){a+=b; return a>=p ? a - p : a;}
inline ll mul(ll a,ll b,ll p){
    ll rst = 0;
    while (b) {
        if (b&1) rst = add(rst, a, p);
        b >>= 1, a = add(a, a, p);
    }
    return rst;
}
inline ll pow2(ll a,ll b,ll p)  //a^b mod p 
{  
    ll rst = 1;
    while (b) {
        if (b&1) rst = mul(rst, a, p);
        b >>= 1, a = mul(a, a, p);
    }
    return rst;
}
inline bool witness2(ll a, ll x) {
    ll m = x - 1, t = 0;
    while (! (m&1)) {
        m>>=1, ++t;
    }
    ll fac = pow2(a, m, x), lst = fac;
    while (t--) {
        fac = mul(fac, fac, x);
        if (fac == 1 and lst != 1 and lst != x - 1) return true;
        lst = fac;
    } return lst != 1;
}
#define CAP (20) 
inline ll range(ll l,ll r){return l+rand()%(r-l+1); }
inline bool MillerRabin(ll x) {
    if (x < 2LL) return false;
    if (x == 2LL) return true;
    if (! (x&1LL)) return false;
    for (int i = 1; i <= CAP; i++) {
        ll a = range(2, x-1);
        if (witness2(a, x)) return false;
    } return true;
}
inline ll PollardRho(ll x, ll a) {
    ll x1 = range(0, x - 1), x2 = x1;
    while(1) {
        x2 = add(mul(x2, x2, x), a, x);
        x2 = add(mul(x2, x2, x), a, x);
        ll d = gcd(abs2(x2 - x1), x);
        if (d != 1 and d != x) return d;
        x1 = add(mul(x1, x1, x), a, x);
        if (x1 == x2) return x;
    } return x; 
}

vector<ll> ls;
inline void findFac(ll x) {
    if (MillerRabin(x))
        return ls.push_back(x), void(0);
    ll p = x;
    while (p >= x) p = PollardRho(p, range(0, x - 1));
    findFac(p), findFac(x / p  );
}

int main()
{
//  freopen("bzoj4802.in","r",stdin);
//  freopen(".out","w",stdout);

    srand(23333);
    ll n=read();
    if (n==1) puts("1");
    else if (MillerRabin(n)) cout<<n-1<<endl;
    else {
        findFac(n);
        sort(ALL(ls));
        ls.erase(unique(ALL(ls)),ls.end());
        ll phi=n;
        vll::iterator it;
        for(it=ls.begin();it!=ls.end();it++) {
            phi=phi/(*it)*((*it)-1);
        }
        cout<<phi<<endl;
    }

    return 0;
}

HDU 1164 Eddy’s research I

Problem Description

Eddy’s interest is very extensive, recently he is interested in prime number. Eddy discover the all number owned can be divided into the multiply of prime number, but he can’t write program, so Eddy has to ask intelligent you to help him, he asks you to write a program which can do the number to divided into the multiply of prime number factor .

Input

The input will contain a number 1 < x<= 65535 per line representing the number of elements of the set.

Output

You have to print a line in the output for each entry with the answer to the previous question.

Sample Input

11

9412

Sample Output

11

2*2*13*181

Author

eddy

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<vector>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define vll vector<ll> 
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10LL+ch-'0';ch=getchar();}
    return x*f;
}
inline ll gcd(ll a,ll b){if (!b) return a;return gcd(b,a%b);} 
inline ll abs2(ll x){if (x>=0) return x;return -x;} 
inline ll add(ll a,ll b,ll p){a+=b; return a>=p ? a - p : a;}
inline ll mul(ll a,ll b,ll p){
    ll rst = 0;
    while (b) {
        if (b&1) rst = add(rst, a, p);
        b >>= 1, a = add(a, a, p);
    }
    return rst;
}
inline ll pow2(ll a,ll b,ll p)  //a^b mod p 
{  
    ll rst = 1;
    while (b) {
        if (b&1) rst = mul(rst, a, p);
        b >>= 1, a = mul(a, a, p);
    }
    return rst;
}
inline bool witness2(ll a, ll x) {
    ll m = x - 1, t = 0;
    while (! (m&1)) {
        m>>=1, ++t;
    }
    ll fac = pow2(a, m, x), lst = fac;
    while (t--) {
        fac = mul(fac, fac, x);
        if (fac == 1 and lst != 1 and lst != x - 1) return true;
        lst = fac;
    } return lst != 1;
}
#define CAP (20) 
inline ll range(ll l,ll r){return l+rand()%(r-l+1); }
inline bool MillerRabin(ll x) {
    if (x < 2LL) return false;
    if (x == 2LL) return true;
    if (! (x&1LL)) return false;
    for (int i = 1; i <= CAP; i++) {
        ll a = range(2, x-1);
        if (witness2(a, x)) return false;
    } return true;
}

inline ll PollardRho(ll x, ll a) {
    ll x1 = range(0, x - 1), x2 = x1;
    while(1) {
        x2 = add(mul(x2, x2, x), a, x);
        x2 = add(mul(x2, x2, x), a, x);
        ll d = gcd(abs2(x2 - x1), x);
        if (d != 1 and d != x) return d;
        x1 = add(mul(x1, x1, x), a, x);
        if (x1 == x2) return x;
    } return x; 
}
// need clear ls
// x>1 
vector<ll> ls;
inline void findFac(ll x) {
    if (MillerRabin(x))
        return ls.push_back(x), void(0);
    ll p = x;
    while (p >= x) p = PollardRho(p, range(0, x - 1));
    findFac(p), findFac(x / p  );
}

int main()
{
//  freopen("hdu1164.in","r",stdin);
//  freopen(".out","w",stdout);

    srand(time(NULL));
    int n;
    while(scanf("%d",&n)!=EOF) {
        int ans=0;
        ls.erase(ALL(ls));
        findFac(n);
        sort(ALL(ls));
        Rep(i,SI(ls)) {
            if (i) putchar('*');
            printf("%lld",ls[i]);
        }
        puts("");
    }

    return 0;
}

HDU D_num

一个数是否有且只有4个因子。

分解质因数，2种情况

  • 1,a,b,a∗b
  • 1,a,a2,a3 
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<vector>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define vll vector<ll> 
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10LL+ch-'0';ch=getchar();}
    return x*f;
}
inline ll gcd(ll a,ll b){if (!b) return a;return gcd(b,a%b);} 
inline ll abs2(ll x){if (x>=0) return x;return -x;} 
inline ll add(ll a,ll b,ll p){a+=b; return a>=p ? a - p : a;}
inline ll mul(ll a,ll b,ll p){
    ll rst = 0;
    while (b) {
        if (b&1) rst = add(rst, a, p);
        b >>= 1, a = add(a, a, p);
    }
    return rst;
}
inline ll pow2(ll a,ll b,ll p)  //a^b mod p 
{  
    ll rst = 1;
    while (b) {
        if (b&1) rst = mul(rst, a, p);
        b >>= 1, a = mul(a, a, p);
    }
    return rst;
}
inline bool witness2(ll a, ll x) {
    ll m = x - 1, t = 0;
    while (! (m&1)) {
        m>>=1, ++t;
    }
    ll fac = pow2(a, m, x), lst = fac;
    while (t--) {
        fac = mul(fac, fac, x);
        if (fac == 1 and lst != 1 and lst != x - 1) return true;
        lst = fac;
    } return lst != 1;
}
#define CAP (20) 
inline ll range(ll l,ll r){return l+rand()%(r-l+1); }
inline bool MillerRabin(ll x) {
    if (x < 2LL) return false;
    if (x == 2LL) return true;
    if (! (x&1LL)) return false;
    for (int i = 1; i <= CAP; i++) {
        ll a = range(2, x-1);
        if (witness2(a, x)) return false;
    } return true;
}

inline ll PollardRho(ll x, ll a) {
    ll x1 = range(0, x - 1), x2 = x1;
    while(1) {
        x2 = add(mul(x2, x2, x), a, x);
        x2 = add(mul(x2, x2, x), a, x);
        ll d = gcd(abs2(x2 - x1), x);
        if (d != 1 and d != x) return d;
        x1 = add(mul(x1, x1, x), a, x);
        if (x1 == x2) return x;
    } return x; 
}
// need clear ls
// x>1 
vector<ll> ls;
inline void findFac(ll x) {
    if (MillerRabin(x))
        return ls.push_back(x), void(0);
    ll p = x;
    while (p >= x) p = PollardRho(p, range(0, x - 1));
    findFac(p), findFac(x / p  );
}

int main()
{
//  freopen("hdu3864.in","r",stdin);
//  freopen(".out","w",stdout);

    srand(time(NULL));
    ll n;
    while(scanf("%lld",&n)!=EOF) {
        if (n==1) {
            puts("is not a D_num");
            continue;
        }
        int ans=0;
        ls.erase(ALL(ls));
        findFac(n);
        sort(ALL(ls));

        if(SI(ls)==2) {
            if (ls[0]!=ls[1]) printf("%lld %lld %lld\n",ls[0],ls[1],n);
            else puts("is not a D_num");
        }
        else if (SI(ls)==3) {
            if (ls[0]==ls[2]) printf("%lld %lld %lld\n",ls[0],ls[0]*ls[0],n);
            else puts("is not a D_num");
        } 
        else puts("is not a D_num");



    }

    return 0;
}

HDU 4910 Problem about GCD

题意：给定一个数n，求所有小于等于n且与n互素的数的乘积再mod n

打表，发现
f(n)={n−1,1,n=1,2,4 or n=pk,2∗pk,p is prime and p≠2 else 

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<vector>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define vll vector<ll> 
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10LL+ch-'0';ch=getchar();}
    return x*f;
}
inline ll gcd(ll a,ll b){if (!b) return a;return gcd(b,a%b);} 
inline ll abs2(ll x){if (x>=0) return x;return -x;} 
inline ll add(ll a,ll b,ll p){a+=b; return a>=p ? a - p : a;}
inline ll mul(ll a,ll b,ll p){
    ll rst = 0;
    while (b) {
        if (b&1) rst = add(rst, a, p);
        b >>= 1, a = add(a, a, p);
    }
    return rst;
}
inline ll pow2(ll a,ll b,ll p)  //a^b mod p 
{  
    ll rst = 1;
    while (b) {
        if (b&1) rst = mul(rst, a, p);
        b >>= 1, a = mul(a, a, p);
    }
    return rst;
}
inline bool witness2(ll a, ll x) {
    ll m = x - 1, t = 0;
    while (! (m&1)) {
        m>>=1, ++t;
    }
    ll fac = pow2(a, m, x), lst = fac;
    while (t--) {
        fac = mul(fac, fac, x);
        if (fac == 1 and lst != 1 and lst != x - 1) return true;
        lst = fac;
    } return lst != 1;
}
#define CAP (20) 
inline ll range(ll l,ll r){return l+rand()%(r-l+1); }
inline bool MillerRabin(ll x) {
    if (x < 2LL) return false;
    if (x == 2LL) return true;
    if (! (x&1LL)) return false;
    for (int i = 1; i <= CAP; i++) {
        ll a = range(2, x-1);
        if (witness2(a, x)) return false;
    } return true;
}

inline ll PollardRho(ll x, ll a) {
    ll x1 = range(0, x - 1), x2 = x1;
    while(1) {
        x2 = add(mul(x2, x2, x), a, x);
        x2 = add(mul(x2, x2, x), a, x);
        ll d = gcd(abs2(x2 - x1), x);
        if (d != 1 and d != x) return d;
        x1 = add(mul(x1, x1, x), a, x);
        if (x1 == x2) return x;
    } return x; 
}
// need clear ls
// x>1 
vector<ll> ls;
inline void findFac(ll x) {
    if (MillerRabin(x))
        return ls.push_back(x), void(0);
    ll p = x;
    while (p >= x) p = PollardRho(p, range(0, x - 1));
    findFac(p), findFac(x / p  );
}

int main()
{
//  freopen("hdu3864.in","r",stdin);
//  freopen(".out","w",stdout);

    srand(time(NULL));
//  For(i,100) {
//      int p=1;
//      For(j,i) if (gcd(i,j)==1) p=p*j%i;
//      cout<<i<<':'<<p<<endl;
//  }
    ll n;
    while(scanf("%lld",&n)!=EOF&& n!=-1) {
        ll ans;
        if (n==1||n==2||n==4) ans=n-1;
        else {
            ls.erase(ALL(ls));
            findFac(n);
            sort(ALL(ls));

            int sz=SI(ls);
            if (ls[0]==ls[sz-1]) {
                ans=(ls[0]==2) ?1: (n-1);
            } else {
                ans=1;
                if (SI(ls)>=2) {
                    if (ls[0]==2&&ls[0]!=ls[1]&&ls[1]==ls[sz-1]) ans=n-1;
                }
            }

        }
        cout<<ans<<endl;
    }
    return 0;
}
nike0good
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【模板】Pollard-Rho算法
C202044zxy的博客
09-18 675
一、题目 电磁看题 二、解法 作者讲的解法只是一个缩略版，要想获得更好的理解请参照大佬的博客。 0x01 生日悖论 如果你不是很了解，点此观看 生日悖论告诉了我们一个道理，也就是我们用组合计算的方式会比单个确切值更有概率得到答案。什么意思呢？举个例子，我们要找到xxx，用r−lr-lr−l的方式比直接randrandrand概率有更大概率得到xxx，这里作者就不给详细证明了，感性理解。 所以与其直...
Pollard Rho算法
思维缜密的博客
01-17 1189
假设一年有n天，房间中有k人，每个人的生日在这n天中，服从均匀分布，两个人的生日相互独立问至少要有多少人，才能使其中两个人生日相同的概率达到p解：考虑k≤n设k个人生日互不相同为事件A，则PAnn​nn−1​⋯nn−k1​由题意PA≤1−p,由1x≤exPA≤e−n1​e−n2​⋯enk−1​e−2nkk−1​≤1−p解得k≥211−2nln。
Pollard Rho因子分解算法
weixin_30335353的博客
02-09 310
　　有一类问题，要求我们将一个正整数x，分解为两个非平凡因子（平凡因子为1与x）的乘积x=ab。 　　显然我们需要先检测x是否为素数（如果是素数将无解），可以使用Miller-Rabin算法来进行测试。 　　Pollard Rho是一个非常玄学的方式，用于在O(n^1/4)的期望时间复杂度内计算合数n的某个非平凡因子。事实上算法导论给出的是O(√p)，p是n的某个最小因子，满足p与n/p互质。...
pollard_rho学习笔记 数论
forever_shi的博客
06-02 1515
pollard_rho算法是用来解决质因数分解问题的。 我们知道，朴素的质因数分解是O(n−−√)O(n)O(\sqrt{n})的，但是如果nnn很大的话应该怎么办呢？这时候就要用到pollard_rho了。 pollard_rho也是一种基于随机的算法，它的思路是先用miller_rabin来判断当前数是否已经是素数了，如果是的话记录并返回。如果不是，我们设要分解的数为nnn,那么我们考虑去...
Pollard’s rho算法
永远相信美好的事情即将发生
06-09 538
给定一正整数N∈N∗，试快速找到它的一个因数。很久很久以前，我们曾学过试除法来解决这个问题。很容易想到因数是成对称分布的：即N1N​和N​N。这个很容易想清楚，我们只需要把区间1N​扫一遍就可以了，此时试除法的时间复杂度为ON​。对于N≥1018的数据，这个算法运行起来无疑是非常糟糕的。我们希望有更快的算法。对于这样的算法，一个比较好的想法是去设计一个随机程序，随便猜一个因数。如果你运气好，这个程序的时间复杂度下限为O1。但对于一个N≥10。
MillerRabin素性检验
01-21
nt MillerRabin(unsigned char Odd[],unsigned char Random[],int s,unsigned char Mp[]){//MillerRabin素性检验//OK unsigned char A[M],B[M]; int t,i; PowerFunction(Random,Mp,Odd,A); t=Judgment1(A); if...
MillerRabin.java
最新发布
11-27
MillerRabin素数检测算法
rsa.rar_mpz_millerrabin_rsa_脉动阵列
09-24
4. 效率与安全性：在优化速度的同时，不能忽视算法的安全性，确保在硬件实现中仍能抵御已知的攻击。 总结来说，"rsa.rar_mpz_millerrabin_rsa_脉动阵列"涉及的是使用GMP库中的`mpz_millerrabin`素数测试，并结合...
MillerRabin素数测试
10-29
三、算法优化与错误率 Miller-Rabin测试虽然有可能产生错误，但错误率可以通过增加测试次数来降低。对于随机选择的a，测试错误的概率是四分之一。通过增加测试次数，可以将总体错误率控制在一个非常小的值，如2^(-t...
素数判定MillerRabin算法
01-01
强素数测试与随机性** Miller-Rabin素数测试是一种基于概率的强素数测试方法，它利用了费马小定理和二次剩余的概念。这个算法不是绝对确定性的，而是具有一定的错误率，即存在一定的概率将合数误判为素数。然而，...
Pollard-Rho算法详解
09-12
Pollard-Rho算法详解
64位以内Rabin-Miller 强伪素数测试和Pollard rho 因数分解算法的实现
10-22
64位以内Rabin-Miller 强伪素数测试和Pollard rho 因数分解算法的实现的C代码
Pollard Rho算法思想
07-28
Pollard Rho算法思想
大整数分解——Pollard Rho算法
xffyjq的博客
09-28 8876
延续上一篇，这次来讲一讲大整数分解算法的应用。要解决的问题很简单，对一个整数进行分解质因数。首先还是效率非常低的暴力算法，相信大家都会，不多提。和上次一样，当数达到非常大的时候，分解将变得非常困难。于是这次又带来一个提升分解速度的“非完美”算法。之所以打引号，是因为这次不完美的不是结果，而是时间效率。Pollard Rho算法分解一个数n的过程大体上是这样子的：1、找到一个数p，使得p|np|n，将
Pollard-Rho算法--大数分解
热门推荐
萌新也要搞算法
04-11 1万+
其实刚开始知道这个算法时，我以为需要字符串操作什么的，毕竟是大数嘛，可这家伙只用了个long long，无语了....long long int，能叫大数吗，连2^100次方都处理不了。说是这样说，这个算法已经不错了，复杂度为O(n^1/4), 貌似目前学界没有找到特别好的算法，据说有什么艾-鲁法，威廉斯夫法，可以分解一个千位素数（需要一周时间）（上面两种算法百度不到QAQ），量子的那个shore...
Pollard-Rho算法证明
总理同学的编程尝试
12-12 1528
我觉得我写得挺模糊的，里面由一些没给出证明的东西，我不会证，欢迎同学们联系我。
pollard-rho&miller-rabin
DJ的博客
02-27 284
今天学习pollard-rho和miller-rabin。 前置： 不管怎么说大力筛一遍2,3,5,7,11,13,17都是错不了的！ Miller-rabin %%% Doggu太神啦！ 前置技能：费马小定理 ap−1≡1(Mod&nbsp;p)p为素数a^{p-1}{\equiv} 1(Mod\ p)p为素数ap−1≡1(Mod&nbsp;p)p为素数 证明： 若gcd(i,j)==1gcd(...
Pollard_rho质因数分解
小蒟蒻wddwjlss的博客
06-04 227
转载自forever_shi神犇点击打开链接
Pollard Rho 大整数分解
I wonder how , I wonder why !!
08-21 1494
http://www.cnblogs.com/jackiesteed/articles/2019910.html 1975年，John M. Pollard提出了第二种因数分解的方法。Pollard rho因数分解方法基于下列几点： (1) 假定有两个整数 和 使得p可以整除-，但是n不能整除 。 (2) 可以证明 。因为p可以整除- ，可以写成 。但是，因为n不能整除-，很
“python”millerrabin素数测试算法
11-25
Miller-Rabin素数测试算法是一种概率性的素数测试算法，用于判断一个给定的大数是否是素数。该算法的原理是基于费马小定理和二次探测定理，通过随机选择一系列的底数进行检测，如果通过所有的底数测试，则该数被认为是素数，否则被认为是合数。 具体的算法步骤如下： 1. 输入一个待判断的大数n。 2. 当n小于等于3时，直接判断为素数，返回True。 3. 将n-1表示为2^s*d的形式，其中d是一个奇数。 - 首先对n-1进行除2操作，直到得到一个奇数d，同时累计计算除2的次数s。 4. 对于给定的底数a，计算b = a^d mod n。 5. 如果b等于1或b等于n-1，则判断为素数，返回True。 6. 重复进行k-1次以下步骤： - 将b的值平方，即b = b^2 mod n。 - 如果b等于n-1，则判断为素数，返回True。 7. 判断为合数，返回False。 在实际运用中，可以通过多次测试底数的方式来提高测试的准确性，每次选择的底数越多，结果越准确。在大多数情况下，该算法可以快速判断一个数是不是合数，但也存在一定的概率误判的风险。 总之，Miller-Rabin素数测试算法是一种基于概率的素数测试方法，能够快速判断一个给定的大数是否为素数。通过随机选择底数并进行测试，可以提高判断的准确性。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值