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+,inf=0x3f3f3f3f;
 ll L,R,isp1[N],isp2[N*],pri[N*],npri;
 void getpri1() {
     for(ll i=; i<N; ++i)isp1[i]=;
     for(ll i=; i<N; ++i)if(isp1[i])for(ll j=; i*j<N; ++j)isp1[i*j]=;
 }
 void getpri2() {
     for(ll i=; i<=R-L; ++i)isp2[i]=;
     for(ll i=; i*i<=R; ++i)if(isp1[i])
             for(ll j=max(2ll,(L-)/i+); i*j<=R; ++j)isp2[i*j-L]=;
 }
 int main() {
     getpri1();
     while(scanf("%lld%lld",&L,&R)==) {
         getpri2();
         npri=;
         for(ll i=L; i<=R; ++i)if(isp2[i-L]&&i!=)pri[npri++]=i;
         if(npri<)puts("There are no adjacent primes.");
         else {
             ll mi=,mx=;
             for(ll i=; i<npri; ++i) {
                 if(pri[i]-pri[i-]<pri[mi]-pri[mi-])mi=i;
                 if(pri[i]-pri[i-]>pri[mx]-pri[mx-])mx=i;
             }
             printf("%lld,%lld are closest, %lld,%lld are most distant.\n",pri[mi-],pri[mi],pri[mx-],pri[mx]);
         }
     }
     return ;
 }