给定两个整数 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.
由于本题数据范围大,所以不能直接线性筛求出所有质数。
但是范围是2^31-1,且r-l小于100万,所以所有合数的最小质因子小于5e5,所以我们可以求出5e5范围内质数,在筛去 l~ r范围内所有合数求出l~r所有质数,再一一枚举求最大最小距离
注意1要筛去
max(k * 2,ll(ceil(l*1.0/k)*k))是求出l-r范围内至少是prime[j]两倍的数,ceil为向上取整函数
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
const int N = 1e6+10;
int p[N],cnt,pr[N],flag;
bool st[N],is_p[N];
ll l,r;
inline void init(){
cnt=0;
int M=N/2;
for(int i(2);i<M;i++){
if(!st[i])p[cnt++]=i;
for(int j=0;p[j]<=M/i;j++){
st[p[j]*i]=1;
if(i%p[j]==0)break;
}
}
}
int main(){
init();
while(~scanf("%lld %lld",&l,&r)){
memset(is_p,0,sizeof is_p);
for(int i(0);i<cnt;i++){
ll k=p[i];
for(ll j=max(k*2,ll(ceil(l*1.0/k)*k));j<=r;j+=k)
is_p[j-l]=1;
}
flag=0;
for(int i=0;i<=r-l;i++)
{
if(!is_p[i]&&l+i>1)pr[flag++]=l+i;
}
if(flag<2)printf("There are no adjacent primes.\n");
else {
int minp=0,maxp=0;
for(int i(0);i<flag-1;i++)
{
if(pr[i+1]-pr[i]<pr[minp+1]-pr[minp])minp=i;
if(pr[i+1]-pr[i]>pr[maxp+1]-pr[maxp])maxp=i;
}
printf("%d,%d are closest, %d,%d are most distant.\n",pr[minp],pr[minp+1],pr[maxp],pr[maxp+1]);
}
}
}