LightOJ 1197（大数范围筛区间内素数个数）

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cassert>
#include <iostream>
#include <string>
#include <vector>
using namespace std;
typedef long long ll;

const int MAXN = 1000000+10;

bool is_prime[MAXN];//true的为合数，false的为素数
int prime[MAXN/10];//存素数
bool num[MAXN];
int tol;
void find_prime()
{
memset(is_prime,false,sizeof(is_prime));
tol=0;
is_prime[1]=1;
for(int i=2; i<MAXN; i++)
{
if(!is_prime[i])
{
prime[tol++]=i;
for(int j=i*2; j<MAXN; j+=i)
is_prime[j]=1;
}
}
}

int main()
{
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
int T;
cin>>T;
find_prime();
for(int cas=1; cas<=T; cas++)
{
int a,b;
scanf("%d %d",&a,&b);
int n=b-a;
memset(num,false,sizeof(num));
for(int i=0;i<tol&&prime[i]*prime[i]<=b;i++)
{
int j=0;
if(a%prime[i]!=0)
j=prime[i]-a%prime[i];//找到第一个需要筛掉的数(j+a)%prime[i]==0
if(a<=prime[i])
j+=prime[i];//(j+a)/prime[i]==1,则（j+a）是素数，向下推一个
for(;j<=n;j+=prime[i])
num[j]=true;
}
int ans=0;
for(int j=0;j<=n;j++)
{
if(!num[j]) ans++;
}
if(a==1) ans--;
printf("Case %d: %d\n",cas,ans);
}
return 0;
}