埃氏筛算法解读（c++)

1.前言

本文将介绍埃氏筛的定义和应用。

我会先用纯数学知识去推导此算法，再回归到编程中。

2.什么是埃氏筛？

埃氏筛，全名为埃拉托斯特尼筛法（Sieve of Eratosthenes），是一种用于找出一定范围内所有质数的古老算法。它是由古希腊数学家埃拉托斯特尼在公元前3世纪发明的。

3.埃氏筛的推导过程

1.列出数字：首先，列出从2开始的所有小于等于给定数的自然数。

2.选择质数：选择列表中的第一个数（它是质数），然后删除它的所有倍数。

3.重复筛选：重复步骤2，选择列表中的下一个未被删除的数（它是下一个质数），然后删除它的所有倍数。

4.结束条件：重复这个过程，直到达到根号 n（ n 为需要筛选的最大数）。

举例（查找2~25之间的质数）：

1.列出2以后的所有序列：

• 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

2.标出序列中的第一个质数，也就是2，序列变成：

• 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

3.将剩下序列中，划掉2的倍数，序列变成：

• 2 3 5 7 9 11 13 15 17 19 21 23 25

4.如果这个序列中最大数小于最后一个标出的质数的平方，那么剩下的序列中所有的数都是质数，否则回到第二步。

5.本例中，因为25大于2的平方，我们返回第二步：

6.剩下的序列中第一个质数是3，将主序列中3的倍数划掉，主序列变成：

• 2 3 5 7 11 13 17 19 23 25

7.我们得到的质数有：2，3

8.25仍然大于3的平方，所以我们还要返回第二步：

9.序列中第一个质数是5，同样将序列中5的倍数划掉，主序列成了：

• 2 3 5 7 11 13 17 19 23

10.我们得到的质数有：2，3，5

11.因为23小于5的平方，跳出循环

结论：2到25之间的素数是：2 3 5 7 11 13 17 19 23。 

4.埃氏筛的优缺点

优点：

1. 相对高效：与其他简单的质数筛选算法相比，埃氏筛在处理大量数据时更加高效。
2. 易于理解和实现：埃氏筛的原理简单直观，易于编程实现。
3. 内存使用：埃氏筛在筛选过程中不需要保存所有非质数，只需要知道它们是否被筛选掉，因此内存使用相对较少。

缺点：

1. 空间复杂度：尽管埃氏筛的内存使用比一些算法少，但当数据量非常大时，它仍需要较大的内存空间来存储筛选过程中的状态。复杂度为O(n)。
2. 时间复杂度：埃氏筛的时间复杂度较高，特别是在处理大量数据时。复杂度为：O(n log log n)。
3. 并行计算限制：埃氏筛的串行特性限制了其在并行计算环境中的效率。虽然可以对其进行并行化改进，但这通常需要复杂的算法调整。

5.代码框架

这段代码实现找到从1到n的所有质数，并将这些质数换行输出的代码：

#include <bits/stdc++.h>
using namespace std;
bool isnp[100005];
int n,m,p[10005];
void init(int n)//埃氏筛
{
	m=0;
	for (int i=2;i<=n;i++)
	{
		if (!isnp[i])//如果i没被标记，说明它是质数
		{
			m++;
			p[m]=i;
			for (int j=2*i;j<=n;j+=i)//把除自身以外的倍数都标记成合数
			{
				isnp[j]=true;
			}
		}
	}
}
int main()
{
	int n;
	cin>>n;
	init(n);
	for (int i=1;i<=m;i++)//输出质数
	{
		cout<<p[i]<<" ";
	}
}

6.后记：

我们从零讲解了埃氏筛的用法，并用数学方法推导出了埃氏筛的过程。

如果喜欢的话请到我主页阅读，后续我会更新更多关于c++算法知识的博客，我以后会做得更好。

我是蒟蒻，请神犇们多加指点。