Prime Distance[POJ2689]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门
POJ传送门(数据更强)

思路

$\mathrm{L},\mathrm{R}$范围很大，所以直接线性筛是过不了的。

但是$\mathrm{L}-\mathrm{R}$的范围只有$10^6$，

对于$[\mathrm{L},\mathrm{R}]$区间内的任意合数$n$，一定存在小于$\sqrt{n}$的约数，证明不难。

所以我们现将$[2,\sqrt{n}]$的素数筛出来，用这些素数去标记$[L,R]$区间内的任意合数，

之后将$[L,R]$中未被标记的数取出，两两比较即可。

细节:

1. 如果用数组存素数的话，需要开long long
2. 注意$L$可能等于$1$,
3. $[2,\sqrt{n}]$素数可能会被标记。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#define ll long long 
using namespace std;
const int N=1e6+10;
const int inf=50010;
bool v[N];int prime[N],m;
void g_p()
{
	m=0;memset(v,false,sizeof(v));
	for(int i=2;i<=inf;i++)
	{
		if(!v[i])prime[++m]=i;
		for(int j=1;j<=m&&i*prime[j]<=inf;i++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int main()
{
	g_p();
	int l,r;
	while(~scanf("%d%d",&l,&r))
	{
		memset(v,true,sizeof(v));
		if(l==1)v[0]=false;
		for(int i=1;i<=m;i++)
		{
			int j=l/prime[i];
			if(l%prime[i])j++;
			j=max(2,j);
			for(;j<=r/prime[i];j++)
			{
				v[j*prime[i]-l]=0;
			}
		}
		int t1=2147483647,t2=0;
		int s1,s2,s3,s4;int cnt=0,last=0;
		for(int i=0;i<=r-l;i++)if(v[i])
		{
			cnt++;
			if(cnt!=1)
			{
				if(i+l-last<t1)t1=i+l-last,s1=last,s2=i+l;
				if(i+l-last>t2)t2=i+l-last,s3=last,s4=i+l;
			}
			last=i+l;
		}
		if(!t2)puts("There are no adjacent primes.");
		else printf("%d,%d are closest, %d,%d are most distant.\n",s1,s2,s3,s4);
	}
	return 0;
}