简介
埃拉托斯特尼筛法(sieve of Eratosthenes ),简称埃氏筛,是一种简单且年代久远的算法,用来找出一定范围内所有的素数。
所使用的原理是从2开始,将每个素数的各个倍数,标记成合数。一个素数的各个倍数,是一个差为此素数本身的等差数列。此为这个筛法和试除法不同的关键之处,后者是以素数来测试每个待测数能否被整除。
算法
给出要筛数值的范围n,。先用2去筛,即把2留下,把2的倍数剔除掉;再用下一个素数,也就是3筛,把3留下,把3的倍数剔除掉;接下去用下一个素数5筛,把5留下,把5的倍数剔除掉;不断重复下去……。
步骤
详细列出算法如下:
列出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 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 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
如果现在这个序列中最大数小于最后一个标出的素数的平方,那么剩下的序列中所有的数都是质数,否则回到第二步。
本例中,因为25大于2的平方,我们返回第二步:
剩下的序列中第一个质数是3,将主序列中3的倍数划出(红色),主序列变成:
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,3
25仍然大于3的平方,所以我们还要返回第二步:
现在序列中第一个质数是5,同样将序列中5的倍数划出,主序列成了:
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 3 5 。
因为25等于5的平方,跳出循环.
结论:去掉红色的数字,2到25之间的质数是:2 3 5 7 11 13 17 19 23。
伪代码
Input: an integer n > 1
Let A be an array of Boolean values, indexed by integers 2 to n,
initially all set to true.
for i = 2, 3, 4, ..., not exceeding √n:
if A[i] is true:
for j = i2, i2+i, i2+2i, i2+3i, ..., not exceeding n :
A[j] := false
Output: all i such that A[i] is true.
以上算法可以得到小于等于n的所有素数,它的复杂度是O(n log log n)。
代码
**python 3.6**
def eratosthenes(n):
P = [i for i in range(2, n+1)]
p = 0
while True:
for i in P[p + 1:]:
if i % P[p] == 0:
P.remove(i)
if P[p]**2 >= P[-1]:
break
p += 1
return P
if __name__ == "__main__":
print (eratosthenes(120))
**C++**
#include <vector>
#include <algorithm>//std::remove_if
#include <numeric>//std::iota
std::vector<int> eratosthenes(int max){
std::vector<int> vi(max+1);//0, 1, 2, ..., max
std::iota(vi.begin(), vi.end(), 0);
if(max>=2){
int prime=2;
while(prime*prime<=max){//2 to sqrt(max)
for(size_t index=prime*2; index<vi.size(); index+=prime){
vi[index]=0;//Rule out this number.
}
for(prime++; prime*prime<=max; prime++){
if(vi[prime]>0){
break;//Jump to next non-zero.
}
}
}
}
vi.erase(std::remove_if(vi.begin(), vi.end(), [](int i)->bool{
return i<=1;
}), vi.end());//Remove all zeros and the one.
return vi;
}
以上内容来自维基百科