luogu【模板】线性筛素数 (Miller-Rabin素数测试模板）

这个感觉还是挺好理解的，就是复杂度证明看不懂~

Code: 

#include <cstdio> 
#include <algorithm> 
#include <cstring> 
#include <cmath>   
#define ll long long 
#define setIO(s) freopen(s".in","r",stdin)       
int array[10]={2,3,5,7,11,13,17,23};  
using namespace std;       
ll qpow(ll base,ll k,ll mod) 
{
	ll tmp=1; 
	for(;k;base=base*base%mod,k>>=1) if(k&1) tmp=tmp*base%mod; 
	return tmp; 
}
int prime(ll x)  
{     
	if(x<=1) return 0; 
	int i,j,k; 
	ll pre,a,cur;  
	for(i=0;i<8;++i) 
	{
		if(x==array[i]) return 1;          
		for(cur=x-1,k=0;cur%2==0;cur>>=1) ++k;              
		pre=a=qpow(array[i],cur,x);              
		for(j=1;j<=k;++j) 
		{  
			a=(a*a)%x;                 
			if(a==1&&pre!=1&&pre!=x-1) return 0; 
			pre=a;   
		}
		if(a!=1) return 0;    
	}   
	return 1;   
}
int main() 
{ 
	// setIO("input"); 
	int n,m; 
	scanf("%d%d",&n,&m); 
	for(int i=1;i<=m;++i) 
	{
		ll a; 
		scanf("%lld",&a),printf("%s
",prime(a)?"Yes":"No");  
	}
	return 0; 
}