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+];
 long long prime[N], cnt;
 bool res[len+];

 void init() { // 筛选50000内的素数。要注意i比较大时 会容易超int的范围
      long long i, j;
      cnt = ;
      memset(isprime, true, sizeof(isprime));
      for (i=; 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, , sizeof(res));
         for (i=; i<cnt; ++i) {
             s = (L-)/prime[i]+;  // 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的时候也不算？？？？？
                     res[j*prime[i]-L] = true;
                 }
             }
         }
         k = -, minn = inf, maxx = -;
         long long dis, m1, m2;
         for (i=; i<=U-L; ++i) { // 求解最优值
             if (!res[i]) {
                 if (k!=-) {
                     dis=i-k;  // 记录两个素数的最远距离。
                     if (dis > maxx) {
                         maxx = dis;
                         m1 = i;
                     }
                     if (dis < minn) {
                         minn = dis;
                         m2 = i;
                     }
                 }
                 if (i+L != ) // 注意1 的时候特殊判断1不是素数。
                     k=i;
             }
         }
         if (maxx == -)
             printf("There are no adjacent primes.\n");
         else printf("%lld,%lld are closest, %lld,%lld are most distant.\n", m2-minn+L, m2+L, m1-maxx+L, m1+L);
     }
     return ;
 }

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