AcWing 868. 筛质数 (线性筛法)

题目链接 ： 点击查看

题目描述 ：

给定一个正整数 n，请你求出 1∼n 中质数的个数。

输入输出格式 ：

输入

共一行，包含整数 n。

输出

共一行，包含一个整数，表示 1∼n 中质数的个数。

输入输出样例 ：

输入

8

输出

4

题目分析 ：

线性筛法相当于是对埃氏筛法求素数个数的一种优化，线性筛是用最小质因数进行素数筛选，而埃氏筛法则是通过每一个数的倍数进行素数筛选。比如6这个数，用线性筛法只会被筛一次，而用埃氏筛法会分别被2和3各筛一次。

用线性筛求素数个数的具体做法为 ： 在void getPrimes(int n)函数中，最外层依旧是从2遍历到n进行循环，st[i]表示数字i是否被筛除。每一次循环开始时，若是if(!st[i]), 则primes[cnt ++ ] = i，表明i是个素数。不管i是否为素数，我们都要再在for (int j = 0; primes[j] <= n / i ; j ++ )这个循环中，对其他数进行筛选。先让st[primes[j] * i ]= true， 在此可以保证prime[j] 一定是prime[j] * i的最小质因数，因为后面还有一行代码，if(i % primes[j] == 0) break ：i % primes[j] == 0 表明primes[j]是i的最小质因数，如果i % primes[j] != 0 ，prime[j]也一定小于i 的最小质因数，因为primes[j]是从大到小存储素数的。至于这里为什么要进行break，因为当primes[j]成为i的最小质因数之后，在下一次循环得到的primes[j + 1] 就不再是i * primes[j + 1] 的最小质因数了，故循环结束。详见如下代码。

代码 ：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 1e6 + 7;
int primes[N], cnt;
bool st[N];
void getPrimes(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;//筛掉i与质数的乘积(被最小质因子筛去) 
			if (i % primes[j] == 0) break;//  primes[j]为i的最小质因子 
		}
	}
}
int main() {
	int n;
	cin >> n;
	getPrimes(n);
	cout << cnt << endl;
	return 0;
}

---

下面给出线性筛法的相关模板

int primes[N], cnt;     // primes[]存储所有素数
bool st[N];         // st[x]存储x是否被筛掉

void 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;
        }
    }
}