题意:
求区间 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;
}