2018.09.01 poj2689 Prime Distance（埃式筛法）

传送门

一道挺有趣的。

第一眼以为每个数都用miller_rabin判一次，但感觉会被卡时间啊。

继续分析发现可以晒出sqrt(r)中的所有素数，然后用类似埃式筛法的方法晒出[l,r]" role="presentation" style="position: relative;">[l,r][l,r]之间的素数。

然后就没了。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 1000005
#define ll long long
using namespace std;
bool f[N];
ll pri[150005],vis[150005],tot=0;
inline void init(ll len){
    for(ll i=2;i<=len;++i){
        if(!vis[i])pri[++tot]=i;
        for(ll j=1;j<=tot;++j){
            ll k=pri[j]*i;
            if(k>len)break;
            vis[k]=1;
            if(i%pri[j]==0)break;
        }
    }
}
inline void check(ll l,ll r){
    memset(f,1,sizeof(f));
    if(l==1)f[0]=0;
    for(ll i=1;pri[i]*pri[i]<=r;++i){
        ll tmp=pri[i];
        for(ll j=(l/tmp)*tmp;j<=r;j+=tmp){
            if(j<l||j==tmp)continue;
            f[j-l]=0;
        }
    }
}
int main(){
    init(50000);
    ll l,r;
    while(~scanf("%lld%lld",&l,&r)){
        check(l,r);
        ll cnt=0,tmp,mx=0,mn=1e9,pos,mx1,mx2,mn1,mn2;
        for(ll i=0;i<=r-l;++i){
            if(f[i]){
                ++cnt;
                if(cnt>=2){
                    tmp=i-pos;
                    if(tmp>mx)mx=tmp,mx1=pos,mx2=i;
                    if(tmp<mn)mn=tmp,mn1=pos,mn2=i;
                }
                pos=i;
            }
        }
        if(cnt>=2)printf("%lld,%lld are closest, %lld,%lld are most distant.\n",mn1+l,mn2+l,mx1+l,mx2+l);
        else printf("There are no adjacent primes.\n");
    }
    return 0;
}