【欧拉筛】筛质数洛谷素数个数含代码讲解

残夏丰秋源

于 2024-07-14 17:08:53 发布

阅读量759

点赞数 24

分类专栏：小算法文章标签：算法 c++

本文链接：https://blog.csdn.net/wayua/article/details/140414732

版权

小算法专栏收录该内容

2 篇文章 0 订阅

订阅专栏

欧拉筛洛谷素数个数

-欧拉筛（Euler’s Sieve），又称为埃氏筛法的优化版，是另一种高效的查找一定范围内质数的方法。它是在线性筛的基础上对部分步骤进行了优化，特别适合寻找较大的质数集合。欧拉筛通常用于大型数据集上，如求解最大不超过某个数n的所有质数的和等问题。

洛谷这道题正可以利用欧拉筛法 👉洛谷P3912素数个数👇

素数个数

题目描述

求 $1,2,\cdots,N$ 中素数的个数。

输入格式

一行一个整数 $N$ 。

输出格式

一行一个整数，表示素数的个数。

样例 #1

样例输入 #1

样例输出 #1

提示

对于 $40\%$ 的数据， $\le N \le 10^6$ 。

对于 $80\%$ 的数据， $\le N \le 10^7$ 。

对于 $100\%$ 的数据， $\le N \le 10^8$ 。

筛法总共大致可以分成三类：

一、无算法筛：
也就是最简单也最慢的筛法，原理是将n个数从头到尾每一个数都判断一遍，（判断要重复至少n/2次），所以基本上做题都会超时，不合适。

二、埃氏筛：
埃拉托斯特尼筛（Sieve of Eratosthenes），又称简单素数筛或欧几里得筛，是一种古老而有效的找出一定范围内所有质数的经典算法。其基本原理如下：

首先创建一个从2开始连续的整数序列，并假设它们都是质数。
然后从第一个质数2开始，将序列中所有2的倍数标记为合数（非质数）。
接下来选择下一个未被标记的数（这里是3），它是新的质数，再将序列中所有3的倍数标记为合数。
重复这个过程，每次选择当前未被标记的最大质数，直到这个质数大于当前范围内的最大数为止。
通过这个过程，剩下的未被标记的数就是该范围内的质数。

埃氏筛非常适合寻找小范围内的质数，但对于大范围，由于需要存储大量信息，可能会占用较多内存。然而，对于优化版本如欧拉筛，可以在节省空间的同时保持效率较高。

三、欧拉筛（线性筛）：
欧拉筛就是埃氏筛的改进版，具体原理如下：

首先，列出从2开始的自然数序列。
将2作为第一个素数，然后剔除掉2的所有倍数。
接着，找到下一个未被剔除的数，将其作为第二个素数，并剔除掉它的所有倍数。
重复上述过程，直到所有小于指定数的素数都被找出。

欧拉筛（Sieve of Euler）和埃氏筛（Sieve of Eratosthenes）都是用来找出一定范围内所有素数的筛选算法，但它们在操作原理上存在一些关键的区别：

筛选过程：
- 埃氏筛：从最小的素数开始，将其所有的倍数标记为合数。对于每个素数，都会重复这个过程，直到达到预定的范围上限。
- 欧拉筛：在筛选过程中，每个合数会被其最小的素因子筛选一次，并且仅筛选一次。这意味着每个合数在筛选过程中只与一个素数相关联。
重复筛选：
- 埃氏筛：可能会多次筛选同一个合数。例如，数字36会被2和3都筛选一次，因为它既是2的倍数也是3的倍数。
- 欧拉筛：每个合数只被其最小的素因子筛选一次。例如，36只会被2筛选，因为2是最小的素因子。
标记方式：
- 埃氏筛：通常使用一个布尔数组来标记合数，素数保持未标记状态。
- 欧拉筛：在筛选过程中，会记录每个合数的最小素因子，并在筛选时使用这个信息。
时间复杂度：
- 埃氏筛：时间复杂度为 (O(n \log \log n))。
- 欧拉筛：时间复杂度为 (O(n))，这是因为每个数只被其最小的素因子筛选一次。
实现复杂性：
- 埃氏筛：实现相对简单，容易理解。
- 欧拉筛：实现较为复杂，需要额外的逻辑来确保每个合数只被筛选一次。

总结来说，埃氏筛虽是一种直观的筛选方法，但可能会对一些合数进行多次筛选，而欧拉筛则通过确保每个合数只被其最小的素因子筛选一次来提高效率。这使得欧拉筛在处理大范围素数时比埃氏筛更加高效。

下面就是欧拉筛的代码实例（内含注释）👇

#include<bits/stdc++.h>
#define IOS std::ios::sync_with_stdio(false),cout.tie(0),cin.tie(0);//加快输入输出
using namespace std;
const int N = 1e8+1;
long long n,suc = 0;
bool vis[N];//定义标记数组
int cnt = 0,arr[N];

void init(int x)//欧拉筛法函数
{
	for(int i = 2;i <= x;i++)//从2开始筛
	{
		if(vis[i] == 0)		arr[++cnt] = i;//判断标记等于0，就把这个质数存进来
		for(int j = 1;j <= cnt;j++)//在存进来的数中，筛掉它的倍数
		{
			if(arr[j] * i >x)	break;//如果越界，直接停止
			vis[arr[j]*i] = 1;//将它的倍数标记为1，即已经筛过
			if(i % arr[j] == 0)	break;//标记完后结束
		}
	}
}

int main()
{
	IOS;//加快输入输出
	cin>>n;//读入n
	init(n);//筛n内的质数
	cout<<cnt;//输出个数
}