Prime Distance(区间筛)

题目链接
题意：输入区间[l,u]，其中l和u为int范围的整数，区间最大为1000000。求出[l,u]中，相邻素数只差最大和最小的素数对。当存在多个时，输出较小的素数对。
思路：
首先看数据量，$l$和$r$是$2^{31}$开不出来这个数组，其次O(n)线性筛也筛不出来，再往后看$r-l<=1e6$这就为我们区间筛做了很大的准备，通过这个条件我们可以想起来把$l->r$转化成$\mathrm{０}－(r-l)$这样数组最大开出$1e6$能开出来,内存的问题解决了，我们来看时间：我们筛$r$的素数需要知道$根r$的素数，那么我们就可以利用根r中的素数来筛$l-r$肯定不能是$1-r$。$l-r$是1e6的长度，$1-r$就是$2^{31}$的长度。那如何来看$l-r$那，我们可以倍增素数的方法来筛。就是对与$2-根r$的p(素数)我们找出$l-r$第一个是p(这个素数)倍数的数，之后就是倍增p(这个素数).直到r。其实在这里我们就能发现这个其实就是区间筛的一种形式，只不过左端点不一样。区间筛左端点是2倍的p(素数)。
找倍数的左右端点

x=(l-1)/prime[i]+1;
///之所以这样，是因为我们要保证是大于l的第一个数
y=r/prime[i];

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1e6+7;
bool isprime[maxn];
ll prime[maxn],p[maxn];
ll cnt;
void Prime(int x) {
    for(int i=2; i<=x; i++) {
        if(!isprime[i]) prime[++cnt]=i;
        for(int j=1; j<=cnt&&i*prime[j]<=x; j++) {
            isprime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int main() {
    Prime(50000);
    ll l,r;
    while(~scanf("%lld%lld",&l,&r)) {
        if(l==1) l++;
        memset(p,0,sizeof p);
        for(int i=1; i<=cnt; i++) {
            ll x=(l-1)/prime[i]+1;///是prime[i]倍数的第一个数。
            ll y=r/prime[i];
            for(int j=x; j<=y; j++) if(j>1)p[j*prime[i]-l]=1; ///是合数　是prime[i]的倍数
        }
        ///l--->r
        ///0--->r-l
        int x=-1,maxx=-1,minn=1e9+7,xi,yi,xm,ym;
        for(int i=0; i<=r-l; i++) {
            if(p[i]==0) {
                if(x==-1) {
                    x=i;
                    continue;
                }
                if(maxx<i-x) {
                    maxx=i-x;
                    xm=x+l;
                    ym=i+l;
                }
                if(minn>i-x) {
                    minn=i-x;
                    xi=x+l;
                    yi=i+l;
                }
                x=i;
            }
        }
        if(maxx==-1)cout<<"There are no adjacent primes."<<endl;
        else cout<<xi<<","<<yi<<" are closest, "<<xm<<","<<ym<<" are most distant."<<endl;
    }
    return 0;
}

才疏学浅，只能写成这样了，如有更清晰明白的，可以分心给俺这个fw.