题目链接 : 点击查看
题目描述 :
给定一个正整数 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;
}
}
}