poj 2689 Prime Distance 二次筛素数

这个是比较经典的二次筛素数的题目。

给出L，U，求在L~U的区间中 相邻两个素数距离最大的和相邻两个素数距离最小的

题目中的L,U的范围是1~2,147,483,647。但是题目中也说了，U-L < 1000000;

解题的思路是，先筛出1~sqrt(2,147,483,647)内的素数，然后用这些素数去筛L~U中的合数。

思路还是很听简单的。

有两点要注意：这个题int可能会越界的，其中有一个 < U的判断，那个地方会越界，血的教训那，我在这个地方错了好久！还有一点是记得处理1哈

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

#define M 47341
#define N 2000010

int prime[N];
bool isPrime[N];
int size;
int L,D;

void getPrime()
{
    size = 0;
    memset(isPrime,1,sizeof(isPrime));
    for(int i=2; i < M; i++)
    {
        if(isPrime[i]) prime[size++] = i;
        for(int j=0; j<size && prime[j]*i < M; j++)
        {
            isPrime[ prime[j]*i ] = 0;
            if(i%prime[j] == 0) break;
        }
    }
}

int prime2[N];
bool isPrime2[N];
int size2;
void getPrime2()
{
    size2 = 0;
    memset(isPrime2,1,sizeof(isPrime2));
    for(int i = 0; i < size && prime[i]*prime[i] <= D; i++)
    {
        int j = L/prime[i]+(L%prime[i] != 0 );
        if(j == 1) j++;
        for(; prime[i] <= D/j; j++) // 这儿是会溢出的！！！教训啊 最后改成了除，要么把这儿改成long long
        {
            isPrime2[ j*prime[i]-L ] = 0;
        }
    }
    for(long long i = L; i <= D; i++)
        if( isPrime2[i-L] ) prime2[size2++] = i;
}

int main()
{
    getPrime();
    while(cin >> L >> D)
    {
        if(L == 1 && D == 2)  {cout << "There are no adjacent primes.\n";continue;}
        else if(L == 1) L++;
        getPrime2();
        int mi = M,ma = -1;
        int t,tt;
        for(int i = 0; i < size2-1; i++)
        {
            if(mi > prime2[i+1]-prime2[i])
            {
                mi = prime2[i+1]-prime2[i];
                t = i;
            }
            if(ma < prime2[i+1]-prime2[i])
            {
                ma = prime2[i+1]-prime2[i];
                tt = i;
            }
        }
        if(size2 <= 1)
            cout << "There are no adjacent primes.\n";
        else cout << prime2[t] << "," << prime2[t+1] << " are closest, " << prime2[tt] << "," << prime2[tt+1] << " are most distant.\n";
    }
    return 0;
}