【区间范围素数筛法】POJ - 2689 Prime Distance

Problem Description

输入a,b代表区间范围[a,b]。问你里面相邻素数距离最近的是那两个素数，相邻素数距离最远的是那两个素数。

思路：

素数筛法复杂度是O(nloglogn)。因为b的范围爆int，直接跑肯定超时。我们可以开两个数组，数组prime[]用来筛[2, sqrt(b))，数组Prime[]用来筛[a, b]。因为b以内的合数的最小质因数一定不超过sqrt(b)。如果有sqrt(b)以内的素数表的话，就可以把埃氏筛法运用到[a,b]上了。也就是说，从[2, sqrt(b))的表中筛得素数的同时，也将其倍数从[a,b]的表中划去，最后剩下的就是区间[a,b)内的素数了。

#include<cstdio>
#include<algorithm>
using namespace std;
#define ll long long
#define nn 1000055
int prime[nn], Prime[nn];
void get_prime(ll a, ll b)
{
//[a,b]区间范围素数筛法。Prime[i - a]=true i是素数
    ll i, j;
    for(i = 0; i*i < b; i++) prime[i] = 1;//初始化
    for(i = 0; i <= b - a; i++) Prime[i] = 1;
    for(i = 2; i*i < b; i++)
    {
        if(prime[i])
        {
            for(j = 2*i; j*j < b; j += i) prime[j] = 0; //划掉区间[2, sqrt(b))的合数
            //核心 j = max(2LL, (a+i-1)/i)*i,让j等于大于等于a的最小的 i的倍数
            for(j = max(2LL, (a+i-1)/i) * i; j <= b; j += i) Prime[j - a] = 0;//划掉区间[a, b]的合数
        }
    }
//求[a,b]区间内相邻素数最近，相邻素数最远
    long long l, r, lt, rt, lc;
    int Min = 0x3f3f3f3f, Max = 0, flag = 0;
    for(i = a; i <= b; i++)
    {
        if(i == 1) continue;
        if(Prime[i - a] && flag)
        {
            if(i - lc < Min)
            {
                Min = i - lc;
                l = lc;
                r = i;
            }
            if(i - lc > Max)
            {
                Max = i - lc;
                lt = lc;
                rt = i;
            }
            lc = i;
            flag++;
        }
        else if(Prime[i - a])
            {lc = i; flag++;}
    }
    if(flag >= 2)//有两个素数，输出
    {
        printf("%lld,%lld are closest, %lld,%lld are most distant.\n", l, r, lt, rt);
    }
    else printf("There are no adjacent primes.\n");
}
int main()
{
    ll a, b;
    while(~scanf("%lld %lld", &a, &b))
    {
        if(a > b) swap(a, b);
        get_prime(a, b);
    }
    return 0;
}