筛质数

给定一个n，求1-n内有多少个质数

从2开始枚举，如果某一个数是质数，那么对所有的2~n上所有的数，只要是该质数是当前数的约数，就把当前数剔除出去。

代码实现

#include<iostream>

using namespace std;

const int N = 1000010;

int n;
int prime[N], cnt;
bool st[N];

void select_prime(int n)
{
    for(int i = 2; i <= n ;i++){
        if(!st[i]) prime[cnt++] = i;
        
        for(int j = i + i; j <= n ; j += i)
            st[j] = true;
    }
}
int main(){
    cin >> n;
    
    select_prime(n);
    
    printf("%d\n", cnt);
    return 0;
}

我们来分析时间复杂度， 一共有n个数，从2开始枚举，循环总数是
n / 2 + n/3 + n/4 … + n/ n, 提出公因式n，可以得到一个调和级数1 + 1/2 + 1/3 + … + 1/n, 调和级数的最值是logn，所以算法的复杂度是
nlogn。
进一步优化：当且仅当i是质数的时候，才开启筛选。
由质数相关定理得到，1~n当中有 n / (ln n)个质数。粗略估计这样整体的时间复杂度就是n ln n / ln n = n。最终的复杂度是n loglogn

这个筛选的方法叫埃氏筛法。

#include<iostream>

using namespace std;

const int N = 1000010;

int n;
int prime[N], cnt;
bool st[N];

void select_prime(int n)
{
    for(int i = 2; i <= n ;i++){
        if(!st[i]) {
            prime[cnt++] = i;
             
            for(int j = i + i; j <= n ; j += i)
            st[j] = true;
        }
    }
}
int main(){
    cin >> n;
    
    select_prime(n);
    
    printf("%d\n", cnt);
    return 0;
}

还有一个线性筛法。

#include<iostream>

using namespace std;

const int N = 1000010;

int n;
int prime[N], cnt;
bool st[N];

void select_prime(int n)
{
    for(int i = 2; i <= n ;i++){
        if(!st[i]) {
            prime[cnt++] = i;
             
        }
         
        for(int j = 0; prime[j] <= n / i ; j ++)    //从小到大枚举所有的质数
        {
            st[prime[j] * i] = true;
            if(i % prime[j] == 0) break;        //prime[j]一定是i的最小质因子， 因此只会筛选一次，线性的
        }
    }
}
int main(){
    cin >> n;
    
    select_prime(n);
    
    printf("%d\n", cnt);
    return 0;
}