20191016 专题：质数

总览：

判断+筛法

筛法：

埃氏筛法：

时间复杂度（$O(nlogn)$）
模板：

int s[N],num;
bool ex[N];

void prepare(){
	ex[1]=1;
	for(int i=1;i<=n;i++){
		if(ex[i])	continue;
		s[++num]=i;
		for(int j=2;i*j<=n;j++)
			ex[i*j]=1;
	}
}

线性筛法：

时间复杂度（$O(n)$）
模板：

int s[N],num;
bool ex[N];

void prepare(){
	ex[1]=1;
	for(int i==1;i<=n;i++){
		if(ex[i])	continue;
		s[++num]=i;
		for(int j=2;j<=num&&i*s[j]<=n;j++){
			ex[i*s[j]]=1;
			if(!(i%s[j]))	continue;
		}
	}
}

判断：

常规：

时间复杂度：（$O（\sqrt{n}）$）

bool pd(int s){
	for(int i=2;i<=sqrt(s);i++)
		if(!(s%i))	return false;
	return true;
}

miller_rabinrabin:

时间复杂度：（$O(mlogn)$）
错误概率：($4^{-s}$ )
注：选择$2，3，5，7，11，13，17，19$做判断底数，$10^{15}$内保证正确性。
模板：

int ran[8]={0,2,3,5,7,11,13,17};

int quick(int p,int n,int mod){
	int r=1;
	p=p%mod;
	while(n){
		if(n&1)	r=r*p%mod;
		p=p*p%mod;
		n>>=1;
	}
	return r;
}

bool wit(int c,int s){
	int t=s-1,j=0;
	while(t%2==0){
		t/=2;
		j++;
	}
	int k=quick(c,t,s);
	if(k==1||k==s-1)	return true;
	while(j--){
		k=k*k%s;
		if(k==s-1)	return true;
	}
	return false;
}

bool m_r(int s){
	if(s==1)	return false;
	for(int i=1;i<=7;i++)
		if(s==ran[i])	return true;
	for(int i=1;i<=7;i++)
		if(s%ran[i]==0)	return false;
	for(int i=1;i<=7;i++)
		if(!wit(ran[i],s))	return false;
	return true;
}