给定两个整数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≤231−1
输入样例:
2 17
14 17
输出样例:
2,3 are closest, 7,11 are most distant.
There are no adjacent primes.
思路:
因为数据范围有点大,时间限制为1s,所以不能开那么大的数组用线性筛法求质数,也不能一个一个的判断是不是素数(会超时),那么有个新思路;
先求出来50000以内所以得质数(用线性筛法);然后直接将50000以内的任意一个质数p扩大到 L — U的范围中;
扩大方法为【L/p】p;方框里为向上取整,也就是等价于(L+p-1)/pp;
然后标记所有的p的倍数即可求出来L——U的质数;
然后对相邻的质数两两比较即可;
(一定要注意过程中的数据范围哪里会爆int,出现段错误!!!!)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxx=1e5+7;
const int N=1e6+7;
int pre[N],cnt;
bool str[N];
void prime(int top)//找到50000以内的素数
{
memset(pre,0,sizeof pre);//以下初始化非常重要!!
memset(str,0,sizeof str);
cnt=0;
for(int i=2;i<=top;i++)
{
if(!str[i]) pre[cnt++]=i;
for(int j=0;pre[j]*i<=top;j++)
{
str[pre[j]*i]=true;
if(i%pre[j]==0) break;
}
}
}
int main()
{
int l,u;
while(cin >>l>>u)
{
prime(50000);
memset(str,0,sizeof str);
for(int i=0;i<cnt;i++)
{
ll p=pre[i];
for(ll j=max((l+p-1)/p*p,2*p);j<=u;j+=p)//通过p的倍数去找属于区间l-u的质数时,临界值会爆int;
{
str[j-l]=true;//巧妙地存储方式;
}
}
cnt=0;//这里的把cnt归零也是个坑;
for(int i=0;i<=u-l;i++)
{
if(!str[i]&&i+l>1) pre[cnt++]=i+l;//把被标记的数找出来,存起来;并且当l=1的时候要特判
}
if(cnt<2) printf("There are no adjacent primes.\n");
else
{
int minn=0,maxx=0;
for(int i=0;i+1<cnt;i++)
{
int d=pre[i+1]-pre[i];
if(d>pre[maxx+1]-pre[maxx]) maxx=i;
if(d<pre[minn+1]-pre[minn]) minn=i;
}
printf("%d,%d are closest, %d,%d are most distant.\n",pre[minn],pre[minn+1],pre[maxx],pre[maxx+1]);
}
}
return 0;
}