poj 2689 Prime Distance(大区间素数)

题目链接：poj 2689 Prime Distance

题意：

给你一个很大的区间（区间差不超过100w），让你找出这个区间的相邻最大和最小的两对素数

题解：

正向去找这个区间的素数会超时，我们考虑逆向思维：

我们先用线性筛 筛出前50000的素数，在int范围内的区间的合数的因子都在我们之前筛出来了，

然后我们就把这个区间的合数标记出来，剩下的就是素数了。

标记合数也是用全部的素数去筛一下，这样就能快速的标记这个区间的合数，然后在暴力枚举一下这个区间，更新一下答案。

#include<cstdio>
#include<algorithm>
#include<cmath>
#define F(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
typedef long long ll;

const int N=5e4+7;
int primes[N],tot=0;
bool vis[N];
void Euler(){
    F(i,2,N-7){
        if(!vis[i])primes[++tot]=i;
        F(j,1,tot){
            if(i*primes[j]>N)break;
            vis[i*primes[j]]=1;
            if(i%primes[j]==0)break;
        }
    }
}

int L,U;
int a[N*20],ed;

int main()
{
    Euler();
    while(~scanf("%d%d",&L,&U))
    {
        if(L==1)L=2;//特别注意
        int en=U-L;
        F(i,0,en)a[i]=0;
        F(i,1,tot)
        {
            int l=(L-1)/primes[i]+1,r=U/primes[i];
            F(j,l,r)if(j>1)a[j*primes[i]-L]=1;
        }
        int an1,dis1=100*N,an2,dis2=-1,pre=-1;
        F(i,0,en)
        {
            if(!a[i])
            {
                if(pre!=-1)
                {
                    if(i-pre<dis1)dis1=i-pre,an1=pre;
                    if(i-pre>dis2)dis2=i-pre,an2=pre;
                }
                pre=i;
            }
        }
        if(dis2==-1)puts("There are no adjacent primes.");
        else printf("%lld,%lld are closest, %lld,%lld are most distant.
",(ll)an1+L,(ll)an1+L+dis1,(ll)an2+L,(ll)an2+L+dis2);
    }
    return 0;
}