欧拉素数筛详解 （洛谷 P3383 【模板】线性筛素数 题解）

什么是欧拉素数筛？

欧拉素数筛是一种十分高效的筛选出素数的的方法，其复杂度为O（n）
有人会问，求一个数是否是素数不也是O（n）复杂度吗？
但是，欧拉筛的作用，是在O（n）复杂度筛出2到n的所有素数
这就是欧拉筛的神奇之处，强行把筛的复杂度降低到O（n）

核心思想

实际上，欧拉筛是一种名叫埃氏筛法的优化版

首先我们说一下埃氏筛法的思想

其实，这是一种预处理的思想
我们知道，一个数的是否是质数，在于它是否有除1和它本身以外的质因数
所以，我们可以通过这样一种方法快速去除不是素数的数：
当我们找到一个素数，便将这个素数乘以倍数 ，产生的数一定是一个非素数，我们便将它标记，就这样一遍一遍筛下去，最后剩下没被标记的数就是素数了
复杂度是O（nloglogn）

埃氏筛法的缺点

我们注意到，埃氏筛法有一个很大的缺点：它将一个非素数重复筛取了
例如：

18
我们在找到素数2的时候，将它乘9的时候已经筛过一次18，18被标记为了非素数
在我们找到素数3的时候，我们同样将3乘6的时候也筛了一次18
这样，一个数筛了多次，导致没有到达线性复杂度

这就是埃氏筛法的缺点

欧拉筛做出的优化

欧拉筛的优化就在于避免了重复筛取，使筛数的复杂度为线性复杂度
那它是如何做到的呢？
我们注意到，当一个素数乘以一个数的时候，如果这个数是这个素数的倍数，那么这个数乘任何数得到的数，也一定是这个素数的倍数
听起来有点绕，我们举个例子：

比如12由26得到，我们知道2一定会是12的质因数
但是，2同样是6的质因数
所以，如果我们继续用6乘一个素数筛取其他数
例如36=18，这个18 因为是由素数3乘6倍得到，但是18的最小因数是2
他应该由2*9筛出来，这样便筛了两次，复杂度增加了

所以，我们到达所乘倍数是素数的倍数的时候，就终止，这样，每个数被筛一次，复杂度便是O（n）

欧拉筛的模板

bool is_prime[MAX_N];		//用来筛取的数组，用作标记
int prime[MAX_N];		//将素数存入其中
int cnt = 0;		//记录素数的个数

void sieve() {
  for (int i = 2; i < MAX_N; i++) {
    if (!is_prime[i]) prime[cnt++] = i;		//没有被标记便存入prime数组
    for (int j = 0; j < cnt; j++) {
      if (i * prime[j] >= MAX_N) break;		//当筛到最大数时，便停止
      is_prime[i * prime[j]] = true;	//标记为非素数
      //停止条件，也是欧拉筛的核心
      if (i % prime[j] == 0) break;
    }
  }
}

模板题

题目描述

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

输入格式

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

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

输出格式

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

输入

100 5
2
3
4
91
97

输出

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。

代码如下：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
using namespace std;

const int MAX_N = 10000005;
bool is_prime[MAX_N];
int prime[MAX_N];
int cnt = 0;
int N, M;

void sieve() {
  for (int i = 2; i < MAX_N; i++) {
    if (!is_prime[i]) prime[cnt++] = i;
    for (int j = 0; j < cnt; j++) {
      if (i * prime[j] >= MAX_N) break;
      is_prime[i * prime[j]] = true;
      if (i % prime[j] == 0) break;
    }
  }
}

void input() {
  cin >> N >> M;
}

void solve() {
  int number;
  for (int i = 0; i < M; i++) {
    cin >> number;
    if (number == 1) cout << "No" << '\n';		//注意1时候的特判
    else if (!is_prime[number])
      cout << "Yes" << '\n';
    else
      cout << "No" << '\n';
  }
}

int main() {
  input();
  sieve();
  solve();
  return 0;
}