《算法竞赛进阶指南 》质数距离

解题思路：

由于这个题给出的l和u的取值范围非常大达到了1e9的数量级，我们不可能从1到2^31进行线性筛来进行筛选素数，但是题目告诉我们L和U的差值不会超过1e6，这里有一个性质，如果一个数x是合数的话，那他就一定能分解出一个小于根号x的质因数，对2的31次方开根号大概是小于50000的级别，所以可以算出1到50000之间的质数，然后用这些质数去排除那些在给定区间之内的合数，那么剩下的数就是区间之内的质数了，所以现在的关键是如何用1到50000内的质数去筛选掉处于区间内的合数，这里有一个公式用来求大于等于L的最小的一个p（1到50000之间的一个质数）的倍数q=(l+p-1)/p*p，这里的(l+p-1)/p表示l/p上取整,然后我们就只需要从q开始枚举p的的倍数并标记一下，这样我们就可以将区间内的质数和合数进行区分，最后根据题意选出距离最近和最远的2个相邻质数对。

代码如下：

#include <iostream>
#include <cstring>
#include<cstdio>
#include <algorithm>
using namespace std;
const int N =1E6+10;
int primes[N],cnt;
bool st[N];
void get_primes(int n)
{
    memset(st,false,sizeof st);//不要忘记初始化
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        primes[cnt++]=i;
        for(int j=0;i*primes[j]<=n;j++)
        {
            st[i*primes[j]]=true;
            if(i%primes[j]==0)
            break;
        }
    }
}
int main()
{
    long long l,r;
    while(cin>>l>>r)
    {
        cnt=0;       //不要忘记初始化
        get_primes(50010);
        memset(st,false,sizeof st);
        for(int i=0;i<cnt;i++)
        {
            long long p=primes[i];
            for(long long j=max((l+p-1)/p*p,2*p);j<=r;j+=p)//从大于等于l的最小的一个p的倍数开始，注意j最小要从2*p开始
            st[j-l]=true;
        }
        cnt=0;
        for(int i=0;i<=r-l;i++)
        if(!st[i]&&i+l>1)//特判当i+l==1时不应该被视为质数
        primes[cnt++]=i+l;
        if(cnt<2)
        cout<<"There are no adjacent primes."<<endl;
        else
        {
            int maxp=0,minp=0;
            for(int i=0;i+1<cnt;i++)
            {
                int d=primes[i+1]-primes[i];
                if(d<primes[minp+1]-primes[minp])
                minp=i;
                if(d>primes[maxp+1]-primes[maxp])
                maxp=i;
            }
            printf("%d,%d are closest, %d,%d are most distant.\n",primes[minp],primes[minp+1],primes[maxp],primes[maxp+1]);
        }
    }
    return 0;
}