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.埃氏筛的优缺点
优点:
- 相对高效:与其他简单的质数筛选算法相比,埃氏筛在处理大量数据时更加高效。
- 易于理解和实现:埃氏筛的原理简单直观,易于编程实现。
- 内存使用:埃氏筛在筛选过程中不需要保存所有非质数,只需要知道它们是否被筛选掉,因此内存使用相对较少。
缺点:
- 空间复杂度:尽管埃氏筛的内存使用比一些算法少,但当数据量非常大时,它仍需要较大的内存空间来存储筛选过程中的状态。复杂度为O(n)。
- 时间复杂度:埃氏筛的时间复杂度较高,特别是在处理大量数据时。复杂度为:O(n log log n)。
- 并行计算限制:埃氏筛的串行特性限制了其在并行计算环境中的效率。虽然可以对其进行并行化改进,但这通常需要复杂的算法调整。
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++算法知识的博客,我以后会做得更好。
我是蒟蒻,请神犇们多加指点。