一本通1619【例 1】Prime Distance

1619： 【例 1】Prime Distance

题目描述

原题来自：Waterloo local，题面详见 POJ 2689

给定两个整数 L,R，求闭区间 [L,R] 中相邻两个质数差值最小的数对与差值最大的数对。当存在多个时，输出靠前的素数对。

输入格式

多组数据。每行两个数 L,R。

输出格式

详见输出样例。

样例

样例输入

2 17
14 17

样例输出

2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

数据范围与提示

对于全部数据 1≤L<R<2³¹,R−L≤10⁶。

sol：不会算复杂度，总之怎么暴力都能过

先筛出1~500000的质数，在 l~r 中暴力去掉合数（真的非常非常非常非常非常非常暴力）

for(i=;i<=*Prim&&Prim[i]<=r;i++)
{
    int aa=(l-)/Prim[i]+;
    aa+=((aa==)?():());
    int bb=r/Prim[i];
    for(j=aa;j<=bb;j++)
    {
        Bo[j*Prim[i]-l+]|=;
    }
}

然后尴尬的过了

Ps：注意如果 l =1 的话要把 l 变成 2，因为 1 不是质数也筛不掉，不直接去掉会当成答案记录

#include <bits/stdc++.h>
using namespace std;
typedef int ll;
inline ll read()
{
    ll s=;
    bool f=;
    char ch=' ';
    while(!isdigit(ch))
    {
        f|=(ch=='-');
        ch=getchar();
    }
    while(isdigit(ch))
    {
        s=(s<<)+(s<<)+(ch^);
        ch=getchar();
    }
    return (f)?(-s):(s);
}
#define R(x) x=read()
inline void write(ll x)
{
    if(x<)
    {
        putchar('-');
        x=-x;
    }
    if(x<)
    {
        putchar(x+'');
        return;
    }
    write(x/);
    putchar((x%)+'');
    return;
}
inline void writeln(ll x)
{
    write(x);
    putchar('\n');
    return;
}
#define W(x) write(x),putchar(' ')
#define Wl(x) writeln(x)
const int N=;
int Prim[N];
bool Bo[N];
bool Used[N];
inline void Pre_Prime()
{
    int i,j;
    for(i=;i<=;i++)
    {
        if(!Bo[i])
        {
            Prim[++*Prim]=i;
        }
        for(j=;j<=*Prim&&Prim[j]*i<=;j++)
        {
            Bo[Prim[j]*i]=;
            if(i%Prim[j]==) break;
        }
    }
}
int Num[N];
int main()
{
    int i,j,l,r;
    Pre_Prime();
    while(~scanf("%d%d",&l,&r))
    {
        l+=(l==)?():();
        *Num=;
        memset(Bo,,sizeof Bo);
        for(i=;i<=*Prim&&Prim[i]<=r;i++)
        {
            int aa=(l-)/Prim[i]+;
            aa+=((aa==)?():());
            int bb=r/Prim[i];
            for(j=aa;j<=bb;j++)
            {
                Bo[j*Prim[i]-l+]|=;
            }
        }
        for(i=;i<=r-l+;i++) if(!Bo[i])
        {
//            printf("i=%d\n",i);
            Num[++*Num]=i+l-;
        }
        if(*Num<)
        {
            puts("There are no adjacent primes.");
        }
        else
        {
            int O1=,O2=;
            for(i=;i<*Num;i++)
            {
                if(Num[i+]-Num[i]<Num[O1+]-Num[O1]) O1=i;
                if(Num[i+]-Num[i]>Num[O2+]-Num[O2]) O2=i;
            }
            printf("%d,%d are closest, %d,%d are most distant.",Num[O1],Num[O1+],Num[O2],Num[O2+]);
            putchar('\n');
        }
    }
    return ;
}
/*
input
2 17
14 17
output
2,3 are closest, 7,11 are most distant.
There are no adjacent primes.
*/