poj2689Prime Distance 素数筛选

最新推荐文章于 2022-07-26 14:51:29 发布

laserss

最新推荐文章于 2022-07-26 14:51:29 发布

阅读量846

点赞数 1

分类专栏： acm pku 数学

本文链接：https://blog.csdn.net/dellaserss/article/details/8746161

版权

acm 同时被 3 个专栏收录

111 篇文章 0 订阅

订阅专栏

pku

25 篇文章 0 订阅

订阅专栏

数学

16 篇文章 0 订阅

订阅专栏

poj2689Prime Distance

题目要求：在给定范围L~R中找出差值最大和最小的两组素数

要在L~R筛选素数，直接算出1~2,147,483,647的话肯定超时

但实际上，筛选1~2,147,483,647中的素数需要的因子并不多，只要找出2~sqrt(2,147,483,647)中的素数就够了，筛选法就是这么干的

之后用这些素数对L~R的区间试除得到L~R内的素数求出结果

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

bool Q[1001000];
bool p[50000];
int P[5000];

int main()
{
	__int64 i,j,k,L,R;
	for(i=2;i<=50000;i+=2)	p[i]=0,p[i+1]=1;	p[1]=0;p[2]=1;P[1]=2;P[0]=1;
	for(i=1;i<=46345;i+=2)
	{
		if(p[i])
		{
			int k=i*i;
			while(k<50000){
				p[k]=0;
				k+=2*i;
			}
			P[++P[0]]=i;
		}
	}
	
	while(scanf("%I64d%I64d",&L,&R)!=EOF)
	{
		
		int Start;
		for(i=L;i<=R;i++)	Q[i-L]=1;
		for(i=1;i<=P[0];i++)
		{	
			if(P[i]>=R)	break;
			Start=L/P[i];
			if(Start*P[i]<L)	Start++;
			if(Start==1)	Start++;
			for(j=Start;j*P[i]<=R;j++)			Q[j*P[i]-L]=0;
		}
		if(L==1)	Q[0]=0,Q[2]=1;
		if(L==2)	Q[0]=1;
		__int64 pre=R,now=0,min=1000000000,max=0,a,b,x,y;
		for(i=L;i<=R;i++){
			if(Q[i-L]){
				pre=i;
				break;
			}
		}
		for(i=pre+1;i<=R;i++)
		{
			if(Q[i-L]){
				now=i;
				if(now-pre<min){
					min=now-pre;
					a=pre,b=now;
				}
				if(now-pre>max){
					max=now-pre;
					x=pre,y=now;
				}
				pre=now;
			}
		}
		if(!now)	printf("There are no adjacent primes.\n");
		else	printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",a,b,x,y);
	}
	return 0;
}
/*

2000000000 2000000001
2146483648 2147483647

*/