我们以洛谷上的这道题为例题来讲:传送门
题目描述
如题,给定一个范围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;
}