poj 2689 Prime Distance

题意 给两个数m,n求m~n间的任意相邻的两个素数的最短距离和最大距离，若不存在两个素数输出

There are no adjacent primes.

这个题坑的我WA了十多次，千万注意数据范围：m-n<1000000;
然而m和n均在int范围内所以必须采用映射比如说数m~n之间的任意一个数i都可以用i-m存到数组里，若求的是2000000~3000000，则int数组b;b[0]记录的就是2000000，b[1]记录的就是2000001；

解题的大致思路：因为两数之差不超过一千万，所以在m~n之间的素因子必定在根号214783648<50000之下，所以一开始只要把小于50000的素数筛选出来即可，再通过这些素数标记m~n之间的合数，剩下的就是这个区间里的素数了。

代码：

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
long long zd,zx,a[50100],b[1000010],c[50100];
int main(){
    long long i,j,k,m,n,end1,end2,sum=0,end3,end4,u;
    long long p=floor(sqrt(50010));
    for(i=2;i<=p;i++){
        if(a[i])continue;
        for(j=i*i;j<=50010;j+=i){
            a[j]=1;
        }
    }
    for(i=2;i<=50010;i++)if(!a[i])c[++sum]=i;
    while(scanf("%lld%lld",&m,&n)!=EOF){
        if(m==1)m=2;
        long long zd=0;
        long long zx=999999999;
        long long cnt=-1;
        memset(b,0,sizeof(b));
        for(i=1;i<=sum;i++)
        {
            long long l,r;
            if(m%c[i]==0)l=m/c[i];
            else l=m/c[i]+1;
            r=n/c[i];
            for(j=l;j<=r;j++)
                if(j>1)b[-1*m+c[i]*j]=1;
        }
        for(i=m;i<=n;i++)
            if(b[i-m]==0)b[++cnt]=i;
        if(cnt<1)printf("There are no adjacent primes.\n");
        else {
            for(i=0;i<cnt;i++){
                if(zd<b[i+1]-b[i]){
                    zd=b[i+1]-b[i];
                    end1=b[i];
                    end2=b[i+1];
                }
                if(zx>b[i+1]-b[i]){
                    zx=b[i+1]-b[i];
                    end3=b[i];
                    end4=b[i+1];
                }
            }
            printf("%lld,%lld are closest, %lld,%lld are most distant.\n",end3,end4,end1,end2);
        }
    }
    return 0;
}