质数距离（poj2689——PrimeDistance）

给定两个整数L和U，你需要在闭区间[L,U]内找到距离最接近的两个相邻质数C1和C2（即C2-C1是最小的），如果存在相同距离的其他相邻质数对，则输出第一对。

同时，你还需要找到距离最远的两个相邻质数D1和D2（即D1-D2是最大的），如果存在相同距离的其他相邻质数对，则输出第一对。

输入格式
每行输入两个整数L和U，其中L和U的差值不会超过1000000。

输出格式
对于每个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.

思路：
因为数据范围有点大，时间限制为1s，所以不能开那么大的数组用线性筛法求质数，也不能一个一个的判断是不是素数（会超时），那么有个新思路；
先求出来50000以内所以得质数（用线性筛法）；然后直接将50000以内的任意一个质数p扩大到 L — U的范围中；
扩大方法为【L/p】p；方框里为向上取整，也就是等价于（L+p-1）/pp;
然后标记所有的p的倍数即可求出来L——U的质数；

然后对相邻的质数两两比较即可；
（一定要注意过程中的数据范围哪里会爆int，出现段错误！！！！）

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int maxx=1e5+7;
const int N=1e6+7;
int pre[N],cnt;
bool str[N];
void prime(int top)//找到50000以内的素数
{
    memset(pre,0,sizeof pre);//以下初始化非常重要！！
    memset(str,0,sizeof str);
    cnt=0;
    for(int i=2;i<=top;i++)
    {
        if(!str[i]) pre[cnt++]=i;
        for(int j=0;pre[j]*i<=top;j++)
        {
            str[pre[j]*i]=true;
            if(i%pre[j]==0) break;
        }
    }
}
int main()
{
    int l,u;
    
    while(cin >>l>>u)
    {
        prime(50000);
        memset(str,0,sizeof str);
        for(int i=0;i<cnt;i++)
        {
            ll p=pre[i];
            for(ll j=max((l+p-1)/p*p,2*p);j<=u;j+=p)//通过p的倍数去找属于区间l-u的质数时，临界值会爆int；
            {
                str[j-l]=true;//巧妙地存储方式；
            }
        }
        cnt=0;//这里的把cnt归零也是个坑；
        for(int i=0;i<=u-l;i++)
        {
            if(!str[i]&&i+l>1) pre[cnt++]=i+l;//把被标记的数找出来，存起来；并且当l=1的时候要特判
        }
        if(cnt<2) printf("There are no adjacent primes.\n");
        else
        {
            int minn=0,maxx=0;
            for(int i=0;i+1<cnt;i++)
            {
                int d=pre[i+1]-pre[i];
                if(d>pre[maxx+1]-pre[maxx]) maxx=i;
                if(d<pre[minn+1]-pre[minn]) minn=i;
            }
            printf("%d,%d are closest, %d,%d are most distant.\n",pre[minn],pre[minn+1],pre[maxx],pre[maxx+1]);
        }
    }
    return 0;
}