POJ 2689 Prime Distance

点击打开链接

题意：

给你一个区间。求这个区间内相邻最近的两个素数和相邻最远的两个素数分别是多少。

题解：

很明显直接找出l，u之间的素数，然后遍历一遍。就能求出结果。

但这个题l，u，范围比较大。直接筛法打表肯定不能实现，

又已知   u-l < 1e6  所以单独求出l，u之间的素数就行。即筛去非质数。

已知1e32内的要么是素数，要么是能被 sqrt（1e32）整除的数，即 l，u中的所有非素数的因子都在sqrt（1e32）之中。

所以将sqr（1e32）内的素数先筛出来。然后用这些数筛去l，u之间的非素数。

 

#include <cmath>
#include <stdio.h>
#include <cassert>
#include <algorithm>
#include <cstring>
#define inf 99999999
typedef long long LL;
const int maxn = 50000;
int tot,cnt, pr[maxn], d[maxn];
bool notprime[maxn*20]={1,1};
LL l,u,prime2[maxn*20+1],i,j,k;
void makeprime2(){
    LL a,b;
    memset(notprime,0,sizeof(notprime));
    for(i=0;i<tot;++i){
        a=(l-1)/pr[i]+1;
        b=u/pr[i];
        for(j=a;j<=b;++j)
            if(j>1) notprime[j*pr[i]-l]=1;
    }
}
int main(){
    tot=0;
	for(i=2; i<maxn; ++i){//快速线性筛法求素数
		if(!notprime[i])
			pr[tot++] = i;
		for(j=0,k;j<tot&&(k=i*pr[j])<maxn;++j){
			notprime[k]=1;
			if(!(i%pr[j]))
				break;
		}
	}
	while(~scanf("%lld%lld",&l,&u)){
         if(l==1) l++;
         makeprime2();
         cnt=0;

         LL maxp=-inf,minp=inf,minl,minr,maxl,maxr;
         for(i=0;i<=u-l;++i)
            if(!notprime[i]){
                prime2[cnt++]=i+l;
            }
//for(i=0;i<cnt;i++) printf("%d ",prime2[i]);printf("\n");
        if(cnt<2)
            printf("There are no adjacent primes.\n");
        else {
            for(i=0;i<cnt-1;++i){
                if(prime2[i+1]-prime2[i]<minp){
                    minp=prime2[i+1]-prime2[i];
                    minl=prime2[i],minr=prime2[i+1];
                }
                if(prime2[i+1]-prime2[i]>maxp){
                    maxp=prime2[i+1]-prime2[i];
                    maxl=prime2[i],maxr=prime2[i+1];
                }
            }
            printf("%lld,%lld are closest, %lld,%lld are most distant.\n",minl,minr,maxl,maxr);
        }
    }
	return 0;
}