埃筛(区间筛)

来源:互联网 发布:erp系统模拟软件 编辑:程序博客网 时间:2024/06/14 21:10

题目>>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);}
原创粉丝点击