素数筛

素数筛

1.朴素筛法O(nlogn)

void get_primes(int n){
    for(int i = 2;i <= n; ++i){
        if(!st[i]){
            prime[cnt++] = i;
        }
        //把i的倍数全部筛掉
        for(int j = i + i;j <= n;j += i)st[j] = true;
    }
}

2.埃氏筛法O($nloglogn$)

void get_primes(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;
        }
    }
}
//当i为质数时，筛掉i的所有倍数
//因为当i为合数时，i的倍数一定已经被i的一个因子筛掉

3.欧拉(线性)筛法O(n)

//核心:n只会被最小质因子筛掉
int cnt; //素数的个数
int v[N]; //存每个数的最小质因子
int prime[N]; //存素数

void get_primes(int n){
    for(int i = 2;i <= n; ++i){
        if(v[i] == 0){
            v[i] = i;
			prime[++cnt] = i;
        }
        for(int j = 1 ; j <= cnt;++j){
            if(prime[j] > v[i] || i * prime[j] > n) break;
            v[prime[j] * i] = prime[j];
        }
    }
}