素数筛（埃拉托斯特尼筛法，欧拉筛，米勒拉宾素数测试，区间筛）

1.埃拉托斯特尼筛法.

//得到自然数n以内的全部素数，依次的所有素数剔除素数的倍数，剩下的就是素数。
//复杂度:O(nloglogn)
int prime[N];  //第i个素数是prime[i],从0开始记录
bool is_prime[N];  //is_prime[i]为true代表i是素数

int sieve(int n)  //筛出小于等于n的所有素数并返回个数
{
    int p=0;  //记录个数
    for(int i=0;i<=n;i++)
        is_prime[i]=true;
    is_prime[0]=false;
    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)
                is_prime[j]=false;
        }
    }
    return p;
}

2.欧拉筛

int prime[N];  //第i个素数是prime[i],从0开始记录
int vis[N];  //vis[i]代表当前合数被筛过没有,保证每个合数只被筛一次
bool is_prime[N];  //is_prime[i]为true代表i是素数

int euler(int n)  //筛出小于等于n的所有素数并返回个数
{
    int cnt=0;  //记录个数
    for(int i=0;i<=n;i++)
    {
        vis[i]=0;
        is_prime[i]=false;
    }
    for(int i=2;i<=n;i++)
    {
        if(vis[i]==0)
        {
            prime[cnt++]=i;
            is_prime[i]=true;
        }
        for(int j=0;j<cnt&&i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
    return cnt;
}

3.米勒拉宾素数测试

ll qpow(ll a,ll b,ll r)//快速幂
{
    ll ans = 1;
    while(b)
    {
        if(b&1)ans=(ans*a)%r;
        a = (a*a)%r;
        b>>=1;
    }
    return ans;
}
bool Miller_Rabbin(int n,int a)//米勒拉宾素数测试
{
    int r=0,s=n-1,j;
    if(!(n%a))
        return false;
    while(!(s&1)){
        s>>=1;
        r++;
    }
    ll k=qpow(a,s,n);
    if(k==1)
        return true;
    for(j=0;j<r;j++,k=k*k%n)
        if(k==n-1)
            return true;
    return false;
}
bool IsPrime(int n)//判断是否是素数
{
    int tab[]={2,3,5,7};
    for(int i=0;i<4;i++)
    {
        if(n==tab[i])
            return true;
        if(!Miller_Rabbin(n,tab[i]))
            return false;
    }
    return true;
}

4.区间筛法

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1000005;
bool is_prime[maxn];
bool is_prime_small[maxn];
ll prime[maxn];
ll prime_num=0;

//对区间[a,b)内的整数执行筛法，is_prime[i-a]=true  ---  表示i是素数 注意这里下标偏移了a，所以从0开始。
void segment_sieve(ll a,ll b) {
    for(ll i=0;i*i<b;++i) is_prime_small[i]=true; //对[2,sqrt(b))的初始化全为质数
    for(ll i=0;i<b-a;++i) is_prime[i]=true; //对下标偏移后的[a,b)进行初始化

    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;  //筛选[2,sqrt(b));
            //(a+i-1)/i得到最接近a的i的倍数，最低是i的2倍，然后筛选
            for(ll j=max(2LL,(a+i-1)/i)*i;j<b;j+=i) is_prime[j-a]=false;
        }
    }
    for(ll i=0;i<b-a;++i)  //统计个数
        if(is_prime[i]) prime[prime_num++]=i+a;
}

int main()
{
    ll a,b;
    while(~scanf("%lld%lld",&a,&b))
    {
        prime_num=0;
        memset(prime,0,sizeof(prime));
        segment_sieve(a,b);
        //for(ll i=0;i<prime_num;++i) printf("%lld\n",prime[i]);
        printf("%lld\n",prime_num);
    }
    return 0;
}