LightOJ 1197 区间素数筛(注意溢出的处理)

点击打开链接

#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
typedef long long LL;
const int maxn=1e5+10;
const int maxc=65536;
bool p[maxn+10];
int cnt=0;
int prime[maxc];
void is_prime()
{
    memset(p,true,sizeof(p));
    p[0]=p[1]=false;
    for(int i=2;i<=maxc;i++)
        if(p[i]) {
           prime[cnt++]=i;
           for(LL j=(LL)i*i;j<=maxc;j+=i)                   //注意开long long避免溢出
                p[j]=false;
       }
}
int main()
{
    LL n,m;
    int T,kase=0;
    is_prime();
    scanf("%d",&T);
    while(T--)
     {
            cin>>n>>m;
            int ans=0;
            memset(p,true,sizeof(p));
            int k=floor(sqrt(m)+0.5);
            int p1=lower_bound(prime,prime+cnt,k)-prime;
            for(int i=0;i<=p1;i++)
            {
                LL st=max((LL)n/prime[i]*prime[i],(LL)prime[i]*prime[i]);
                for(LL j=st;j<=m;j+=prime[i])
                    if(j>=n) p[j-n]=false;
            }
            for(int i=0;i<=m-n;i++)
                if(p[i]) ans++;
            if(n==1) ans--;
            printf("Case %d: %d\n",++kase,ans);
     }
    return 0;
}