POJ 2689 Prime Distance（筛法求素数）

来源：互联网发布：屏幕触摸检测软件编辑：程序博客网时间：2024/05/18 18:53

http://poj.org/problem?id=2689

题意：

给定一个范围[l,r]求区间[l,r]中距离最近的一对素数和距离最远的一对素数，l,r的值非常打，但是区间比较小，大概在100w左右。

思路：

这题只能采用素数筛法了，但是由于l，r的值实在太大，开不了如此大的空间。但是[l,r]区间长度比较小，我们换个思路，先把2到根号r范围内的素数筛出来，打个表，这样一来，[l,r]区间内的质因子数都在表中了，再筛一次，就能求得[l,r]区间内的素数表，还挺难写的。。。。

代码：

#include <cstdio>#include <cstring>#define LL long longusing namespace std;int prime1[1000000 + 10],prime2[1000000 + 10];bool f[1000000 + 10];void getPrime(int n){    prime1[0] = 0;    memset(f, true, sizeof(f));    for(int i = 2; i <= n; i++)    {        if(f[i])        {            prime1[0]++;            prime1[prime1[0]] = i;            //printf("%d %d\n", i, n);        }        for(int j = 1; ((j <= prime1[0]) && (i * prime1[j] <= n)); j++)        {            f[i * prime1[j]] = false;            if(i % prime1[j] == 0)                break;        }    }}void Prime(int l, int r){    if(l == 1) l = 2;    memset(f, true, sizeof(f));    for(int i = 1; (i <= prime1[0]) && (LL)prime1[i] * prime1[i] <= r; i++)    {        int x = l / prime1[i] + (l % prime1[i] > 0);        if(x == 1) x = 2;        for(int j = x; (LL)j * prime1[i] <= r; j++)            if((LL)j*prime1[i] >= l)                f[j*prime1[i] - l] = false;    }    prime2[0] = 0;    for(int i = 0; i <= r - l; i++)        if(f[i])        {            //printf("%d\n", i + l);            prime2[++prime2[0]] = i + l;        }    if(prime2[0] < 2)    {        printf("There are no adjacent primes.\n");        return;    }    else    {        int x1 = 0, x2 = 0, y1 = 0, y2 = 1000000 + 10;//        for(int i = 1; i <= prime2[0]; i++)//            printf("%d\n", prime2[i]);        for(int i = 2; i <= prime2[0]; i++)        {            if(prime2[i] - prime2[i - 1] > x2 - x1)            {                x2 = prime2[i];                x1 = prime2[i - 1];            }            if(prime2[i] - prime2[i - 1] < y2 - y1)            {                y2 = prime2[i];                y1 = prime2[i - 1];            }        }        printf("%d,%d are closest, %d,%d are most distant.\n", y1, y2, x1, x2);    }}int main(){    int l,r;    getPrime(50000);    while(~scanf("%d %d", &l, &r))    {        Prime(l, r);    }    return 0;}

0 0