0
点赞
收藏
分享

微信扫一扫

Prime Distance(区间筛)


​​题目链接​​​ 题意:输入区间[l,u],其中l和u为int范围的整数,区间最大为1000000。求出[l,u]中,相邻素数只差最大和最小的素数对。当存在多个时,输出较小的素数对。
思路:
首先看数据量,开不出来这个数组,其次O(n)线性筛也筛不出来,再往后看这就为我们区间筛做了很大的准备,通过这个条件我们可以想起来把转化成这样数组最大开出能开出来,内存的问题解决了,我们来看时间:我们筛的素数需要知道的素数,那么我们就可以利用根r中的素数来筛肯定不能是是1e6的长度,就是的长度。那如何来看那,我们可以倍增素数的方法来筛。就是对与的p(素数)我们找出第一个是p(这个素数)倍数的数,之后就是倍增p(这个素数).直到r。其实在这里我们就能发现这个其实就是区间筛的一种形式,只不过左端点不一样。区间筛左端点是2倍的p(素数)。
找倍数的左右端点

x=(l-1)/prime[i]+1;
///之所以这样,是因为我们要保证是大于l的第一个数
y=r/prime[i];

#include<bits/stdc++.h>
using namespace std;
#define
const int maxn=1e6+7;
bool isprime[maxn];
ll prime[maxn],p[maxn];
ll cnt;
void Prime(int x) {
for(int i=2; i<=x; i++) {
if(!isprime[i]) prime[++cnt]=i;
for(int j=1; j<=cnt&&i*prime[j]<=x; j++) {
isprime[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
int main() {
Prime(50000);
ll l,r;
while(~scanf("%lld%lld",&l,&r)) {
if(l==1) l++;
memset(p,0,sizeof p);
for(int i=1; i<=cnt; i++) {
ll x=(l-1)/prime[i]+1;///是prime[i]倍数的第一个数。
ll y=r/prime[i];
for(int j=x; j<=y; j++) if(j>1)p[j*prime[i]-l]=1; ///是合数 是prime[i]的倍数
}
///l--->r
///0--->r-l
int x=-1,maxx=-1,minn=1e9+7,xi,yi,xm,ym;
for(int i=0; i<=r-l; i++) {
if(p[i]==0) {
if(x==-1) {
x=i;
continue;
}
if(maxx<i-x) {
maxx=i-x;
xm=x+l;
ym=i+l;
}
if(minn>i-x) {
minn=i-x;
xi=x+l;
yi=i+l;
}
x=i;
}
}
if(maxx==-1)cout<<"There are no adjacent primes."<<endl;
else cout<<xi<<","<<yi<<" are closest, "<<xm<<","<<ym<<" are most distant."<<endl;
}
return 0;
}

才疏学浅,只能写成这样了,如有更清晰明白的,可以分心给俺这个fw.


举报

相关推荐

0 条评论