#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
今日推荐
周排行