给定两个整数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.
思路:l和r的范围是 2^31−1,这个范围超出了线性筛素数的范围,因为线性素数筛选法的时间复杂度为o(n),2^31−1大概为2e9,所以把这个范围内的素数全筛出来超出了时间和空间的限制。
但是,对于l~r这个区间的长度最大为1e6,显然这个范围是可以的,这个区间内的每个合数x,都满足性质:能找到一个小于sqrt(x)的一个质因子,而最大的x为2^31−1,而sqrt(x)最多50000就够了,所以只要找出50000以内的所有素数,然后把l~r这个区间内的所有50000以内的所有素数的倍数都筛掉,剩下的就是这个区间内的所有素数了,在保存这些素数时,由于数组下标可能会越界,所以整体偏移l个长度即可;
最后用两指针分别找距离最近和最远的两对素数的较小者的下标即可。
完整代码:
#include <iostream>
#include <algorithm>
#include <cstring>
#define int long long
using namespace std;
const int maxn=1e6+5;
int primes[maxn],cnt;
bool vis[maxn];
void get_primes(int n)
{
memset(vis,false,sizeof vis);
cnt=0;
for(int i=2;i<=n;i++){
if(!vis[i]) primes[cnt++]=i;
for(int j=0;primes[j]*i<=n;j++){
vis[primes[j]*i]=true;
if(i%primes[j]==0) break;
}
}
}
signed main()
{
int l,r;
while(cin>>l>>r){
get_primes(50000);
memset(vis,false,sizeof vis);
//用50000内的所有素数将[l,r]内的合数筛掉:
for(int i=0;i<cnt;i++){
int p=primes[i];
for(int j=max(2*p,(l+p-1)/p*p);j<=r;j+=p){
vis[j-l]=true;//数向左偏移l
}
}
//用数组存储[l,r]内所有的素数:
cnt=0;
for(int i=0;i<=r-l;i++){//所有数向左偏移了l,所以从0遍历r-l个即为[l,r]全部r-l个数
if(!vis[i]&&i+l>1) primes[cnt++]=i+l;//i+l向右偏移还原成原来的素数
}
if(cnt<2) cout<<"There are no adjacent primes."<<endl;
else{
int minp=0,maxp=0;
for(int i=0;i+1<cnt;i++){
if(primes[i+1]-primes[i]<primes[minp+1]-primes[minp]) minp=i;
if(primes[i+1]-primes[i]>primes[maxp+1]-primes[maxp]) maxp=i;
}
cout<<primes[minp]<<","<<primes[minp+1]<<" are closest, "<<primes[maxp]<<","<<primes[maxp+1]<<" are most distant."<<endl;
}
}
return 0;
}