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;

}