POJ2689Prime Distance 素数筛选+枚举

最新推荐文章于 2024-03-13 21:53:20 发布

Qianzhez。

最新推荐文章于 2024-03-13 21:53:20 发布

阅读量109

点赞数

分类专栏： poj

本文链接：https://blog.csdn.net/qianzhezhi1127/article/details/97023696

版权

poj 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

题目大意

给一个左右区间L,U（1<=L< U<=2,147,483,647），但是U-L<1000000，找出两对连续的素数，一对的差最大（相距最远），一对的差最小（相距最近），没有就输出There are no adjacent primes.有就输出。
样例
Sample Input2 17
14 17
Sample Output2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

解题思路

1.因为给的范围很大，不可能全部本地存下来，所以就存根号那个最大的，大概5w的样子,离线得到前五万个质数
2.然后根据所给的区间的到要枚举的区间（（l-1)/p+1 , u/p )的数
3.离散化把素数挑出来
4.然后就开始枚举挑出来的素数，找出满足条件的

解题代码

#include <iostream>
#include <cstdio>
#include <string.h>
#include <algorithm>
#define INF 0x7fffffff
#define n 50004
#define m 1000001
#define ll long long
using namespace std;
ll p[m],coun=0,l,u;
bool st[m*20];
void is_prime()
{
     memset(st,false,sizeof(st));
    for(ll i=2;i<=n;i++)
{

    if(!st[i])
        p[coun++]=i;
    for(ll j=0;j<coun&&i*p[j]<=n;j++)
    {
        st[i*p[j]]=true;
        if(i%p[j]==0)
            break;
    }

}

}
ll p2[m];
void is_prime2()
{
    if(l<2)
        l=2;
  memset(st,false,sizeof(st));
  for(ll i=0;i<coun;i++)
  {
      ll r1=(l-1)/p[i]+1;
      ll r2=u/p[i];
      for(ll j=r1;j<=r2;j++)
      {
          if(j>1)
            st[p[i]*j-l]=true;
      }

  }



}
int main()
{
    is_prime();

    while(~scanf("%lld%lld",&l,&u))
    {
        is_prime2();


        ll maxn=-INF,minn=INF,coun2=0;
        ll maxnl,maxnu,minnl,minnu;
        for(ll i=0;i<=u-l;i++)
        {
            if(!st[i])
                p2[coun2++]=l+i;


        }

        if(coun2<=1)
            printf("There are no adjacent primes.\n");
        else
        {
            for(int i=0;i<coun2-1;i++)
        {
            if(p2[i+1]-p2[i]>maxn)
            {
                maxnl=p2[i];
                maxnu=p2[i+1];
                maxn=p2[i+1]-p2[i];
            }
            if(p2[i+1]-p2[i]<minn)
            {
                minnl=p2[i];
                minnu=p2[i+1];
                minn=p2[i+1]-p2[i];
            }
        }
        printf("%lld,%lld are closest, %lld,%lld are most distant.\n",minnl,minnu,maxnl,maxnu);
        }


    }
    return 0;
}