0
点赞
收藏
分享

微信扫一扫

米勒拉宾素数测试

流沙雨帘 2022-08-31 阅读 54


ll quick_pow(ll a,ll b,ll r)//快速幂 
{
ll ans = 1;
while(b)
{
if(b&1)ans=(ans*a)%r;
a = (a*a)%r;
b>>=1;
}
return ans;
}
bool Miller_Rabbin(int n,int a)//米勒拉宾素数测试
{
int r=0,s=n-1,j;
if(!(n%a))
return false;
while(!(s&1)){
s>>=1;
r++;
}
ll k=qpow(a,s,n);
if(k==1)
return true;
for(j=0;j<r;j++,k=k*k%n)
if(k==n-1)
return true;
return false;
}
bool IsPrime(int n)//判断是否是素数
{
int tab[]={2,3,5,7};
for(int i=0;i<4;i++)
{
if(n==tab[i])
return true;
if(!Miller_Rabbin(n,tab[i]))
return false;
}
return true;
}


举报

相关推荐

0 条评论