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;

    }

}