【Ybtoj 第28章例2】质数距离【质数和约数】

---

解题思路

首先我们发现：$R-L$ 的范围很小，我们应该要能够快速求出 $L\sim R$ 之间的质数。

显然有推论：任意一个合数 x 必定包含一个不超过 $\sqrt{x}$的质因子。

所以我们可以筛出 $[1,\sqrt{R}]$ 之间的所有质数，对于每个质数 pp，把 $[L,R]$中能被 $p$ 整除的数标记为合数。最终没有被标记的数就是质数，对相邻的质数两两比较，找出差值最小和最大的即可。

---

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<iomanip>
#include<cstring>
#include<cmath>
#include<map>
#include<queue>
#define ll long long
#define ldb long double
using namespace std;

const int N=50010;
ll l,r,cnt,pri[10000010],v[10000010],ans[10000010];

void solve(){
	if(l==1)l++;
	ll s=0;
	for(ll i=1; i<=cnt; i++) {
		for(ll j=l/pri[i]; j<=r/pri[i]; j++)
			if(j>1)
				v[j*pri[i]-l]=1;
			}
	for(int ll i=l;i<=r;i++)
		if(!v[i-l])
			ans[++s]=i;
			
	ll x1,x2,y1,y2,maxn=0,minn=2147482600;
	for(int i=1;i<s;i++)
	{
		ll x=ans[i+1]-ans[i];
		if(x>maxn)maxn=x,x1=ans[i],x2=ans[i+1];
		if(x<minn)minn=x,y1=ans[i],y2=ans[i+1];
	}
	if(!maxn)printf("There are no adjacent primes.\n");
	else printf("%d,%d are closest, %d,%d are most distant.\n",y1,y2,x1,x2);
}

int main() {
	for(ll i=2; i<=N; i++) {
		if(!v[i]) {
			v[i]=i;
			pri[++cnt]=i;
		}
		for(ll j=1; j<=cnt; j++) {
			if(pri[j]>v[i]||pri[j]>N/i)	break;
			v[i*pri[j]]=pri[j];
		}
	}
	while(scanf("%lld%lld",&l,&r)!=EOF) {
		memset(v,0,sizeof(v));
		solve();
	}
}