lightoj1197Help Hanzo(二次筛)

来源：互联网发布：淘宝单号充值编辑：程序博客网时间：2024/05/17 06:10

原题链接：http://lightoj.com/volume_showproblem.php?problem=1197
题目大意：查询一个区间[a,b]内的质数个数,其中1≤a≤b<231，但是有a−b≤1e5，那么首先我想到用miller-rabin冲一下，结论是数据卡了miller-rabin，反正我冲不过去，那就先打出来一个小范围的素数表用筛法去筛就行了。

#include <stdio.h>#include<iostream>#include<cmath>#include<cstring>using namespace std;typedef long long LL;LL pri[100010];LL vis[100010];int cnt=0;void init(){    for(int i=2;i<=100000;i++){        if(!vis[i])pri[++cnt]=(LL)i;        for(int j=1;j<=cnt&&pri[j]*i<=100000;j++){            vis[pri[j]*i]=1;            if(i%pri[j]==0)break;        }    }}int main(void){    int t,cas=0;//  freopen("1197.in","r",stdin);//  freopen("1197.out","w",stdout);    init();    scanf("%d",&t);    while(t--){        LL a,b;        scanf("%lld%lld",&a,&b);        if(a==1)a++;        LL mx=b-a;        for(LL i=0;i<=mx;i++)vis[i]=0;        for(int i=1;i<=cnt&&pri[i]<b;i++){            LL k=(a-1)/pri[i]+1;            if(k==1)k++;            for(LL j=k*pri[i];j<=b;j+=pri[i]){                vis[j-a]=1;            }        }        LL ans=0;        for(LL i=0;i<=mx;i++)if(!vis[i])ans++;        printf("Case %d: %lld\n",++cas,ans);    }    return 0;}

0 0