POJ2689

最新推荐文章于 2020-02-17 18:16:57 发布

AC_王道

最新推荐文章于 2020-02-17 18:16:57 发布

阅读量793

点赞数

分类专栏： acm解题报告文章标签： acm 数论素数

本文链接：https://blog.csdn.net/oqqtzb/article/details/8929418

版权

acm解题报告专栏收录该内容

10 篇文章 0 订阅

订阅专栏

题目连接：http://poj.org/problem?id=2689
转载请注明出处：http://hi.baidu.com/pub/show/createtext
题意：求某个区间内素数最近距离和素数最远距离

思路：二次筛选法
开始根据素数定理可知在最大范围的素数分布在50000的区间内，我们可以求出在这范围内的素数，这是第一次筛选
void doprime1()
{
    numprime1=0;
    memset(isprime,true,sizeof(isprime));
    isprime[1]=false;
    for(long long i=2; i<maxn; i++)
    {
        if(isprime[i])
        {
            prime1[++numprime1]=i;
            for(long long j=i*i; j<maxn; j+=i)
            {
                isprime[j]=false;
            }
        }
    }
}
由于数据较大  我们必须在[l,u]之间在把这个区间的素数筛选出来，通过开始的素数表吧[l,u]的非素数排除掉
void doprime2()
{
    long long b;
    long long j;
    memset(isprime,true,sizeof(isprime));
    for(int i=1; i<=numprime1; i++)
    {
        b=l/prime1[i];
        while(b*prime1[i]<l||b<=1)
            b++;
        for(j=b*prime1[i]; j<=u; j+=prime1[i])
        {
            if(j>=l)
            {
                isprime[j-l]=false;
            }
        }
    }
    if(l==1)
    {
        isprime[0]=0;
    }
AC Code
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<climits>
using namespace std;
const int maxn=50000;
const int maxm=99999999;

bool isprime[maxn*20];
int prime1[maxn];
int numprime1;
void doprime1()
{
    numprime1=0;
    memset(isprime,true,sizeof(isprime));
    isprime[1]=false;
    for(long long i=2; i<maxn; i++)
    {
        if(isprime[i])
        {
            prime1[++numprime1]=i;
            for(long long j=i*i; j<maxn; j+=i)
            {
                isprime[j]=false;
            }
        }
    }
}

long long l,u;
long long prime2[1000001];
int numprime2;
void doprime2()
{
    long long b;
    long long j;
    memset(isprime,true,sizeof(isprime));
    for(int i=1; i<=numprime1; i++)
    {
        b=l/prime1[i];
        while(b*prime1[i]<l||b<=1)
            b++;
        for(j=b*prime1[i]; j<=u; j+=prime1[i])
        {
            if(j>=l)
            {
                isprime[j-l]=false;
            }
        }
    }
    if(l==1)
    {
        isprime[0]=0;
    }
}

void solve()
{
    long long minr,minl,maxr,maxl;
    long long mmin=maxm,mmax=0;
    doprime2();
    numprime2=0;
    for(int i=0; i<=u-l; i++)
    {
        if(isprime[i])
        {
            prime2[++numprime2]=i+l;
        }
    }
    if(numprime2<=1)cout<<"There are no adjacent primes."<<endl;
    else
    {
        for(int i=1; i<numprime2; i++)
        {
            if(prime2[i+1]-prime2[i]>mmax)
            {
                mmax=prime2[i+1]-prime2[i];
                maxl=prime2[i];
                maxr=prime2[i+1];
            }
            if(prime2[i+1]-prime2[i]<mmin)
            {
                mmin=prime2[i+1]-prime2[i];
                minl=prime2[i];
                minr=prime2[i+1];
            }
        }
        cout<<minl<<","<<minr<<" are closest, "<<maxl<<","<<maxr<<" are most distant."<<endl;
    }
}

int main()
{
    doprime1();
    while(cin>>l>>u)
    {
        solve();
    }
    return 0;
}