埃筛（区间筛）

题目>>http://codevs.cn/problem/3223/

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int MAXN = 3e6;
ll p[MAXN],l,r,ans = 0,cnt = 0;
bool used[MAXN],used2[MAXN];
int main()
{
    scanf("%lld%lld",&l,&r);
    for(ll i = 2;i <= sqrt(r);i ++) //先筛出2~根号r的素数 
        if(!used[i])
        {
            p[++cnt] = i;
            for(ll j = i * i;j <= sqrt(r);j += i)//根据唯一分解定理所有的数都是由质数构成的。 
                used[j] = 1;                    //对于每一个数只要被一个因数素数筛到那他就会被筛掉。 
        }
    for(int i = 1;i <= cnt;i ++)
    {
        if(l % p[i] == 0&&l != p[i]) //特判开头 
            used2[1] = 1;
        for(ll j = l / p[i] * p[i] + p[i];j <= r;j += p[i])//区间筛：先找到该素数所能筛到的第一个数。 
            if(j != p[i])//根号r可能远大于l，因此之前的素数不能计入答案。 
                used2[j - l + 1] = 1;
    }
    for(ll i = l;i <= r;i ++)
        if(!used2[i - l + 1])
            ans ++;
    printf("%lld",ans);
}