POJ 2689 筛法求素数

DES：给出一个区间[L, U]。找出这个区间内相邻的距离最近的两个素数和距离最远的两个素数。其中1<=L<U<=2147483647 区间长度不超过1000000.

思路：因为给出的U的范围超出了int的最大。所以不能直接打出1-U的素数表。【我们知道。利用素数筛法时是把所有的合数都筛掉了。那么我们也想办法把L-U内的合数都筛掉就可以了。那么只要用sqrt(U)内的素数去筛掉L-U的合数就可以了。】只是大概理解。。。。。。。

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

#define N 500000
#define len 10000000
#define inf 0x7fffffff
bool isprime[N+5];
long long prime[N], cnt;
bool res[len+5];

void init() { // 筛选50000内的素数。要注意i比较大时 会容易超int的范围
     long long i, j;
     cnt = 0;
     memset(isprime, true, sizeof(isprime));
     for (i=2; i<=N; ++i) {
        if (isprime[i]) {
            prime[cnt++] = i;
            if (i*i <= N) {
                for (j=i*i; j<=N; j+=i) {
                    isprime[j] = false;
                }
            }
        }
     }
}

int main() {
    long long L, U, i, j, k, minn, maxx, s, t;
    init();
    while(~scanf("%lld%lld", &L, &U)) {
        memset(res, 0, sizeof(res));
        for (i=0; i<cnt; ++i) {
            s = (L-1)/prime[i]+1;  // s和t表示了当前范围里的数对这个素数的最小倍数和最大倍数。 比如说 区间[3, 9] 对素数2 的 s 和 t 应是 2， 4
            t = U/prime[i];   // 由这个例子可以知道为什么 L-1 最后再+1。t就可以直接除。
            for (j=s; j<=t; ++j) { // 筛选出区间内的素数
                if (j>1) { // 不知道为什么==1的时候也不算？？？？？
                    res[j*prime[i]-L] = true;
                }
            }
        }
        k = -1, minn = inf, maxx = -1;
        long long dis, m1, m2;
        for (i=0; i<=U-L; ++i) { // 求解最优值
            if (!res[i]) {
                if (k!=-1) {
                    dis=i-k;  // 记录两个素数的最远距离。
                    if (dis > maxx) {
                        maxx = dis;
                        m1 = i;
                    }
                    if (dis < minn) {
                        minn = dis;
                        m2 = i;
                    }
                }
                if (i+L != 1) // 注意1 的时候特殊判断1不是素数。
                    k=i;
            }
        }
        if (maxx == -1)
            printf("There are no adjacent primes.
");
        else printf("%lld,%lld are closest, %lld,%lld are most distant.
", m2-minn+L, m2+L, m1-maxx+L, m1+L);
    }
    return 0;
}

代码可以理解。裸敲好像。。。还是有一点难度。明天再敲一遍。