POJ 2689 Prime Distance （区间大数素数）

题目链接：http://poj.org/problem?id=2689

题目大意：给两个数L，R，问【L，R】之间的相邻两个素数最大相差多少，最少相差多少。1<= L< R<=1e9      1<=R - L<=1e6

样例：

2 17

2-17素数有 2，3，5，7，11，13，17 。 2和3相差最小，7和11相差最大为4，13和17也为4，题目要求取前边的。

 

解题思路：

        注意到L，R的范围很大所以任何筛都不可能在规定时间筛出所有的数，但是注意到R-L的范围只有1e6 ，所以可以讲素数的范围偏移也就是说数组ans[0]代表了L这个数是不是素数，0则是，1则不是，ans【1】代表了L+1这个数.....ans【r-l】代表了r这个数；

       那么问题来了怎么判断ans数组代表的数字是不是素数，首先我们要明确一个事实就是任何一个合数x他都有至少一个小于的因子，所以我们可以筛出2~sqrt（r）的全部素数，sqrt（r）的素数可以快速筛选出，然后他们的倍数都是合数，也就是说 j*p要被筛掉 （p为枚举到的素数）就是L/p<=j<=R/p，和埃氏筛一个思想，不过这里要注意到的是，这里不能像埃氏筛一样直接莽j*pri【i】因为他可能超出了数组的范围，所以一定要加判断如果j*pri【i】大于了r就要跳出循环不然会RE到自闭，ans[j*pri[i]-L]=1;这样就解决了数据过大，偏移后下标和真实代表的数据差L。

          这个题还有一个坑点就是，L=1的时候要ans【0】=1；因为1不是质数，偏移之后他不会被筛掉，还有就是如果这个枚举到的素数刚刚好在这个区间，那就比较骚了，特殊判断一下就好了，这里的L/p<=j<=R/p,要注意L/p要向上取整，因为一旦是向下取整，那么j*pri【i】-L就有可能小于0 ，从而又无情的RE。

笔者不才，代码极其丑陋。。。。各位可以按照思路打。代码就别看了

 

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
#define ll long long
using namespace std;
ll ans[1000005];
ll pri[50005];
bool is_prime[50005];
int main()
{
    ll num=0;
    memset(is_prime,0,sizeof(is_prime));
    ll r,l;
    for(ll i=2; i<=50000; i++)
    {
        if(is_prime[i])continue;
        pri[num++]=i;
        for(ll j=i; j<=50000/i; j++)
        {
            is_prime[i*j]=1;
        }
    }
    while(cin>>l>>r)
    {
        memset(ans,0,sizeof(ans));
       if(l==1)ans[0]=1;
        for(ll i=0;i<num;i++){
                bool ff=0;
            if(pri[i]>sqrt(r))break;
            ll left=ceil(l*1.0/pri[i]);
            ll right=r/pri[i];
            if(right==0)break;
            for(ll j=left;j<=right;j++){
                if((j*pri[i])>(r))break;
                if(j==left&&left==1)continue;
                ans[j*pri[i]-l]=1;
            }
        }

        ll aa=-1,bb=-1,Min_dis=0x7fffffff,Max_dis=-1;
        for(ll i=0;i<=r-l;i++)
        {
            if(!ans[i]&&aa!=-1&&bb!=-1){aa=bb;bb=i;}
            if(!ans[i]&&aa==-1){aa=i;continue;}
            if(!ans[i]&&bb==-1)bb=i;
            if(bb!=-1&&aa!=-1&&bb-aa>Max_dis)Max_dis=bb-aa;
            if(bb!=-1&&aa!=-1&&bb-aa<Min_dis)Min_dis=bb-aa;
        }
        ll ans_a=-1,ans_b=-1,ans_c=-1,ans_d=-1;
        bool f1=0,f2=0;
        aa=-1,bb=-1;
        for(ll i=0;i<=r-l;i++)
        {
            if(!ans[i]&&aa!=-1&&bb!=-1){aa=bb;bb=i;}
            if(!ans[i]&&aa==-1){aa=i;continue;}
            if(!ans[i]&&bb==-1)bb=i;
            if(bb!=-1&&aa!=-1&&bb-aa==Max_dis&&!f2){ans_a=aa+l;ans_b=bb+l;f2=1;}
            if(bb!=-1&&aa!=-1&&bb-aa==Min_dis&&!f1){ans_c=aa+l;ans_d=bb+l;f1=1;}
        }
        if(ans_a!=-1&&ans_b!=-1&&ans_c!=-1&&ans_d!=-1){
            printf("%lld,%lld are closest, %lld,%lld are most distant.\n",ans_c,ans_d,ans_a,ans_b);
        }
        else printf("There are no adjacent primes.\n");
    }
    return 0;
}