Help Hanzo LightOJ - 1197（素数筛法）

 Help Hanzo LightOJ - 1197

题目链接：

https://cn.vjudge.net/problem/26909/origin

题意：

求两个数a，b之间有多少个素数。

思路：

因为这两个数太大了，用任何办法都无法一次全部筛选出1到b的素数，所以对于这个题，我们只能够先筛选出比较小的一部分素数，再用这一部分的素数来筛选a，b之间的素数。

代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+20;
ll prime[maxn], psize;
//int factors[100][2], fsize;
bool isprime[maxn];
//long long MAX=1LL<<31+5;
//bitset<maxn>isprime;
void initPrimes()
{
    memset(isprime, true, sizeof(isprime));
    //isprime.set();
    isprime[0]=isprime[1]=0;
    for (int i=2; i<=maxn; i++)
    {
        if(isprime[i]) prime[psize++]=i;
        for (int j=0; j<psize && i*prime[j]<=maxn; j++)
        {
            isprime[prime[j]*i]=0;
            if (i%prime[j]==0) break;
        }
    }
}
//bitset<(maxn+100)/5>isprime2;
bool isprime2[100005];
int solve(long long a,long long b)
{
    //memset(isprime2,0,sizeof(isprime2));
    int res=0;
    if(a<=maxn-20&&b<=maxn-20)
    {
        for(int i=a; i<=b; i++)
        {
            if(isprime[i])
                res++;
        }
        return res;
    }
    for(int i=0;i<psize&&prime[i]*prime[i]<=b;i++)
    {
        for(ll j=a/prime[i]*prime[i];j<=b;j+=prime[i])
        {
            if(j>=a&&j>prime[i])
                isprime2[j-a]=1 ;
        }
    }
    for(int i=0;i<=b-a;i++)
    {
        if(!isprime2[i])
            res++;
        else
            isprime2[i]=0;
    }
    return res;
}
int main()
{
    initPrimes();
    int t;
    scanf("%d",&t);
    for(int cas=1; cas<=t; cas++)
    {
        long long a,b;
        scanf("%lld%lld",&a,&b);
        printf("Case %d: %d\n",cas,solve(a,b));
    }
}