zoj 1842 Prime Distance

// 数论题，增强的筛法，回想素数筛法
// 只要筛到最大数的开方，剩下的就是素数
// 于是这里，开一个 sqrt(2^31) 大约 65536 的素数表，然后
// 对于每个 L~U 的区间，筛掉 sqrt(U) 之内的素数的倍数即可
#include <iostream>
#include <math.h>
#include <map>
#include <stack>
#include <queue>
#include <vector>
#include <algorithm>
#include <stdio.h>
#include <string.h>
using namespace std;
#define MOD 1000000007
#define maxm 10010
#define maxn 50010
int prim[maxn/3],p;
bool f[maxn];
int gcd(int a,int b){
int r;
while(r=a%b){a=b;b=r;}
return b;
}
bool isp(int n){
if(n==2) return true;
if(n%2==0||n==1) return false;
int m=(int)(sqrt(n+1.0));
for(int i=3;i<=m;i+=2)
if(n%i==0) return false;
return true;
}
void getprime(){
int i,j;
f[1]=true;
for(i=4;i<=maxn;i+=2)
f[i]=true;
int m=(int)(sqrt(maxn+1.0));
for(i=3;i<=m;i+=2){
for(j=i*i;j<=maxn;j+=i)
f[j]=true;
}
for(i=1;i<=maxn;i++)
if(!f[i]) prim[p++]=i;
}
int sum[maxn];
void sum_divisor(int n){
int i,j;
for(i=2;i<=n/2;i++)
for(j=i+i;j<=n;j+=i)
sum[j]+=i;
sum[1]=-1;
}
int v[maxn*20];
int L,U,len;
void shaixuan(){
int i=0;
long long j;
long long l=L,r=U;
while(prim[i]<47000){
j=(l-1+prim[i])/prim[i];
// printf("%lld ",j);
if(j==1) j++;
for(j=prim[i]*j;j<=r;j+=prim[i])
v[j-l]=1;//,printf("%lld ",j);
i++;
}
}
int main()
{
int n;
int m;
int i,k;
int l,r;
//  int tp=~(1<<31);
//  m=sqrt(tp+1.0);
//  printf("%d",m);
getprime();
int maxl,maxr,maxlen,minl,minr,minlen;
while(scanf("%d %d",&L,&U)!=EOF){
len=U-L;

for(i=0;i<=len;i++) v[i]=0;
shaixuan();
l=-1;
maxlen=0;
minlen=MOD;
if(L==1) v[0]=1;
for(i=0;i<=len;i++)
if(!v[i]){// printf("%d ",i+L);
if(l==-1)l=i;
else{
if(i-l>maxlen){
maxlen=i-l;
maxl=l+L;
maxr=i+L;
}
if(i-l<minlen){
minlen=i-l;
minl=l+L;
minr=i+L;
}
l=i;
}
}
if(maxlen){
printf("%d,%d are closest, %d,%d are most distant.\n",minl,minr,maxl,maxr);
}else printf("There are no adjacent primes.\n");
}

return 0;
}