AcWing868. 筛质数 (朴素筛法、埃氏筛法、线性筛法)
1、题目:
给定一个正整数n,请你求出1~n中质数的个数。
输入格式
共一行,包含整数n。
输出格式
共一行,包含一个整数,表示1~n中质数的个数。
数据范围
1≤n≤106
输入样例:
8
输出样例:
4
2、基本思想:
(1)朴素筛法:不管是合数还是质数,筛掉2~n中每个数的倍数 O(nlogn);
(2)埃氏筛法:因为每个合数都是由有限个质数相乘得到的,仅筛掉质数的倍数 O(nloglogn);
(3)线性筛法:在埃式筛法中,有的合数会被不同的质数筛多次,会浪费时间,而每个合数的最小质数是唯一的(否则就不是最小质数),所以只需要用每个合数的最小的质数将它筛掉一次 即可O(n)。
3、核心步骤:
(1)朴素筛法
void GetPrimes(int n)
{
for (int i = 2; i <= n; i ++ )
{
if (!st[i]) p[cnt ++ ] = i;
for (int j = i + i; j <= n; j += i ) st[j] = true;
}
}
(2)埃氏筛法
void GetPrimes(int n)
{
for (int i = 2; i <= n; i ++ )
{
if (st[i]) continue;//是合数就直接进下一层循环
p[cnt ++ ] = i;
for (int j = i + i; j <= n; j += i ) st[j] = true;//筛掉质数的倍数
}
}
(3)线性筛法
void GetPrimes(int n)
{
for (int i = 2; i <= n; i ++ )
{
if (!st[i]) p[cnt ++ ] = i;
for (int j = 0; p[j] <= n / i; j ++ )//p[j] <= n / i移项就是p[j] * i <= n,大于n的自然没必要筛。此处不写成p[j] <= i是为了防止算p[j] * i时越界;不写成p[j] <= sqrt(n)是因为一开始p[j]中还没有存那么多质数,j可能会大于cnt。
{
st[p[j] * i] = true;
if (i % p[j] == 0) break;
//①为什么这里要break?
//因为内层循环是从最小质数开始枚举,此条件成立表明此时p[j]是i的最小质因数(就是说i * p[j]中有至少两个相同的最小质因数p[j],一个是p[j]还有一个在i中),如果这个条件成立还要继续下一次循环,则在下一次循环中p[j + 1]必不可能是i的最小质因数,所以要break。
//②为什么顺序是st[p[j] * i] = true写在前面?
//因为就算此条件不成立,则表明此时p[j]是比i的最小质因数还要小的质数,所以此时的p[j]一定是前面p[j] * i的最小质因数,所以这个条件成立与否都不妨碍前面p[j] * i被它的最小质因数筛掉。
}
}
}
比如,如果是36在i = 18的时候被删,35在i = 7的时候被删,8在i = 4的时候被删。
4、C ++ 代码如下(该代码引用AcWing网站代码)
#include <iostream>
using namespace std;
const int N = 10000010;//比题给多一个0,防越界
int cnt;
int p[N];
bool st[N];
//线性筛法:
void GetPrimes(int n)
{
for (int i = 2; i <= n; i ++ )
{
if (!st[i]) p[cnt ++ ] = i;
for (int j = 0; p[j] <= n / i; j ++ )
{
st[p[j] * i] = true;
if (i % p[j] == 0) break;
}
}
}
//埃氏筛法:
//void GetPrimes(int n)
//{
// for (int i = 2; i <= n; i ++ )
// {
// if (st[i]) continue;//是合数就直接进下一层循环
// p[cnt ++ ] = i;
// for (int j = i + i; j <= n; j += i ) st[j] = true;//筛掉质数的倍数
// }
//}
//朴素筛法:
//void GetPrimes(int n)
//{
// for (int i = 2; i <= n; i ++ )
// {
// if (!st[i]) p[cnt ++ ] = i;
// for (int j = i + i; j <= n; j += i ) st[j] = true;
// }
//}
int main()
{
int n;
cin >> n;
GetPrimes(n);
cout << cnt;
return 0;
}//该代码引用AcWing网站代码
注意事项:
数据量比较小的时候埃氏算法比线性要快一点点,数据量大的时候线性要比埃氏快很多。