【阶段1】【198. 质数距离】线性筛+巧妙优化

【题意】
给定两个整数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≤2^31−1 

 【输入样例】
2 17
14 17 

 【输出样例】
2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

分析

• 因为U<=2^31所以我们没法直接求【1，U】之间的素数，但是题目中L与U的差不会超过1000000，给了我们一些提示
• 对于任何数是否为素数，只需要判断它是否能被sqrt(这个数)内的任何数整除即可，进一步优化，就是是否能被sqrt（这个数）内的任何素数整除即可，如果它能被一个合数整除，它必然能被这个合数的质约数整除，而显然这个质约数是包含在这个范围内的。而sqrt（U）的范围较小，我们可以用线性筛选出此范围内的素数，这样我们就将题目优化为：sqrt+线性筛
• 如何证明（对于任何数是否为素数，只需要判断它是否能被sqrt(这个数)内的任何数整除即可）？
• 若a是X的约数，必有X/a是X的约束，如果a>sqrt(X)，必有X/a<sqrt(X)，如果X/a>sqrt(X),必有X<sqrt(X)，得证

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define INF 0xf0f0f0f
using namespace std;
typedef long long ll;
inline ll read()
{
	char c=getchar();ll s=0,f=1;
	while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)){s=s*10+c-'0';c=getchar();}
	return f*s;
}
ll l,r,pri[510000];
bool bk[51000000];
int tot;
inline void pre()
{
	memset(bk,true,sizeof(bk));
	ll tt=sqrt(r+1);
	for(int i=2;i<=tt;i++)
	{
		if(bk[i])pri[++tot]=i;
		for(int j=1;j<=tot && pri[j]*i<=tt;j++)
		{
			bk[i*pri[j]]=false;
			if(i%pri[j]==0)break;
		}
	}
}
int main()
{
	while(scanf("%lld%lld",&l,&r)!=EOF)
	{
		tot=0;
		pre();
		memset(bk,true,sizeof(bk));
		ll last=-1,ans1=INF,t1,t2,ans2=0,t3,t4,ttt=0;
		if(l==1)bk[0]=false;//如有需要，必须将1标记为非素数 
		for(int i=1;i<=tot;i++)
		{
			for(ll j=(l+pri[i]-1)/pri[i];j<=r/pri[i];j++)
			{
				//因为约数一般是成对出现的，我们只需要用小于sqrt(X)的那个约数去标记就可以了
				//如果不这样做就会很有问题，  例如素数被等于自身的素数标记为非素数 
				if(pri[i]>sqrt(j*pri[i]))continue;
				bk[j*pri[i]-l]=false;//这里有巧妙优化，同一减去一个l，避免数太大标记不了 
			}
		}
		for(ll i=l;i<=r;i++)
		{
			if(bk[i-l]==true)
			{
				if(last!=-1)
				{
					if(ans1>i-last)
					{
						ans1=i-last;
						t1=last;t2=i;
					}
					if(ans2<i-last)
					{
						ans2=i-last;
						t3=last;t4=i;
					}
				}
				last=i;ttt++;
			}
		}
		if(ttt>=2)printf("%d,%d are closest, %d,%d are most distant.\n",t1,t2,t3,t4);
		else printf("There are no adjacent primes.\n");
	}
	return 0;
}