Help Hanzo

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/A7_RIPPER/article/details/81912593

题目大意：

求区间a,b内素数的数量

解题思路：

由于b极大，所以打表会爆内存。但并不意味着放弃打表，我们可以先打一个小点的素数表出来，如果b在这个表内直接二分找一下a,b就可以了。否则利用到b-a<=100000这个性质，可以开一个这么大的桶下标表示为j-a来筛选a-b内的素数，这样就用到我们之前的小素数表来筛选了。
代码如下：

#include<iostream>
#include<cstdio>
#include<fstream>
#include<set>
#include<cmath>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<iomanip>
#include<cstdlib>
#include<list>
#include<queue>
#include<stack>
#include<algorithm>
#define inf 0x3f3f3f3f
#define MOD 1000000007
#define mem0(a) memset(a,0,sizeof(a))
#define mem1(a) memset(a,-1,sizeof(a))
#define meminf(a) memset(a,inf,sizeof(a))
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
#define MAXN 1000010
bool vis1[MAXN],vis2[MAXN];
int prime[MAXN];
int getprime()//欧拉筛一个小的素数表
{
  int cnt=0;
  for(int i=2;i<MAXN;i++)
  {
    if(vis1[i]==0)
    {
      prime[cnt++]=i;
      vis1[i]=1;
    }
    for(int j=0;j<cnt&&prime[j]*i<MAXN;j++)
    {
      vis1[i*prime[j]]=1;
      if(i%prime[j]==0)
      break;
    }
  }
  return cnt;
}
int main()
{
  std::ios::sync_with_stdio(false);
  cin.tie(0);
  //freopen("test.txt","r",stdin);
  //freopen("output.txt","w",stdout);
  int len=getprime();
  int t,cas=1;;
  cin>>t;
  while(t--)
  {
    ll a,b,result=0;
    cin>>a>>b;
    if(b<=MAXN)//b在素数表内
    {
      int loab=lower_bound(prime,prime+len,b)-prime;//二分找a和b
      if(prime[loab]>b)loab-=1;
      int loaa=lower_bound(prime,prime+len,a)-prime;
      result=loab-loaa+1;//答案
    }
    else//不在素数表内
    {
        memset(vis2,0,sizeof(vis2));
        for(int i=0;i<len;i++)//用已经得到的素数表筛去a-b内的其他合数
        {
            ll k=a/prime[i];
            if(k*prime[i]<a) k++;
            for(ll j=k*prime[i];j<=b;j+=prime[i]) vis2[j-a]=1;
        }
        for(ll i=a;i<=b;i++)
        {
            if(!vis2[i-a])
                result++;//统计区间内素数得数量
        }
    }
    cout<<"Case "<<cas++<<": "<<result<<endl;
  }
}