UVA10140 Prime Distance

https://www.luogu.com.cn/problem/UVA10140

题意翻译

• 给定两个正整数 l,rl,r，求 [l,r][l,r] 间 相邻 的两个差最大的质数和 相邻 的两个差最小的质数。如果区间内质数个数 \le 1≤1，输出 There are no adjacent primes.。
• 1\le l<r\le 2^{31}1≤l<r≤2^31，r-l\le 10^6r−l≤106。

输入输出样例

输入 #1复制

2 17
14 17

输出 #1复制

2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

---

这个r的上界是2^32>=1e9，所以线性筛也没办法预处理。（当然我知道大佬们还有其他的一个神奇筛法可以过）但是我觉得做题还是锻炼为主。

我们看这个题，2^32次没法预处理，那么往r-l的范围在1e6去思考。这里就要知道一个结论了：任何一个合数n必定包含至少一个不超过sprt(n)的质因子

我的理解如下：

任何一个数n有因子的话那么因子在sqrt(n)两端任意分布

若n是合数，那么n有因子同样在sqrt（n）两端任意分布，按照算术基本定理（任何一个大于1的正整数都能唯一分解成有限个质数的乘积）

那么，n为合数的因子也满足sqrt（n）两端对称分布且能质因数分解，那么至少有一个<=sqrt(n)的质因子。

那么也就是说：对于一个数n，一定存在<=sqrt(n)的因子，而这个因子能质因数分解的话，满足任何一个合数n必定包含至少一个不超过sprt(n)的质因子，如果不能质因数分解，那么这个数就是质数，也满足任何一个合数n必定包含至少一个不超过sprt(n)的质因子。

那么数学书上的解释：有没有好心人帮忙解释解释..或者等我哪天想不开了学初等数论解释

回到这道题。对于询问的区间[L,R]里面的数只有质数和合数两种，我们知道了里面的任意一个数n，如果是质数，就是所求；如果是合数，那么必定包含至少一个不超过sprt(n)的质因子 。既然这样，我们预处理[2,sqrt(2^34)]的质数，然后通过把【L，R】里面能被这些质因子整除的给标记掉，剩下的就是质数了。最后再O（n）判断就好了。

#include<iostream>
#include<vector>
#include<queue>
#include<map> 
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=1e6+10;
typedef long long LL;
map<LL,LL> b;//标记[l,r]中被质数p整除的数 
LL c[maxn];//[l,r]中剩下的合数 
//埃式筛
LL v[maxn];//合数标记 
LL a[maxn];//质数集存[1,sqrt(2^31)]的质数 
LL k=0;
LL x1,y3,x2,y2;//最后找相邻质数 查了下y1和c++内置变量名冲突... 
void primes(LL n)
{
	memset(v,0,sizeof(v));
	for(LL i=2;i<=n;i++)
	{
		if(v[i]) continue;
		a[++k]=i;
		for(LL j=i;j<=n/i;j++) v[i*j]=1;
	}
}
//
int main(void)
{
  cin.tie(0);std::ios::sync_with_stdio(false);
  primes(sqrt(2147483648));
  LL l,r;
  while(cin>>l>>r)
  {
  		b.clear();
  		LL m=0;
  		if(l==1) l=2;//处理 '1' 
  		for(LL i=1;i<=k;i++)
  		{
			for(LL j=l/a[i];j<=r/a[i];j++)
			{
				if(a[i]*j<l||a[i]*j>r) continue;
				if(j>1)	b[a[i]*j]=1;
			}
		}
		for(LL i=l;i<=r;i++)
			if(!b[i]) c[++m]=i;
		if(m<=1) 
		{
			cout<<"There are no adjacent primes."<<endl;
			continue;
		}	 
		LL t1=2147483647;LL t2=0;
		for(LL i=1;i<=m-1;i++)
		{
			if(t2<c[i+1]-c[i]){t2=c[i+1]-c[i];x1=c[i];y3=c[i+1];}
			if(t1>c[i+1]-c[i]){t1=c[i+1]-c[i];x2=c[i];y2=c[i+1];}
		}
		cout<<x2<<','<<y2<<" are closest, "<<x1<<','<<y3<<" are most distant."<<endl;
  }
return 0;
}