快判素数

#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long ll;
ll ran[3]={2,7,61};
ll powl(ll a,ll b,ll mod)
{
    ll res=1,sum=a%mod;
    while(b)
    {
        if(b&1)
            res=(res*sum)%mod;
        b>>=1;
        sum=sum*sum%mod;
    }
    return res;
}//快速幂,用来加快计算和防止爆long long 
bool witness(ll a,ll n)//Miller Rabin算法的重点 
{
    ll t=n-1;
    int j=0;
    while(t%2==0)
    {
        t/=2;
        j++;
    }//把n-1拆分为2^s*r的形式,这里j为s,tem为r 
    ll x=powl(a,t,n);
    if(x==1||x==n-1) return true;
    while(j--)
    {
        x=x*x%n;
        if(x==n-1) return true;
    }//否则判断等式a^(2jr) ≡-1 mod n 看是否有满足的 j  
    return false;
}
bool Miller_Rabin(ll n)
{
    if(n==2||n==7||n==61) return true;
    if(n<2||n%2==0||n%7==0||n%61==0) return false;//因为要满足a与n是互素的任何随机整数,所以要特判 
    for(int i=0;i<3;i++)
    {
        ll a=ran[i];
        if(a!=n)
            if(!witness(a,n)) return false;            
    }
    return true;
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    while(m--)
    {
        scanf("%d",&n);
        if(Miller_Rabin(n))printf("Yes\n");
        else printf("No\n");
    }
    return 0;
} 

猜你喜欢

转载自blog.csdn.net/qq_42756710/article/details/81269894