题解 P3383 【【模板】线性筛素数】

本质上不是很难想，这里给大家列两个筛法

一、埃式筛法： 1.能在n * log（logn）的时间内筛出所有质数

2.其基本思路为利用给定范围n内的质数将所有合数筛出来

3.目前NOIP范围内掌握这个筛法即可

上代码：

int m = sqrt(n) + 1;
memset(vis,0,sizeof(vis));//先清空数组
for(int i = 2;i <= m;i++)
if(!vis[i])//如果没被标记过
for(int j = i * i;j <= n;j += i)
vis[j] = true;//有合数就标记，输出时进行判断就好
二、欧拉筛： 1.线性筛出质数

2.原理在于每个非素数只会被最小的素数因子筛一次

3.欧拉筛还可以用来筛很多其他积性函数（我并不会…）

上代码：

long long n;//n为要筛的范围
int prime[n + 5],vis[n + 5],cnt;
void find()
{
int i,j;
for(int i = 2;i <= n;i++)
{
if(!vis[i])
prime[++cnt] = i;
for(int j = 1;prime[j] * i <= n && j <= cnt;j++)
{
vis[i * prime[j]] = 1;
if(i % prime[j] == 0)
break;
}
}
}
欧拉筛有能力的话尽力了解吧，埃式筛是必须要掌握的！！！

上本题代码（我用的埃式筛）：

#include
#include
#include
#include
#include
#define rep(i,a,n) for(int i=a;i<=n;i++)//定义循环，个人习惯
typedef long long ll;
using namespace std;

ll read(){
ll ans=0;
char last=’ ‘,ch=getchar();
while(ch<‘0’ || ch>‘9’)last=ch,ch=getchar();
while(ch>=‘0’ && ch<=‘9’)ans=ans*10+ch-‘0’,ch=getchar();
if(last==’-’)ans=-ans;
return ans;
}//快读，可忽略，个人习惯，同上

ll num[100005];
bool vis[10000005];
int main()
{
long long n,m;
n = read(),m = read();
long long k = sqrt(n) + 1;
memset(vis,0,sizeof(vis));
rep(i,2,k)
{
if(!vis[i])
for(int j = i * i;j <= n;j += i)
vis[j] = true;
}
rep(i,1,m)
{
num[i] = read();
}
rep(i,1,m)
{
if(vis[num[i]] == true || num[i] == 1) printf(“No\n”);//注意要特判1，1为非素数
else printf(“Yes\n”);
}
return 0;