四种素数筛法：朴素素数筛，埃氏筛，欧拉筛和区间筛详解

预备知识：

1.算术基本定理：任何一个大于1的自然数 N,如果N不为质数，那么N可以唯一分解成有限个质数的乘积
2. 若一个数可以进行因数分解，则得到的两个数一定是有一个>=sqrt(x),另一个<=sqrt(x).

朴素算法

复杂度：O(n√n)
介绍：这个算法是最简单的素数判断算法+遍历素组，耗时长

#include<bits/stdc++.h>
#define maxn 1000010
#define ll long long
using namespace std;
int prime[maxn];
bool judge(ll x)
{   
    if(x==2) return true;
    if(x<2||x%2==0)
        return false;   
    for(int i=3;i<=sqrt(x+1);i+=2)//这个地方可以优化一下，一个是范围，一个是步长 （素数一定是奇数） 
        {
            if(x%i==0)
                return false;
        }
    return true;
}
int sieve(int n)
{
    int p=0;
    for(int i=0;i<=n;i++)
        if(judge(i))
            {
            prime[p++]=i;}
    return p;   
}
int main()
{
    int n;
    cin>>n;
    cout<<sieve(n)<<endl;

}

2.埃氏筛法

复杂度：O(nloglogn)
介绍：埃拉托斯特尼筛法，利用当前已经找到的素数，从后面的数中筛去当前素数的倍数，由预备知识一可知，当前素数已经是筛去数的质因子，如此下去能筛除所有之后的合数，是一种比较快的筛法

#include<bits/stdc++.h>
#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;

}

3.欧几里得筛法

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

#include<bits/stdc++.h>
#define maxn 1000010
#define ll long long
using namespace std;
int prime[maxn];
int factor[maxn];//用来记录最小素因子 
int sieve(int n)
{   
    int p=0;
    for(int i=2;i<=n;i++)
    {
        if(!factor[i])//如果没有找到素因子，当然就是素数啦 
        {
            prime[p++]=i;
            factor[i]=i;            
        }
        for(int j=0;j<p&&prime[j]*i<=n;j++)
        {
            factor[prime[j]*i]=prime[j];
            if(!(i%prime[j]))//筛的数已经被前面的数筛过了 
                break;      
        }   
    }
    return p;   
}
int main()
{
    int n;
    cin>>n;
    cout<<sieve(n)<<endl;

}

区间筛法

介绍：对于于超大区间，如a<’b<=10^12,b-a<=10^6,这种情况下，数组是开不到这么大的，所以要用偏移来筛选，具体是先打一个2-√b的小表,然后用此表来筛除a到b的素数

#include<bits/stdc++.h>
#define ll long long
#define maxn 1000005 
using namespace std;
bool is_prime[maxn];
bool is_prime_small[maxn];
ll ans=0;
void segment_sieve(ll a,ll b)
{   
    for(ll i=0;i*i<b;i++)   is_prime_small[i]=true;
    for(ll i=0;i<b-a;i++)   is_prime[i]=true;

    for(ll i=2;i*i<b;i++)
    {
        if(is_prime_small[i])
            {
                for(ll j=2*i;j*j<b;j+=i) is_prime_small[j]=false;
                for(ll j=max((long long)2,(a+i-1)/i)*i;j<b;j+=i) is_prime[j-a]=false;                   

            }       

    }
} 
int main()
{

    ll a,b;
    cin>>a>>b;
    segment_sieve(a,b);
    for(ll i=0;i<b-a;i++)   
        if(is_prime[i])
            ans++;  

    cout<<ans<<endl;    

}