Acwing196. 质数距离

题目

给定两个整数L和U，你需要在闭区间[L,U]内找到距离最接近的两个相邻质数C1和C2（即C2-C1是最小的），如果存在相同距离的其他相邻质数对，则输出第一对。

同时，你还需要找到距离最远的两个相邻质数D1和D2（即D1-D2是最大的），如果存在相同距离的其他相邻质数对，则输出第一对。

输入格式

每行输入两个整数L和U，其中L和U的差值不会超过1000000。

输出格式

对于每个L和U ，输出一个结果，结果占一行。

结果包括距离最近的相邻质数对和距离最远的相邻质数对。（具体格式参照样例）

如果L和U之间不存在质数对，则输出“There are no adjacent primes.”。

数据范围

1 < L < U < 2^31-1

输入样例：

2 17
14 17

输出样例：

2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

题解

L和U的差值不会超过1000000,那就对于 x,用 [2,sqrt(x)] 间的素数把[L, U]中的素数筛出来,然后for一遍看哪两个之间差最大
技巧:
for(ll j = max(2ll*p,(l+p-1)/p * p) ; j<=r ; j+=p ){
w[j - l] = 1;
}
1. max(2p,(l + p -1) / p * p) 表示大于等于2p的第一个p的倍数
2. 因为区间比较大,可以用 j - l 把j变小,最后算的时候再加上就行

int v[N],prime[N]; 
int w[N],q=0;

int init(int n){   //线性筛选1 ~ N之间的质数
	 memset(v,0,sizeof(v));
	 int m = 0;
	 for(int i=2;i<=n;i++){
	 	if(!v[i]){
	 		v[i] = i;prime[++m] = i;
		 }
		 for(int j=1;j<=m;j++){
		 	if(prime[j] > v[i] || prime[j] > n/i)break;
		 	v[i*prime[j]] = prime[j];
		 }
	 }
	 return m;
}

int main(){
	while(cin >> l >> r){
		int m = init(60000);
		memset(w,0,sizeof(w));
	 	for(int i=1;i<=m;i++){  //l...r 筛一遍 
	 		int p = prime[i];
	 		for(ll j = max(2ll*p,(l+p-1)/p*p);j<=r;j+=p){
	 			w[j - l]  = 1;
			 }
	 }	
		q = 0;
		for(int i=0;i<=r-l;i++){
			if(!w[i] && i+l > 1)prime[++q] = i+l;
		}
		if(q <= 1)printf("There are no adjacent primes.\n"); 
		else {
			int minx = 1000000,maxn = 0,x1,y1,x2,y2;
			for(int i=2;i<=q;i++){
				if(prime[i] - prime[i-1] < minx){
					minx = prime[i] - prime[i-1];
					x1 = prime[i-1];y1 = prime[i];
				}
				if(prime[i] - prime[i-1] > maxn){
					maxn = prime[i] - prime[i-1];
					x2 = prime[i-1];y2 = prime[i];
				}
			}
			printf("%d,%d are closest, %d,%d are most distant.\n",x1,y1,x2,y2);
			
		}
	}
	return 0;
}