区间筛

题目链接

题意：

求区间 L, R 内的相邻的最接近的素数和相邻的最远的素数，区间的端点很大，但是区间长度不大，输出最接近的一对素数和最远的一对相邻素数

区间筛即可，第一次接触区间筛，所以wa了好多次

区间筛是埃式筛的一种变形，回顾埃式筛原理，首先是从小到大的，对于一个数字 n ，如果它是素数就把对应的倍数筛掉，遇到的第一个没被筛掉的数就是素数，其实相当于遍历的1 到 n - 1，发现没有哪个数字能把 n 整除，说明 n 是素数。

区间筛也是同样的原理，首先埃式筛得到 2 到 b ^ 0.5 的全部素数，把素数在区间 [L，R]内的一倍以上的倍数筛掉，最后没有被筛掉的数就是素数

#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <vector>
using namespace std;
bool is[1000010], c[1000010];
typedef long long ll;
ll a, b;
ll cnt(ll n){
	for (ll i = a; i <= b; i++){
		if(i % n == 0 && i != n)return i;
	}
	return -1;
}
vector<ll> ans;
int main()
{
	while (~scanf("%lld %lld", &a, &b)){
		ans.clear();
		memset(is, true, sizeof is);
		memset(c, true, sizeof c);
		ll sq = sqrt(b);
		if(a <= 1)for (ll i = a; i <= 1; i++)c[i - a] = false;
		is[1] = false;
		is[0] = false;
		for (ll i = 2; i <= sq; i++){
			if(is[i]){
				for (ll j = i + i; j <= sq; j += i){
					is[j] = false;
				}
				ll t;
				if(a == i){
					t = i + i;
				}
				else{
					if(a > i){
						if(a % i == 0)t = a;
						else t = (a / i + 1) * i;
					}
					else{
						t = i + i;
					}
				}
				for (ll j = t; j <= b; j += i){
					c[j - a] = false;
				}
			}
		}
		for (ll j = a; j <= b; j++){
			if(c[j - a])ans.push_back(j);
		}
		if(ans.size() < 2)puts("There are no adjacent primes.");
		else{
			ll a1, a2, b1, b2;
			a1 = -99999999;
			a2 = 999999999;
			b1 = b2 = 0;
			int len = ans.size();
			for (int i = 1; i < len; i++){
				if(a2 - a1 > ans[i] - ans[i - 1]){
					a2 = ans[i];
					a1 = ans[i - 1];
				}
				if(b2 - b1 < ans[i] - ans[i - 1]){
					b2 = ans[i];
					b1 = ans[i - 1];
				}
			}
			printf("%lld,%lld are closest, %lld,%lld are most distant.\n", a1, a2, b1, b2);
		}
	}
	return 0;
}