质数距离

质数距离

因为U和R的范围在2^31 - 1内，数据太大无论时间和空间都会超，我们需要找到其它的性质
如果一个数d能整除n，则n / d也能整除n，这就让我们发现了一个性质，我们只用筛出2 ~ sqrt(n)中的所有质数，然后再用埃氏筛法筛除[L, R]内的所有约数，即之前筛出的质数的倍数，这样我们只用筛到sqrt(2 ^ 31)内的质数也就是5w之内的质数，大大降低消耗的时间和空间

这里有一个小技巧就是求[L, R]中第一个是p的倍数，即l / p上取整的数用
(l + p - 1) / p * p

1、l是p的倍数
p - 1 / p = 0，等于没有加
2、l不是p的倍数
l mod p >= 1
l + p - 1 = l - 1 + p
p / p = 1，l - 1之后除p不影响下取整的结果

综合上面两点可以得出(l + p - 1) / p就是l / p上取整的结果

---

具体代码如下

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N = 1000010;

int n;
int primes[N], cnt;
bool st[N];

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

int main(){
    
    int l, r;
    while(cin >> l >> r){
        
        init(50000);
        
        memset(st, 0, sizeof st);
        
        for(int i = 0; i < cnt; ++ i){
            
            LL p = primes[i];
            //至少要从2p开始筛，因为p是一个质数所以不用被筛除
            
            for(LL j = max(2 * p, (l + p - 1) / p * p); j <= r; j += p){
                st[j - l] = true;//数据很大但偏移量很小
            }
        }

        cnt = 0;
        for(int i = 0; i <= r - l; ++ i)
            if(!st[i] && i + l >= 2)
                primes[cnt ++] = i + l;


        if(cnt < 2) puts("There are no adjacent primes.");
        else{
            
            int minp = 0, maxp = 0;
            for(int i = 0; i + 1 < cnt; ++ i){
                
                int d = primes[i + 1] - primes[i];
                if(d < primes[minp + 1] - primes[minp]) minp = i;
                if(d > primes[maxp + 1] - primes[maxp]) maxp = i;
            }
            
            printf("%d,%d are closest, %d,%d are most distant.\n",
                    primes[minp], primes[minp + 1], 
                    primes[maxp], primes[maxp + 1]);
        }
    }
    
    return 0;
}