196. 质数距离 - AcWing题库（素数筛，区间，筛两遍）

196. 质数距离 - AcWing题库

分析：

• 要查询的是一段区间的素数情况，筛两遍即可，与Sum of a Function （素数筛，思维转换）类似

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

const int N=1e5+5;

int b[N], pri[N], tot;
void isp(int n=N-5)
{
    b[0]=b[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!b[i]) pri[++tot]=i;
        for(int j=1;pri[j]<=n/i;j++)
        {
            b[pri[j]*i]=1;
            if(i%pri[j]==0) break; 
        }
    }
}
int vis[N<<5],c[N];
signed main()
{
    isp();
    int l,r;
    while(cin>>l>>r)
    {
        for(int i=0;i<=r-l+1;i++) vis[i]=0; 
        for(int i=1;i<=tot && pri[i]<=sqrt(r);i++)
        {
            int st=l/pri[i]*pri[i];
            st=max(st,pri[i]*2);
            if(st<l) st+=pri[i];
            for(int j=st;j<=r;j+=pri[i]) vis[j-l+1]=1;
        }
        if(l==1) vis[1]=1;
        int cnt=0;
        for(int i=l;i<=r;i++)
        {
            if(!vis[i-l+1]) c[++cnt]=i; 
        }
        if(cnt<2) cout<<"There are no adjacent primes."<<endl;
        else 
        {
            int mn=1e9,mx=0,pn,px;
            for(int i=2;i<=cnt;i++)
            {
                int d=c[i]-c[i-1];
                if(d<mn) { mn=d; pn=i; }
                if(d>mx) { mx=d; px=i; }
            }
            printf("%d,%d are closest, %d,%d are most distant.\n",c[pn-1],c[pn],c[px-1],c[px]);
        }
    }
    return 0;
}