线性筛素数 算法笔记整理+模版

我们以洛谷上的这道题为例题来讲：传送门

题目描述

如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）

输入输出格式

输入格式：

第一行包含两个正整数N、M，分别表示查询的范围和查询的个数。

接下来M行每行包含一个不小于1且不大于N的整数，即询问该数是否为质数。

输出格式：

输出包含M行，每行为Yes或No，即依次为每一个询问的结果。

输入输出样例

输入样例#1： 

100 5
2
3
4
91
97

输出样例#1： 

Yes
Yes
No
No
Yes

说明

时空限制：500ms 128M

数据规模：

对于30%的数据：N<=10000，M<=10000

对于100%的数据：N<=10000000，M<=100000

样例说明：

N=100，说明接下来的询问数均不大于100且不小于1。

所以2、3、97为质数，4、91非质数。

故依次输出Yes、Yes、No、No、Yes。

相信大家都会试除法筛素数（时间复杂度O(n*sqrt(n))）普通筛素数（时间复杂度O(n*logn)）但是如果想在O(n)范围内筛出所有的素数怎么办呢？这里我们要引入一个叫做欧拉筛法（怎么哪都有欧拉......）的神奇算法。大家有没有考虑过为什么普通筛素数是O(n*logn）时间复杂度的呢？因为每个合数被自己的所有素因子都筛了一遍，等于这个合数被访问了很多次，因此时间复杂度就很高。如果我们能保证每个合数只被筛了一次，即可做到O(n)解决问题。这里我们需要明白一个原理：一个合数一定至少有一个能整除它的素因子。我们只要找到最小的那个素因子，保证每个合数只被它最小的那个素因子筛了一次，也就是只被访问一次，就能做到线性筛素数。

AC代码：

#include<cstdio>
#include<iostream>
using namespace std;
const int Maxn=10000005;
bool isp[Maxn];
int pri[Maxn/10],cnt,n,m,tx;
void makepri(){//筛素数 始终贯彻的要求是只被最小的素因数筛 
	for(int i=2;i<=Maxn;i++){
		if(!isp[i])pri[++cnt]=i;//如果它之前未被筛过，那么它一定是素数，放入素数表中 
		for(int j=1;j<=cnt&&pri[j]*i<=Maxn;j++){
			isp[pri[j]*i]=1;
			if(i%pri[j]==0)break;//能找到更小的素数来筛pri[j+1]*i 
		}
	}
}
inline int read(){//快速读入 
	char c=getchar();int x=0;
	while(c<'0'||c>'9')c=getchar();
	while(c>='0'&&c<='9'){x=x*10;x=x+c-'0';c=getchar();}
	return x;
}
int main(){
	n=read();m=read();
	makepri();isp[1]=1;
	for(int i=1;i<=m;i++){
		tx=read();
		if(!isp[tx])printf("Yes\n");
		else printf("No\n"); 
	}
	return 0;
}