AcWing868. 筛质数 (朴素筛法、埃氏筛法、线性筛法）（C++算法）

AcWing868. 筛质数 (朴素筛法、埃氏筛法、线性筛法）

1、题目：
给定一个正整数n，请你求出1~n中质数的个数。

输入格式
共一行，包含整数n。

输出格式
共一行，包含一个整数，表示1~n中质数的个数。

数据范围
1≤n≤10⁶
输入样例：
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网站代码

注意事项：
数据量比较小的时候埃氏算法比线性要快一点点，数据量大的时候线性要比埃氏快很多。