标签：const int res ll long Miller 罗宾 mod Rabin

参考博客：博客1  博客2

先贴模板

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define IO ios::sync_with_stdio(false)
#define pb push_back
#define mk make_pair
const int N = 1e5+10;
const int mod = 1e9+7;

ll a, b;

const long long S=20;

ll mult_mod(ll a,ll b,ll mod)//快速乘
{
	ll res=0;
	for(;b;b>>=1){
		if(b&1) res=(res+a)%mod;
		a=(a+a)%mod;
	}
	return res;
}
ll pow_mod(ll a,ll b,ll mod)//快速幂
{
	ll res=1;
	a%=mod;
	for(;b;b>>=1){
		if(b&1) res=mult_mod(res,a,mod);
		a=mult_mod(a,a,mod);
	}
	return res;
}
int check(ll a,ll n,ll x,ll t){
    ll ret=pow_mod(a,x,n);//费马小定理 a^(p-1)%p==1
    ll last=ret;
    for(ll i=1;i<=t;i++){//二次检测定理 如果p是一个素数，则x^2%p==1的解为，则x=1或者x=n-1。
        ret=mult_mod(ret,ret,n);
        if(ret==1&&last!=1&&last!=n-1) return 1;
        last=ret;
    }
    if(ret!=1) return 1;
    return 0;
}
int Miller_Rabin(ll n){
    if(n<2)return 0;
    if(n==2)return 1;
    if((n&1)==0) return 0;
    ll x=n-1;
    ll t=0;
    while((x&1)==0){x>>=1;t++;}
    for(ll i=0;i<S;i++){
        ll a=rand()%(n-1)+1;
        if(check(a,n,x,t))  return 0;
    }
    return 1;
}

int main(){
	for(ll i=100000000000;i<=100000000100;++i){
        if(Miller_Rabin(i)) printf("%lld 是素数\n",i);
        else printf("%lld 不是素数\n",i);
	}
}

根据概率论学 当一件事发生的概率 低于0.01% 时  这件事就不可能发生，即使在竞赛中倒了 ACM历史以来最大的霉运，出现了判断失误的情况 你就再交一次，总不会次次出现低于0.01%的事件吧，比买彩票中大奖、被闪电劈中的概率低多了。

标签：const,int,res,ll,long,Miller,罗宾,mod,Rabin
来源： https://blog.csdn.net/qq_41286356/article/details/106825009

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。