对于线性筛法筛质数的理解

开了两个数组，一个是st[N],初始值全为false，在筛选过程中如果i是合数则对应的st[i]变为true。primes[N]数组中存放1~n中的质数。

从线性筛法的代码中可以看出，在第一重循环中，会先通过st[i]的逻辑值判断该数是不是质数，然后i会乘于目前的primes数组中的各个数，直到i会被primes数组中的某个数整除。

可以从反证法来理解这个算法。

首先假设1~n中的某个质数a,st[a]从false改为true，这显然是不可能的，因为a是质数，a前的任何数乘于primes数组中的任何数都不可能等于a,所以st[a]不可能为true.

然后就是假设1~n的某个合数b,假设st[b]并没有从false变为true。那么也就是代码中的

st[primes[j] * i] = true 中的 primes[j] * i 从来没有等于b过。可我们可以对b进行质因数分解进行分析。假设b分解后的质因数从小到大依次为 b1 b2 b3 b4 b5 ...bn,由上面第一个反证可以知道b1质数会被放在prime数组中，那么当第一重循环i = b2*b3*b4*b5*...*bn时，当此时的i进入第二个循环时，i绝对不会被primes数组中b1之前的元素break掉（因为b2,b3,b4...bn都是质数，它们的乘积自然不会被别的质数整除），那么就一定会执行到st[b1*i] = true这个语句，也就是st[b]一定会变为true.反证成功。

//线性筛法
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1000010;
bool st[N];
int cnt, primes[N];

int get_primes(int n)
{
    for(int i = 2; i <= n; i++)
    {
        if(!st[i]) primes[cnt++] = i;
        for(int j = 0; primes[j] <= n / i; j++)
        {
            st[primes[j] * i] = true;
            if(i % primes[j] == 0) break;
        }
    }
    return cnt;
}

int main()
{
    int n;
    cin >> n;
    cout << get_primes(n);
    return 0;
}