YBTOJ B. 【例题2】质数距离

 题面很简单，但当时怎么都做不出（又看题解，又看博客的，现在看好水）

拿到数据后，其实不难发现我们不必要去处理所有的质数，我们只需要将R开根号后的数内的质数求出来然后再进行一遍埃式筛，就可以得到答案了！

Code

#include<bits/stdc++.h>
#define int long long

using namespace std;
int p[10000100];
int cnt;
bool vis[1000010];
void init(int n){
    memset(vis,false,sizeof vis);
    cnt=0;
    for(int i=2;i<=n;i++){
        if(!vis[i]) p[cnt++]=i;
        for(int j=0;p[j]*i<=n;j++){
            vis[p[j]*i]=true;
            if(i%p[j]==0) break;
        }
    }
}
int l,r,sum;
signed main(){
	while(cin>>l>>r){
		//scanf("%lld%lld",&l,&r);
		init(50000);
		memset(vis,0,sizeof(vis));
		for(int i=0;i<cnt;i++){
            int pri=p[i];
            for(int j=max(2*pri,(l+pri-1)/pri*pri);j<=r;j+=pri){
                vis[j-l]=true;
            }
        }
		cnt=0;
		for(int i=0;i<=r-l;i++){
			if(!vis[i]&&i+l>1) p[cnt++]=i+l;
		}
		if(cnt<2) puts("There are no adjacent primes.");
		else{
			int minn=0,maxn=0;
			for(int i=0;i+1<cnt;i++){
				if(p[i+1]-p[i]>p[maxn+1]-p[maxn]) maxn=i;
				if(p[i+1]-p[i]<p[minn+1]-p[minn]) minn=i;
			}
			printf("%d,%d are closest, %d,%d are most distant.\n",p[minn],p[minn+1],p[maxn],p[maxn+1]);
		} 
	}
	return 0;
}