浅谈线性筛

众所周知，Eratosthenes筛的复杂度是(O(nloglog n))的，并非严格线性。于是有时候就要用到线性的欧拉筛。

---

欧拉筛的原理是这样的：

对于每一个数，我们只用它最小的质因子去筛掉它，因此每一个数只会被筛掉一次，从而保证该筛法为线性的。

在实现时，我们需要用数组(flag_{1,n})去表示每一个数是否为质数，并把所有质数从小到大放入数组(pri_{1,cnt})中。若在区间([1,i])共有(cnt)个质数，对于(i)，我们用(i)和(pri[j])去筛掉(x = i imes pri[j]),其中(j in [1,min(cnt,k)], i mod pri[k] = 0)。在(i mod pri[k] = 0) 时停止继续往下筛，是因为当(j > k)时，(pri[j]>pri[k])，而(pri[k])是(i)的一个因子，也就意味着(pri[j])不再是(x = i imes pri[j])的最小的质因子了，故我们不再往下继续筛，以避免重复计算。

int pri[maxn];    // 存筛出来的素数 

bool flag[maxn];  // flag[i] = false 表示i是素数，flag[i] = true 表示i不是素数 

void work() {
	flag[1] = true;
	for(int i = 2; i <= n; i++) {
		if(!flag[i]) pri[++cnt] = i;
		for(int j = 1; j <= cnt && i * pri[j] <= n; j++) {
			flag[i * pri[j]] = true;
			if(i % pri[j] == 0) break;   // pri[j + 1] 一定不是 i * pri[j + 1] 最小的因子了，停止继续往下筛 
		}
	}
}