线性筛素数

int prime_len;
int primes[n];
bool flag[n];
void init() {
    memset(flag , 0 , sizeof(flag));   
    prime_len = 0;
    for (int i=2 ; i<n ; i++) {
        if (!flag[i]) primes[prime_len++] = i;
        for (int j=0 ; i*primes[j] < n; j++) {
            flag[i*primes[j]] = true;     
            if (i%primes[j]==0) break;
        }
    }
}

证明为什么这样是对的：

     ① 对于任何一个素数p，我们都不可能使用两个数乘积进行标记，所以对于每一个素数，我们都是能找到的。

     ② 对于任何一个合数m = p₁^a1p₂^a2…p_m^am,这里p₁< p₂ < … <p_m，我们都能使用p₁^a1-1p₂^a2…p_m^am  *  p₁进行筛选，所以对于每一个合数，我们都能筛选到做上标记。

每个数只去筛以它为最大因数的数，故每个合数只会被一个数筛（不会被重复筛），就是它的最大因数。