筛指定区间的素数[区间偏移二次筛法]

引入问题：

---

给定两个整数L和U，你需要在闭区间[L,U]内找到距离最接近的两个相邻质数C1和C2（即C2-C1是最小的），如果存在相同距离的其他相邻质数对，则输出第一对。
同时，你还需要找到距离最远的两个相邻质数D1和D2（即D1-D2是最大的），如果存在相同距离的其他相邻质数对，则输出第一对。
$1\le L<U\le 2^{31}-1,U-L\in [0,1e6]$

---

解题思路：

---

(1)$2^{31}-1\approx 1e9$无论是欧拉筛还是埃氏筛法都是从$[1,n]$开始筛的，这个数据范围都是无法接受的
(2)很明显如果这个一个数的是合数：那么这个数一定有一个质因子$d<=\sqrt{n}$
(3)那么我们就先把$[1,\sqrt{n}]$内的素数筛出来$\sqrt{1e9}\approx 5e4$可以接受
(4)那么我们就把这些素数去筛$[L,U]$的素数
(5)根据埃氏筛法我们要找到$[L,U]$这个区间内第一个大于质数$p$,并且是$p$的倍数之后就每次$+p$一次筛下去；
(6)第一个大于$p$并且是$p$的倍数$\lceil \frac{L}{p}\rceil \ast p\approx \lfloor \frac{L+p-1}{p}\rfloor \ast p$
(7)因为数很大我们可以让$[$区间整体$-L]\in [0,1e6]$

---

代码

---

题目链接(qwq)戳我

#include <bits/stdc++.h>

using namespace std;
const int N = 6e6 + 10;
long long prime[N], cnt;
bool vis[N];

 void init(int n)
 {
     cnt = 0;
     memset(vis,0,sizeof (vis));
     for(int i = 2; i <= n; ++ i)
     {
         if(!vis[i]) prime[cnt ++] = i;
         for(int j = 0; j < cnt && i * prime[j] <= n; ++ j)
         {
            vis[i * prime[j]] = true;
            if(i % prime[j] == 0) break;
         }
     }
 }

int main()
{
    long long l, r;
    while(cin >> l >> r)
    {
        init(50000);//Eula筛
        memset(vis,0,sizeof (vis));
        for(int i = 0; i < cnt; ++ i)
        {
            long long p = prime[i];
            for(long long j = max(2 * p, (l + p - 1) / p * p); j <= r; j += p)
              vis[j - l] = true;//减去偏移量
        }
        cnt = 0;
        for(long long i = 0; i <= r - l; ++ i)
          if(!vis[i] && i + l >= 2)
            prime[cnt ++] = i + l;
        int maxp = 0, minp = 0;
        for(int i = 0; i < cnt - 1; ++ i)
        {
            int d = prime[i + 1] - prime[i];
            if(d > prime[maxp + 1] - prime[maxp]) maxp = i;
            if(d < prime[minp + 1] - prime[minp]) minp = i;
        }
        if(cnt < 2) puts("There are no adjacent primes.");
        else 
        {
             cout << prime[minp] << "," << prime[minp + 1] << " are closest, ",
             cout << prime[maxp] << "," << prime[maxp + 1] << " are most distant.\n";
        }
    }
    return 0;
}