演算法筆記--米勒-羅賓素數測試

標籤：快速   演算法   prime   als   pac   isp   for   class   while   

由於偽素數的存在，費馬素數測試有了極大的缺陷，於是有了米勒-拉賓素性測試。

模板：

#include <iostream>#include <cmath>using namespace std;long long qpow(int a,int b,int r)//快速冪{    long long ans=1,buff=a;    while(b)    {        if(b&1)ans=(ans*buff)%r;        buff=(buff*buff)%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++;    }    long long 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;}int main(){    long long n;    while(1)    {       cin >> n;    cout << IsPrime(n)<< endl;    }    return 0;}

 

演算法筆記--米勒-羅賓素數測試