poj2689(素數區間篩法模板)

題目鏈接: http://poj.org/problem?id=2689

 

題意: 給出一個區間 [l, r] 求其中相鄰的距離最近和最遠的素數對 . 其中 1 <= l <  r <= 2,147,483,647, r - l <= 1e6 .

 

思路: 素數區間篩

要找到 [l, r] 中相鄰最近和最遠的素數對肯定是需要找出 [l, r] 內所有素數 . 但是無論是直接線性打表還是暴力都處理不了這么大的數據 .

可以先給 sqrt(r) 內的素數打個表, 再用 sqrt(r) 內的素數去篩選 [l, r] 內的合數, 然后再遍歷一次 [l, r], 記錄答案即可 .

 

代碼:

#include <iostream>
#include <stdio.h>
#include <string.h>
#define ll long long
using namespace std;

const int MAXN = 1e6 + 10;
const int MAX = 1e5;
int prime[MAX], tag[MAX], vis[MAXN], tot;

void get_prime(void){
    for(int i = 2; i < MAX; i++){
        if(!tag[i]){
            prime[tot++] = i;
            for(int j = 2; j * i < MAX; j++){
                tag[j * i] = 1;
            }
        }
    }
}

ll Max(ll a, ll b){
    return a > b ? a : b;
}

int main(void){
    get_prime();
    ll l, r;
    while(~scanf("%lld%lld", &l, &r)){
        memset(vis, 0, sizeof(vis));
        for(int i = 0; i < tot; i++){
            ll a = (l + prime[i] - 1) / prime[i];
            ll b = r / prime[i];
            for(int j = Max(2, a); j <= b; j++){ // 篩[l, r]內的合數
                vis[prime[i] * j - l] = 1; //減個l方便標記,輸出答案時加回去即可
            }
        }
        if(l == 1) vis[0] = 1; // 注意這個1並不是素數
        ll cnt = -1, sol1 = MAXN, sol2 = 0, x1, y1, x2, y2;
        for(int i = 0; i <= r - l; i++){
            if(vis[i] == 0){
                if(cnt != -1){
                    if(sol1 > i - cnt){
                        x1 = cnt;
                        y1 = i;
                        sol1 = i - cnt;
                    }
                    if(sol2 < i - cnt){
                        x2 = cnt;
                        y2 = i;
                        sol2 = i - cnt;
                    }
                }
                cnt = i;
            }
        }
        if(sol2 == 0) puts("There are no adjacent primes.");
        else printf("%lld,%lld are closest, %lld,%lld are most distant.\n", x1 + l, y1 + l, x2 + l, y2 + l);
    }
    return 0;
}