传送门
一道挺有趣的。
第一眼以为每个数都用miller_rabin判一次,但感觉会被卡时间啊。
继续分析发现可以晒出sqrt(r)中的所有素数,然后用类似埃式筛法的方法晒出[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;
}