质数筛 （学习笔记）

埃氏筛法

埃氏筛法的思想就是：利用当前已经找到的素数，从后面的数中筛去当前素数的倍数，由预备知识一可知，当前素数已经是筛去数的质因子，如此下去能筛除所有之后的合数，是一种比较快的筛法 时间复杂度：O(nloglogn）

（已知求得a的平方根b，a的因式分解所得因数c，d必有c>b&&d<b,或 c<b&&d>b)

实现代码如下：

#include<iostream>
#define maxn 1000010
#define ll long long
using namespace std;
int prime[maxn];
bool is_prime[maxn];
int sieve(int n)
{   
    memset(is_prime,true,sizeof(is_prime));//注意初始化 
    int p=0;
    is_prime[0]=is_prime[1]=false;
    for(int i=2;i<=n;i++)
    {
        if(is_prime[i])
        {   
            prime[p++]=i;
            for(int j=2*i;j<=n;j+=i)//美中不足是会筛除多次比如8和16，同时被2和4筛去
                is_prime[j]=false;              
        }
 
    }
return p;   
}
int main()
{
    int n;
    cin>>n;
    cout<<sieve(n)<<endl;
 
} 

欧几里得筛法

这是一种很好的线性筛法，和埃氏筛法的区别是对于每一个要筛除的数，欧拉筛法只筛除一次，而埃氏筛法会重复筛除，比如8和16同时被2和4筛去，推荐使用欧拉筛法

部分实现

memset(st,0,sizeof(st));
st[1]=0;
for(i=2;i<=n;i++)
{
    if(st[i])
    primes[cnt++]=i;
    for(j=0;primes[j]*i<=n&&j<=cnt;j++)
    {
        st[primes[j]*i]=0;
        if(i%primes[j]==0)
        break;
    }
}