POJ 2689 Prime Distance

最新推荐文章于 2020-02-02 16:02:40 发布

RaAlGhul

最新推荐文章于 2020-02-02 16:02:40 发布

阅读量407

点赞数

分类专栏： ACM_Water

本文链接：https://blog.csdn.net/RaAlGhul/article/details/51789280

版权

ACM_Water 专栏收录该内容

19 篇文章 0 订阅

订阅专栏

题目链接：http://poj.org/problem?id=2689

一道看起来很简单的题目，就因为没考虑边界的特殊情况，真的是把我坑死了。

首先我们在考虑的时候如果L为1的话，那么我们应该把L+1，因为这个都要WA哭了，因为我们预处理的质数从2开始，那么如果边界里面包含了1，我们后面的处理方式会把1也算成是质数TT。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
typedef long long LL;
using namespace std;
const int maxn = 1000005;
LL prime[maxn>>1],isprime[maxn],fprime[maxn],len;
void init()
{
	for(LL i=2; i<47000; i++)
		if(!isprime[i])
		{
			prime[len++] = i;
			for(LL j=2; i*j<=47000; j++)
				isprime[i*j] = 1;
		}
}
int main()
{
	init();
	LL L,R;
	while(scanf("%I64d%I64d", &L, &R) !=EOF)
	{
		if(L == 1) L++;
		memset(isprime,0,sizeof(isprime));
		for(LL i=0; i<len && prime[i]*prime[i] <= R; i++)
		{
			LL tmp = L/prime[i]+(L%prime[i]>0);
			if(tmp == 1) tmp++;
			for(LL j=tmp*prime[i]; j<=R; j+=prime[i])
				isprime[j-L] = 1;
		}
		int count = 0;
		for(LL i=L; i<=R; i++)
			if(isprime[i-L] == 0)
				fprime[count++] = i;
		if(count < 2) printf("There are no adjacent primes.\n");
		else
		{
			LL M1=fprime[0],M2=fprime[1],m1=fprime[0],m2=fprime[1];
			for(LL i=2; i<count; i++)
			{
				if(fprime[i]-fprime[i-1] < M2-M1) M1=fprime[i-1],M2=fprime[i];
				if(fprime[i]-fprime[i-1] > m2-m1) m1=fprime[i-1],m2=fprime[i];
			}
			printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",M1,M2,m1,m2);
		}
	}
}