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;
}