acwing-质数距离（筛选区间质数）

质数距离

由于线性筛只能从头开始筛起，所以对于比较数值比较大的一个区间，采用线性筛无法在较好的时间内筛出结果。

可以筛出 1 - sqrt n 中的所有质数，对于每个质数，将他在区间内的倍数标记，所有未被标记的则是区间内所有质数。

关键：
区间左端点上取整，右端点下取整。
区间内的数过大，数组存不了可以减去区间左端点。

代码：

#include<bits/stdc++.h>
using namespace std;

const int N = 45000;

int prime[N],idx,l,r;
int arr[1000010];
int a,b,c,d;
bool st[N];

void get_prime(int n)
{
    for(int i = 2; i <= n; i ++ )
    {
        if(!st[i])
            prime[idx ++ ] = i;

        for(int j = 0; prime[j] * i <= n; j ++ )
        {
            st[prime[j] * i] = true;

            if(i % prime[j] == 0)
                break;
        }
    }
}

int main()
{
    get_prime(N);

    while(cin>>l>>r)
    {
        if(l == 1)//防止出现1，2数据
            l = 2;

        memset(arr,false,sizeof(arr));

        for(int i = 0; i < idx; i ++ )//筛出区间内的合数
        {
            int p = prime[i];

            for(int j = (l - 1) / p + 1; j <= r / p; j ++ )
                if(j > 1)
                    arr[j * p - l] = true;
        }

        int p = -1,minn = 0x3f3f3f3f,maxx = 0;

        for(int i = 0; i <= r - l; i ++ )
        {
            if(!arr[i])
            {
                if(p == -1)//找到第一个数
                {
                    p = i;
                    continue;
                }

                if(i - p < minn)
                {
                    minn = i - p;
                    a = p,b = i;
                }

                if(i - p > maxx)
                {
                    maxx = i - p;
                    c = p,d = i;
                }

                p = i;
            }
        }

        if(maxx == 0)
            cout<<"There are no adjacent primes."<<endl;
        else
            cout<<a + l<<","<<b + l<<" are closest, "<<c + l<<","<<d + l<<" are most distant."<<endl;
    }
    return 0;
}