Eratosthenes筛法
原文链接:https://blog.csdn.net/too_handsome_to_live/article/details/98793838
埃拉托斯特尼筛法,简称埃氏筛或爱氏筛,是一种由希腊数学家埃拉托斯特尼所提出的一种简单检定素数的算法。
用于求得[1, n]区间内的全部素数。
算法流程:
第一步,将[2, n]区间排成一列。
第二步:标出列表中的第一个数,筛去其所有的倍数。
不断重复第二步,被标记的就是[1, n]区间的所有素数。
例如:求[1,10]区间内的所有素数。
- 将[2, 10]排成一列
- 标出列表中的第一个数——2,筛去所有2的倍数
- 标出列表中的第一个数——3,筛去所有3的倍数
- 标出列表中的第一个数——5,筛去所有5的倍数
- 标出列表中的第一个数——7,筛去所有7的倍数
- 被标记的就是[1, 10]区间内的所有素数:2,3,5,7
从这个过程中可以看出,每次筛完后剩下的区间内第一个数一定是素数,因为合数一定是前面某一个素数的倍数,一定在前几次筛选的过程中被筛走了。
#include<iostream>
bool numlist[100000001];
int prime[20000001], cnt;
void work(int n)
{
for (int i = 2; i <= n; i++)
{
if (numlist[i] == false)
{
prime[++cnt] = i;
for (int j = i; i * j <= n; j++)
numlist[i * j] = true;
}
}
return;
}
埃氏筛的思想是:要想得到 n以内的质数,就要把不大于 sqrt(n) 的质数的倍数全部剔除,剩下的就是质数。从 2开始,把 2 的倍数(不包括本身)标记为合数,然后向后枚举,查到一个未标记为合数的,就把它的倍数(不包括本身)标记为合数。以此类推,查到 n为止。
有一个小优化,把 t的倍数标记为合数时,不是从 2t开始,而是从 t^2 开始(因为小于 t^ 的 t的倍数在枚举到 t前就被标记过了)。
当然,埃氏筛效率还是低了些。例如一个数 24,它会被 2, 3, 4 三个数标记,这就重复了两次,更大的数同理。时间复杂度 O(n log log n),对于 10^8的数据,会超时。
欧拉筛
原文链接:https://www.luogu.com.cn/problem/solution/P3383
现在要避免重复筛,要用到另一种筛法:欧拉筛。
避免重复筛,应找到筛合数的一种原则:这个合数只会被它的最大非自身因数(对应最小质因数)筛。这样能保证每个合数只会被筛一次。
(运行过程)从 2 开始,2 加入 prime 数组,再从小到大枚举质数(现在只有 2),筛掉质数与 2 的乘积(4 被筛掉)。
到了 3,3 加入 prime 数组,从小到大枚举质数(此时有 2,3),筛掉质数与 3 的乘积(6,9 被筛掉)。 到了 4,4 没加入 prime 数组,枚举质数(有2,3),筛掉 8 后,因为 4 mod 2 = 0,触发退出条件。(若不触发,就会筛掉 12,而 12 = 4x3 = 6x2,又会被 6*2筛一次)
以此类推,可做出一张表:
i 的值 | 质数表 | 筛去的数 |
---|---|---|
2 | 2 | 4 |
3 | 2, 3 | 6, 9 |
4 | 2, 3 | 8 |
5 | 2, 3, 5 | 10, 15, 25 |
`6 | 2, 3, 5 | 12 |
7 | 2, 3, 5, 7 | 14, 21, 28, 35 |
… | … | … |
每个质数只被筛一次,复杂度变为 O(n)O(n) ,可以AC。 |
#include<iostream>
using namespace std;
//求第k小的素数是多少
//欧拉筛
bool numlist[100000001]; //通过这个数组寻找素数
int prime[20000001]; //存放素数
int cnt; //计数,标记是第几个(第i小的)素数
void getprime(int n) //通过此函数从小到大依次找到小于n的素数,存到prime[]数组中
{
for (int i = 2; i <= n; i++)
{
if (numlist[i] == 0) // numlist[i]等于0代表i是素数
{
cnt++;
prime[cnt] = i; //第cnt个素数是i
}
//从第一个素数开始,直到第cnt个素数时停止,并且筛掉的素数要在查找范围n内
for (int j = 1; j <= cnt && i * prime[j] <= n; j++)
{
numlist[i * prime[j]] = 1; //把所有合数标记为1(即i*prime[j]不是素数)
if (i % prime[j] == 0) //避免重复筛
break;
}
}
}