质数距离

题目:

题解:
转化为筛出$[l,r]$区间内的质数

因为对于每一个数$n$,都存在一个小于等于$\sqrt{n}$的约数

首先求出$\sqrt{2^{23}-1}$ $\approx$ $50000$,所以首先筛出$50000$以内的素数,

然后对于每一个质数$pi$,找到第一个大于等于$l$的倍数$p$,并经行标记,注意最小是$pi$的两倍,

之后没有被标记的点即为区间内的质数,最后循环遍历质数找到最大最小差值即可

#include <bits/stdc++.h>
using namespace std;
long long pri[1000005];
long long st[1000005];
long long primes[1000005];
int con;
void prime()
{
     for(int i=1;i<=50005;i++) pri[i]=1;
     pri[1]=0;
     for(int i=2;i<=50005;i++)
     {
         if(pri[i])
         {
             primes[++con]=i;
             for(int j=2;j*i<=50005;j++)
             {
                 pri[j*i]=0;
             }
         }
     }
}
int main()
{
    long long l,r;
    while(cin>>l>>r)
    {
        prime();
        memset(st,0,sizeof(st));
        for(int i=1;i<=con;i++)
        {
            long long p=primes[i];
            for(long long j=max(p*2,(p+l-1)/p*p);j<=r;j+=p)
            {
                st[j-l]=1;
            }
        }
        int c=0;
        long long cha=r-l;
        //cout<<cha<<endl;
        for(long long i=0;i<=cha;i++)
        {
            if(!st[i]&&i+l>=2) primes[++c]=i+l;
            //cout<<c<<endl;
        }
        //cout<<c<<endl;
        if(c<2) cout<<"There are no adjacent primes."<<endl;
        else 
        {
            long long ma=-1,mi=1e12;
            long long a1,a2;
            for(long long i=1;i<=c-1;i++)
            {
                long long d=primes[i+1]-primes[i];
                if(d>ma) 
                {
                    ma=d;
                    a1=i;
                }
                if(d<mi)
                {
                    mi=d;
                    a2=i;
                }
            }
            cout<<primes[a2]<<","<<primes[a2+1]<<" are closest, ";
            cout<<primes[a1]<<","<<primes[a1+1]<<" are most distant."<<endl;
        }
    }
    return 0;
}