质数距离（二次筛选）

给定两个整数 L 和 U，你需要在闭区间 [L,U] 内找到距离最接近的两个相邻质数 C1 和 C2（即 C2−C1 是最小的），如果存在相同距离的其他相邻质数对，则输出第一对。
同时，你还需要找到距离最远的两个相邻质数 D1 和 D2（即 D1−D2 是最大的），如果存在相同距离的其他相邻质数对，则输出第一对。
输入格式
每行输入两个整数 L 和 U，其中 L 和 U 的差值不会超过 106。
输出格式
对于每个 L 和 U，输出一个结果，结果占一行。
结果包括距离最近的相邻质数对和距离最远的相邻质数对。（具体格式参照样例）
如果 L 和 U 之间不存在质数对，则输出 There are no adjacent primes.。
数据范围
1≤L<U≤231−1
输入样例：
2 17
14 17
输出样例：
2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

主要思想就是筛选法求素数,因为n,m过于大，需要筛选两次
解题步骤：
1、先用一个数组存下一定范围内的素数
2、把n-m之间的所有非素数都做标记（因为数太大了，需要把 j 放到 第 i-n 这个位置）
值得学习的是：

在 l - r 之间的第一个p的倍数是（ l / p) * p ,但 （ l / p） 是下取整 , 所以再往后一个 (l + p - 1) / p * p
但是第1个 x 的倍数 (l + p - 1) / p * p 可能包含了它本身，本身当然不能被删去 , 所以用 2 * x 来限制
3、然后把素数再存入数组中
（2，3两步就是1中预处理的步骤）
4、最后求所有距离，然后比较

#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdio>
#include <iomanip>

using namespace std;
typedef long long ull;
const int N=1000001;
ull n,t,k,sum=0,p,h,m,i,j;
ull ge;
int a[N];
int b[N],c[N],d[N]; //d数组存开始的素数

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    for(i=2; i*i<=N; i++) //预处理
    {
        for(j=i*i; j<=N; j+=i)
        {
            a[j]=1; //非素数
        }
    }
    ge=0;
    for(i=2; i<=N; i++)
    {
        if(a[i]==0) d[ge++]=i;//存下所有素数
    }
    
    while(cin>>n>>m)
    {
        memset(b,0,sizeof(b));
        for(i=0; i<ge; i++) //把区间内的合数用他们的最小质因子筛掉
        {
            p=d[i];
            for(j=max(2*p,(n+p-1)/p*p); j<=m; j+=p)
            {
                b[j-n]=1;
            }
        }
        k=0;
        for(i=0; i<=m-n; i++)
        {
            if(b[i]==0 && i+n>1)// i+n>1特殊处理1
            {
                c[k++]=i+n; //原数存上
            }
        }
        if(k<2) cout<<"There are no adjacent primes."<<endl;
        else
        {
            int maxx=0;
            int minn=0;
            for(int i=0; i<k-1; i++)
            {
                int d = c[i+1]-c[i]; //两个数之间的距离
                //cout<<c[i+1]<<' ' <<c[i]<<endl;
                //cout<<d<<endl;
                if(d < c[minn+1] - c[minn]) minn = i;//因为要记录是哪个数，就记下他的坐标
                if(d > c[maxx+1] - c[maxx]) maxx = i;

            }
            cout<<c[minn]<<","<<c[minn+1]<<" are closest, "<<c[maxx]<<","<<c[maxx+1]<<" are most distant."<<endl;
        }

    }
    return 0;
}