poj 2689 Prime Distance（大数区间素数筛法）

题目链接：http://poj.org/problem?id=2689

题意：给定区间[L,R]，求区间内距离最近的相邻素数对和距离最远的相邻素数对，区间长度不超过1e6。

解题方案：用筛法求出[L,R]的所有素数——利用“合数n一定有小于或等于sqrt(n)的素数因子“这条性质，先预处理出sqrt(2,147,483,647)范围内的所有素数，然后用它们筛掉所有在区间[L,R]内的合数。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>

using namespace std;

#define FOR(i,k,n) for(int i=k;i<n;i++)
#define FORR(i,k,n) for(int i=k;i<=n;i++)
#define scan(a) scanf("%d",&a)
#define scann(a,b) scanf("%d%d",&a,&b)
#define scannn(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define mst(a,n)  memset(a,n,sizeof(a))
#define ll long long
#define N 100005
#define mod 1000000007
#define INF 0x3f3f3f3f

const double eps=1e-8;
const double pi=acos(-1.0);

ll L,R;
int vis[N*10];
int prime[N],prime1[N*10];
int cnt,cnt1;
int Min,Max,minl,minr,maxl,maxr;

void Init()
{
    mst(vis,0);
    cnt=0;
    FOR(i,2,N)
    {
        if(!vis[i])
        {
            prime[cnt++]=i;
            for(int j=i+i;j<N;j+=i)
                vis[j]=1;
        }
    }
}

void range_prime()//大数区间筛素数
{
    mst(vis,0);
    FOR(i,0,cnt)
    {
        ll b=L/prime[i];
        if(b*prime[i]<L)  //将b*prime[i]定位到第一个大于或等于L的合数的位置
            b++;          //注意b==0时和b*prime[i]==L并且b==1时的情况
        if(b<=1)          //这两条if语句也可以合并为 while(b*prime[i]<L||b<=1) b++;
            b++;
        for(ll j=b*prime[i];j<=R;j+=prime[i]) //b>=2
            vis[j-L]=1; //节省空间，否则空间复杂度太高
    }
    if(L==1) vis[0]=1; //如果左边界到了1，需要特判，1不是素数，所以要被筛掉

    cnt1=0;
    for(ll i=L;i<=R;i++)
    {
        if(!vis[i-L])
            prime1[cnt1++]=i;
    }
}

void solve()
{
    range_prime();
    Min=INF; Max=0;
    minl=minr=maxl=maxr=-1;
    FOR(i,1,cnt1)
    {
        int d=prime1[i]-prime1[i-1];
        if(d<Min) Min=d,minl=prime1[i-1],minr=prime1[i];
        if(d>Max) Max=d,maxl=prime1[i-1],maxr=prime1[i];
    }
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);

    Init();
    while(cin>>L>>R)
    {
        solve();
        if(cnt1<2) printf("There are no adjacent primes.\n");
        else printf("%d,%d are closest, %d,%d are most distant.\n",minl,minr,maxl,maxr);
    }
    return 0;
}