算法竞赛——进阶指南 POJ2689

L,R取值1-2e9  没有直接的算法求出1-R之间的所有质数。

然我们知道合数x一定是由小于x的一个质数乘上一个数得来。

我们可以算出1-sqrtR 中的质数,然后用这些质数(类比埃筛)筛去L-R之间的合数

复杂度:基本就是埃筛的复杂度

Osqrt(R)*loglogsqrt(r)+(R-L)*loglogR;

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a,b,c,d;
const int M =1e5+7;
ll vs[M],p[M],m;
ll vv[M*10];
bool gao(ll L,ll R)
{
	memset(vv,0,sizeof(vv));
	memset(vs,0,sizeof(vs));
	ll u=sqrt(R);m=0;
	for(ll i=2;i<=u;i++)
	{
		if(vs[i])continue;
		for(int j=i+i;j<=u;j+=i)vs[j]=1;
		p[++m]=i;
		ll x=i*2;
		if(x<L)x=L/i*i;if(x<L)x+=i;
//		cout<<i<<" "<<x<<endl;
		for(ll j=x;j<=R;j+=i)
		{
			ll tp=j-L;
			vv[tp]=1;
		}
	}
	m=0;
	for(ll i=L;i<=R;i++)
	{
		ll tp=i-L;
		if(!vv[tp]&&i!=1)
		{
			p[++m]=i;
		//	cout<<i<<endl;
		}
		
	}
	if(m<2)return false;
	ll mi=p[2]-p[1],mx=p[2]-p[1];
	a=c=p[1],b=d=p[2];
	for(int i=3;i<=m;i++)
	{
		int tp=p[i]-p[i-1];
		if(tp<mi)mi=tp,a=p[i-1],b=p[i];
		if(tp>mx)mx=tp,c=p[i-1],d=p[i];
	}
	return true;
}
int main()
{
	ll L,R;
	while(~scanf("%lld%lld",&L,&R))
	{
		if(gao(L,R))printf("%lld,%lld are closest, %lld,%lld are most distant.\n",a,b,c,d);
		else printf("There are no adjacent primes.\n");
	}
	return 0;
}
发布了284 篇原创文章 · 获赞 13 · 访问量 1万+

猜你喜欢

转载自blog.csdn.net/bjfu170203101/article/details/104362416