Help Hanzo （数论（素数（打表进阶（任意段区域打表）））））

题目来源：https://vjudge.net/problem/LightOJ-1197
【题意】
求a，b区间里的素数个数。
【思路】
第一时间，想到打表，心想到又遇见一道水题，，然后，很自然的去看数据范围，mmp，极限值，傻眼了。。。同时又看到b-a<=100000，一般这样的题，都是从小范围入手，于是就开始了探索。
结果，还是不会写，，于是就翻了博客，看了一半，迷迷糊糊，心想着学长的教诲，要把别人的变成自己的，于是凭着一知半解的感悟自己想，猛然间发现，打表的数字10^5不是随意定下的，因为可以看到他的平方刚好就是极限（差不多。。）。
然后想是不是根据1~10^5的素数表可以推出来想要得到的a，b之间的素数表，然后就发现可以，因为哦，假设a是1000005，那么对于素数2来说，1000006不是素数，1000008等等都不是素数，也就是用1000005/2==x，判断2*x是否为1000005，若是，则1000005不是素数，然后就这样一直往后推，若不是，也可以继续往后推。
但是呢，有一个问题就出现了，众所周知，打表的时候需要一个vis数组进行判断的，但是如果a超级大，开不了那么大的数组怎么办？
然后，就想啊，我们要求的只是a，b这一段的vis，为什么不把他们的数值同时减去a呢？
以下是两种代码：
都差不多啦。。。
【代码1】

#include<set>#include<map>#include<stack>#include<cmath>#include<queue>#include<cstdio>#include<string>#include<cstring>#include<iostream>#include<algorithm>#define mem(a,b) memset(a,b,sizeof(a))using namespace std;const int INF=1e9;typedef unsigned long long ll;typedef long long LL;const int maxn=1e5+10;int prime[maxn];bool vis1[maxn],vis2[maxn];int l=0;void init(){    mem(vis1,0);    vis1[1]=1;    for(int i=2; i<=maxn; i++)    {        if(!vis1[i])        {            prime[l++]=i;            for(int j=2*i; j<=maxn; j+=i)                vis1[j]=1;        }    }}int main(){    init();    int T,cases=1;    scanf("%d",&T);    while(T--)    {        LL a,b,ans=0;        scanf("%lld%lld",&a,&b);        if(b<=maxn-1)//相当于第二种代码里的a=1的特判        {            for(LL i=a; i<=b; i++)            {                if(!vis1[i])                    ans++;            }            printf("**\n");        }        else        {            mem(vis2,0);            for(LL i=0; i<l&&prime[i]<=b; i++)            {                LL x=a/prime[i];                if(x*prime[i]<a)                    x++;                for(LL j=x*prime[i]; j<=b; j+=prime[i])                {                    if(j!=prime[i])                    vis2[j-a]=1;                }            }            for(LL i=a; i<=b; i++)                if(!vis2[i-a])                   ans++;        }        printf("Case %d: %lld\n",cases++,ans);    }}

【代码2】

#include<set>#include<map>#include<stack>#include<cmath>#include<queue>#include<cstdio>#include<string>#include<cstring>#include<iostream>#include<algorithm>#define mem(a,b) memset(a,b,sizeof(a))using namespace std;const int INF=1e9;typedef unsigned long long ll;typedef long long LL;const int maxn=1e5+10;int prime[maxn];bool vis1[maxn],vis2[maxn];int l=0;void init(){    mem(vis1,0);    vis1[1]=1;    for(int i=2; i<=maxn; i++)    {        if(!vis1[i])        {            prime[l++]=i;            for(int j=2*i; j<=maxn; j+=i)                vis1[j]=1;        }    }}int main(){    init();    int T,cases=1;    scanf("%d",&T);    while(T--)    {        LL a,b,ans=0;        scanf("%lld%lld",&a,&b);        mem(vis2,0);        for(LL i=0; i<l&&prime[i]<=b; i++)        {            LL x=a/prime[i];            if(x*prime[i]<a)                x++;            for(LL j=x*prime[i]; j<=b; j+=prime[i])                if(j!=prime[i])                    vis2[j-a]=1;        }        for(LL i=a; i<=b; i++)            if(!vis2[i-a])                ans++;        if(a==1)            ans-=1;        printf("Case %d: %lld\n",cases++,ans);    }}