笔记
详情见: https://www.acwing.com/solution/content/68127/
代码
#include<iostream>
using namespace std;
const int N = 1e6+10;
bool st[N]; //st[]= false表示素数
int primes[N];
int n;
int cnt;
void getPrimeNumber1(int n) { //朴素筛
for(int i = 2; i <= n; i++) {
if(!st[i]) primes[cnt++] = i; //若是素数
for(int j = i + i; j <= n; j += i) st[j] = true; //若是i的倍数,则不可能是素数
}
}
void getPrimeNumber2(int n) { //埃式筛
for(int i = 2; i <= n; i++) {
if(!st[i]) {
primes[cnt++] = i;
for(int j = i + i; j <= n; j +=i) st[j] = true;
}
}
}
void getPrimeNumber3(int n) { //欧拉筛
for(int i = 2; i <= n; i++) {
if(!st[i]) primes[cnt++] = i;
for(int j = 0; primes[j] <= n / i; j++) { //primes[j]一定是primes[j]*i的最小质因子
st[primes[j] * i] = true;
if(i % primes[j] == 0) break;
}
}
}
int main() {
cin>>n;
getPrimeNumber3(n);
cout<<cnt;
return 0;
}