POJ - 2689 Prime Distance (区间筛)

题意：求[L,R]中差值最小和最大的相邻素数(区间长度不超过1e6)。

由于非素数$n$必然能被一个不超过$\sqrt n$的素数筛掉，因此首先筛出$[1,\sqrt R]$中的全部素数，然后用这些素数去筛$[L,R]$。

注意1不是素数，但不会被筛掉，需要特殊处理(巨坑~~)

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll N=1e5+10,inf=0x3f3f3f3f;
ll L,R,isp1[N],isp2[N*10],pri[N*10],npri;
void getpri1() {
    for(ll i=2; i<N; ++i)isp1[i]=1;
    for(ll i=2; i<N; ++i)if(isp1[i])for(ll j=2; i*j<N; ++j)isp1[i*j]=0;
}
void getpri2() {
    for(ll i=0; i<=R-L; ++i)isp2[i]=1;
    for(ll i=2; i*i<=R; ++i)if(isp1[i])
            for(ll j=max(2ll,(L-1)/i+1); i*j<=R; ++j)isp2[i*j-L]=0;
}
int main() {
    getpri1();
    while(scanf("%lld%lld",&L,&R)==2) {
        getpri2();
        npri=0;
        for(ll i=L; i<=R; ++i)if(isp2[i-L]&&i!=1)pri[npri++]=i;
        if(npri<2)puts("There are no adjacent primes.");
        else {
            ll mi=1,mx=1;
            for(ll i=2; i<npri; ++i) {
                if(pri[i]-pri[i-1]<pri[mi]-pri[mi-1])mi=i;
                if(pri[i]-pri[i-1]>pri[mx]-pri[mx-1])mx=i;
            }
            printf("%lld,%lld are closest, %lld,%lld are most distant.\n",pri[mi-1],pri[mi],pri[mx-1],pri[mx]);
        }
    }
    return 0;
}