1246.素数筛选

http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1246
Time Limit: 2000 MS Memory Limit: 65536 K
Total Submissions: 749 (119 users) Accepted: 89 (58 users)

Description

请求出区间[A,B]之间的素数的个数。

Input

输入有多组数据。
第一行为正整数T(T<10),表示输入的数据组数。
第二行到第T+1行，每行2个数字A,B(0<A<=B<2^31,且B-A<2^20)表示区间的范围。

Output

输出T行数字，每行只有一个数字，表示对应的素数的个数。

Sample Input

2
100 200
2 1000000

Sample Output

21
78498

#include <iostream>
#include <cstring>
using namespace std;
const int N = 46341;
bool isprime[1048576+10];
unsigned int prime[5000];
int Prime()
{
int i,j,k=0;
memset(isprime,1,sizeof(isprime));
for(i=2;i<=N;i++)
{
if(isprime[i])
{
prime[k++]=i;
for(j=2*i;j<=N;j+=i)
isprime[j]=0;
}
}
return k;
}
int main ()
{
unsigned int j,t,n,m,ans,i,k,p;
p=Prime();
cin>>t;
while (t--)
{
cin>>n>>m;
if(n<2) n=2;
ans=m-n+1;
memset(isprime,1,sizeof(isprime));
for(i=0;i<p&&prime[i]*2<=m;i++)
{
k=n/prime[i];
if(n%prime[i]) k++;
if(k==1) k++;
for(j=k*prime[i];j<=m;j+=prime[i])
{
if(isprime[j-n])
{
ans--;
isprime[j-n]=0;
}
}
}
cout<<ans<<endl;
}
}