埃氏筛法
时间复杂度 O ( n l o g l o g n ) O(nloglogn) O(nloglogn)
#include <iostream>
using namespace std;
const int N = 1e6+5;
int primes[N];
bool st[N];
int k = 0;
void get_primes(int n){
for(int i = 2; i <= n; i++){
if(!st[i]) {
primes[k++] = i;
for(int j = i+i; j <= n; j += i){
st[j] = true;
}
}
}
}
int main(){
int n;
cin>>n;
get_primes(n);
cout<<k<<endl;
return 0;
}
线性筛法
用最小质因子来筛每个合数,而且每个合数只会被筛一次,所以时间是线性的。
#include <iostream>
using namespace std;
const int N = 1e6+5;
int primes[N];
bool st[N];
int k = 0;
void get_primes(int n){
for(int i = 2; i <= n; i++){
if(!st[i]) {
primes[k++] = i;
}
for(int j = 0; primes[j] <= n/i; j++){
st[primes[j]*i] = true;
if(i % primes[j] == 0) break;//primes[j]是i的最小质因子
}
}
}
int main(){
int n;
cin>>n;
get_primes(n);
cout<<k<<endl;
return 0;
}