A - Prime Distance（区间素数筛[a,b]，注意不是[a,b)）

最新推荐文章于 2021-05-22 22:05:17 发布

Forward in time

最新推荐文章于 2021-05-22 22:05:17 发布

阅读量210

点赞数 1

分类专栏：数论

本文链接：https://blog.csdn.net/qq_44555205/article/details/99732701

版权

数论专栏收录该内容

40 篇文章 0 订阅

订阅专栏

在这里插入图片描述

题意：找[L,U]区间里面的递增的素数序列；求相邻差值最小和差值最大的两对素数；如果没有就输出There are no ahjacent primes;
很明显区间素数筛；但是注意一点：这里是两边都是闭区间，白书上是半闭半开；
所以AC代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
typedef long long ll;
const ll maxn=2e6+10;
const ll INF=0x3f3f3f3f;
vector<ll> P;
ll is_prime[maxn],is_small_prime[maxn];
void segment_sieve(ll a,ll b){//区间素数筛 
	  for(ll i=0;i*i<=b;i++)is_small_prime[i]=1;//利用合数+唯一分解定理就可以知道这里为什么这样筛了
	  for(ll i=0;i<=b-a;i++)is_prime[i]=1;
	  for(ll i=2;i*i<=b;i++){
	  	  if(is_small_prime[i])
	  	  {
	  	  	   for(ll j=2*i;j*j<=b;j+=i)is_small_prime[j]=0;//筛
	  	  	   for(ll j=max(2ll,(a+i-1)/i)*i;j<=b;j+=i)is_prime[j-a]=0;//降区间，因为j至少是i的2倍，找到最接近a的第一个i的倍数
			}
	  }
	  for(ll i=0;i<=b-a;i++){//求这个区间素数有哪些 
	  	  if(is_prime[i])P.push_back(i+a);
	  }
}
void Find(ll a,ll b){
	ll Min=INF,Max=0;
	ll f1=0,f2=0;
	 for(ll i=1;i<P.size();i++){
	 	if(P[i]-P[i-1]<Min){
	 		  f1=i;
	 		  Min=P[i]-P[i-1];
		 }
		 if(P[i]-P[i-1]>Max){
		 	Max=P[i]-P[i-1];
		 	f2=i;
		 }
	 }
	 if(P.size()<=1){//如果这个区间只有一个素数 或者没有   
	 	puts("There are no adjacent primes.");
	 }else{
	 	printf("%lld,%lld are closest, %lld,%lld are most distant.\n",P[f1-1],P[f1],P[f2-1],P[f2]);
	 }
}
int main(){
	ll a,b;
	ios_base::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	while(~scanf("%lld %lld",&a,&b)){
		if(a<2)a=2;///这里需要注意因为1不是素数
			segment_sieve(a,b);
			Find(a,b);
			P.clear();
	}
	return 0;
}