poj 2689Prime Distance(区间素数）埃氏筛法

这道题的L和R都很大，所以如果直接开一个1~R的数组明显会超时。但是R-L并不大，所以我们考虑把这个区间（L--R）移动到（1--（R-L+1））这个区间再开数组（就是把每个数减L再加1）。接下来先用埃氏筛分（可以自行百度）求出【2，√R】区间的素数，并存在prime数组里。对于prime数组里的每一个质数，求出其在区间（L--R）的倍数并且标记成false（非素数）。那么剩下的区间（L--R）里的数就都是素数咯~然后相邻的比较，求出差最大的和差最小的即可。

注意的细节：1.判断素数的数组（prime是用来储存素数数据的，用long long）isprime[]，和 a[] 最好用bool，节省空间。

　　　　　　2.这道题非bool的数据类型最好都开long long的，不然大的数据点会RE。

　　　　　　3.注意左区间是“1”的情况，“1”是非素数。

下面给出代码：（必要的有注释）

 #include<cstdio>
 #include<vector>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 typedef long long ll;
 ll prime[],cnt;//prime存储【2，√R】区间的素数，cnt记录该区间素数个数
 bool isprime[];//isprime:初始判断 【2，√R】区间的素数
 bool a[];//原区间左右至（1~（R-L+1）） 后，判断素数
 ll l,r;
 void ini()
 {
     memset(isprime,,sizeof(isprime));
     isprime[]=isprime[]=;
     memset(a,,sizeof(a));
 }
 void sushu()//初始用埃氏筛法，筛选 【2，√R】区间的素数
 {
     for (ll i=;i*i<=r;i++)
     {
         if (isprime[i])
         {
             prime[cnt++]=i;
             for (ll j=i*i;j*j<=r;j+=i)
               isprime[j]=;
         }
     }
 }
 void sift()//筛选区间内的素数
 {
     for (ll i=;i<cnt;i++)
     {
         ll bm=l/(ll)prime[i];
         for (ll j=bm*(ll)prime[i];j<=r;j+=(ll)prime[i])
             if ((j)!=(ll)prime[i]) a[j-l+]=;
     }
     if (l==) a[]=;//注意“1”不是素数
 }
 void minus()//相邻素数求差最大和差最小
 {
     ll pre=,minans=,maxans=,x1,x2,y1,y2;
     for (ll i=;i<=(r-l+);i++)
     {
         if (a[i])
         {
             if (pre)
             {
                 if (i-pre > maxans)
                   maxans=i-pre,x1=pre,x2=i;
                 if (i-pre < minans)
                   minans=i-pre,y1=pre,y2=i;
             }
             pre=i;
         }
     }
     if (maxans && minans!=)
       printf("%lld,%lld are closest, %lld,%lld are most distant.\n",(ll)y1+l-,(ll)y2+l-,(ll)x1+l-,(ll)x2+l-);
     else printf("There are no adjacent primes.\n");
 }
 int main()
 {
     while (scanf ("%lld%lld",&l,&r)!=EOF)
     {
         ini();
         sushu();
         sift();
         minus();
     }
     return ;
 }