#数论，线性筛，eratosthenes#poj 2689 Prime Distance

最新推荐文章于 2023-11-01 17:11:55 发布

原创最新推荐文章于 2023-11-01 17:11:55 发布 · 224 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#poj 2689 #Prime Distance

模拟专栏收录该内容

131 篇文章

订阅专栏

本文介绍了一种高效算法，用于在指定区间内寻找差值最大和最小的相邻质数对。通过试除法和线性筛算法，实现了对大量数字的快速筛选，确保了算法的时间效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目

在一个区间 $[l . . . r]$ ，求两对相邻的质数，一对差值最大，一对最小(可以相同) $(1≤l≤r≤2147483648(unsigned))(1\leq l\leq r\leq 2147483648(unsigned))$

分析

然而直接跑绝对会TLE，但是可以用试除法的思想可以发现可以跑2到 $r\sqrt r$ 筛一筛（我选离线），然后每次询问筛一筛，就没有什么了，时间复杂度 $O(47000+(r−l)log⁡log⁡(r−l))O(47000+(r-l)\log \log(r-l))$

代码

#include <cstdio>
#include <bitset>
typedef unsigned uit;
uit prime[4793],v[46341],m,l,r;
int main(){
	for (register uit i=2;i<=46340;i++){//线性筛
		if (!v[i]) v[i]=i,prime[++m]=i;
		for (register uit j=1;j<=m;j++){
			if (prime[j]>v[i]||prime[j]>46340/i) break;
			v[i*prime[j]]=prime[j];
		}
    }
	while (scanf("%u%u",&l,&r)==2){//unsigned
		if (l==1) l=2; if (r==1) r=2;
	    std::bitset<1000010>v1; v1.reset();
		for (register uit j=1;j<=m;j++)
		for (register uit k=l/prime[j];k<=r/prime[j];k++)
		if (prime[j]*k>=l&&prime[j]*k<=r&&k>1) v1[prime[j]*k-l]=1; //欧式筛
		uit pri2=l,pri1,ans1=0,ans2=r-l<<1,a1,a2,b1,b2;
		while (v1[pri2-l]&&pri2<=r) pri2++;//找到较小的质数
		for (register uit j=pri2;j<=r;j=pri1){
			pri1=j+1;
			while (v1[pri1-l]&&pri1<=r) pri1++;//找较大的质数
			if (pri1<=r&&pri2<=r&&!v1[pri1-l]&&!v1[pri2-l]){
				if (ans1<pri1-pri2) ans1=pri1-pri2,a1=pri2,a2=pri1;
				if (ans2>pri1-pri2) ans2=pri1-pri2,b1=pri2,b2=pri1;
			}
			pri2=pri1;//更新较小的质数
		}
		if (!ans1) printf("There are no adjacent primes.\n");//不存在答案
		else printf("%d,%d are closest, %d,%d are most distant.\n",b1,b2,a1,a2);
	}
	return 0;
}