质数筛(暴力优化+埃氏筛+线性筛)

Hello大家好，今天我们来讲讲质数筛法

看到 质数两个字，给大家的第一个反应是不是这个数，除了1和它本身没有别的数能被ta整除( 真的好孤单)
通过ta的定义，我们就有了第一种写法

1.暴力

我们可以从2~n-1枚举一遍，若其中有数可以被n整除就说明n非质数，没有就是质数

#include<iostream>
using namespace std;
int n;
bool check(int x){
   for(int i=2;i<x;i++){
 	  if(x%i==0)return 0;
   }
   return 1;
}
int main(){
	cin>>n;
	if(check(n))cout<<"YES";
	else cout<<"NO";
	return 0;
}

但是，如果要判断多组数据这个方法就需要优化
众所周知，我们枚举到sqrt(n)即可
因为再往后枚举也是重复
例如：
2 * 3 = 6
3 * 2 = 6
很明显，一个小于sqrt(n)的因子必定有一个大于sqrt(n)的因子

于是就有了如下优化

#include<iostream>
using namespace std;
int n;
bool check(int x){
   for(int i=2;i*i<=x;i++){
 	  if(x%i==0)return 0;
   }
   return 1;
}
int main(){
	cin>>n;
	if(check(n))cout<<"YES";
	else cout<<"NO";
	return 0;
}

现在，如果有t组数据，复杂度O(t*sqrt(n))，还行

2.埃氏筛

例题 洛谷U208323 埃氏筛质数
埃氏筛主要处理n以内的质数，可以一次性处理多个质数，其原理也十分简单
众所周知，一个质数的倍数一定不是质数(除了他本身)，于是就用这个来筛出非质数，剩下的就是质数

#include<bits/stdc++.h>
using namespace std;
bool a[50000005];
int main(){
	memset(a,1,sizeof(a));//先全部都认为是质数
	int m;
	cin>>m;
	for(int i=2;i<m;i++){
		if(a[i]!=0){
			for(int j=i+i;j<m;j+=i){//质数的倍数不是质数(如果这个数是i乘以一个比i更小的数字，那么就一定被那个比i更小的数筛过，循环可以从i*i开始，但是i*i可能会爆int，i和j都需要开long long)
				a[j]=0;
			}
		}
	}
	for(int i=2;i<m;i++){
		if(a[i]!=0)cout<<i<<" ";
	}
	return 0;
}

但是，即使用了ii的优化但仍然有一些重复
例如：
n=12
首先被2筛到，接着会被3筛
虽然不会被6筛了，但重复还是有啊(此时的复杂度为O(nlog(log(n))) )
但我们需要一种真正的线性
那我们有没有一种办法来避免重复呢？

3.于是我们的线性筛闪亮登场！

线性筛的原理举个例子就是
10可以被2筛，就不会把10给5筛，因为这样既不会重复
因为每个合数必有一个最小质因子。每个合数仅被它的最小素因子筛去正好一次。所以为线性时间。代码中体现在：if(i%p[j]==0)break;p数组中的素数是递增的,当 i 能整除 p[j]，那么 i*p[j+1] 这个合数肯定也可以被 p[j] 筛掉，因为i中含有p[j], p[j] 比 p[j+1] 小。接下去的素数同理，所以不用筛下去了。在满足i%p[j]==0这个条件之前以及第一次满足该条件时,p[j]必定是p[j]*i的最小因子

上代码

#include<bits/stdc++.h>
using namespace std;
const int MAX=1e8+5;
bool f[MAX];
int p[MAX];
void a(int n){
	memset(f,true,sizeof(f));
	int cnt=0;
	f[1]=0;
	for(int i=2;i<=n;i++){
		if(f[i])p[++cnt]=i;
		for(int j=1;j<=cnt && i*p[j]<=n;j++){
			f[i*p[j]]=0;
			if(i%p[j]==0)break;
		}
	}
}
int main(){ 
	int n,q,l;
	cin>>n>>q;
	a(n);
	for(int i=0;i<q;i++){
		cin>>l;
		cout<<p[l]<<"\n";
	}
	return 0;
}

END