L,R取值1-2e9 没有直接的算法求出1-R之间的所有质数。
然我们知道合数x一定是由小于x的一个质数乘上一个数得来。
我们可以算出1-sqrtR 中的质数,然后用这些质数(类比埃筛)筛去L-R之间的合数
复杂度:基本就是埃筛的复杂度
Osqrt(R)*loglogsqrt(r)+(R-L)*loglogR;
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a,b,c,d;
const int M =1e5+7;
ll vs[M],p[M],m;
ll vv[M*10];
bool gao(ll L,ll R)
{
memset(vv,0,sizeof(vv));
memset(vs,0,sizeof(vs));
ll u=sqrt(R);m=0;
for(ll i=2;i<=u;i++)
{
if(vs[i])continue;
for(int j=i+i;j<=u;j+=i)vs[j]=1;
p[++m]=i;
ll x=i*2;
if(x<L)x=L/i*i;if(x<L)x+=i;
// cout<<i<<" "<<x<<endl;
for(ll j=x;j<=R;j+=i)
{
ll tp=j-L;
vv[tp]=1;
}
}
m=0;
for(ll i=L;i<=R;i++)
{
ll tp=i-L;
if(!vv[tp]&&i!=1)
{
p[++m]=i;
// cout<<i<<endl;
}
}
if(m<2)return false;
ll mi=p[2]-p[1],mx=p[2]-p[1];
a=c=p[1],b=d=p[2];
for(int i=3;i<=m;i++)
{
int tp=p[i]-p[i-1];
if(tp<mi)mi=tp,a=p[i-1],b=p[i];
if(tp>mx)mx=tp,c=p[i-1],d=p[i];
}
return true;
}
int main()
{
ll L,R;
while(~scanf("%lld%lld",&L,&R))
{
if(gao(L,R))printf("%lld,%lld are closest, %lld,%lld are most distant.\n",a,b,c,d);
else printf("There are no adjacent primes.\n");
}
return 0;
}