POJ2689 Prime Distance

用线性筛打出质数，在打标记

引用：

首先我们发现：R-LR−L 的范围很小，我们应该要能够快速求出 L\sim RL∼R 之间的质数。

显然有推论：任意一个合数 xx 必定包含一个不超过 \sqrt xx​ 的质因子。

所以我们可以筛出 [1,\sqrt R][1,R​] 之间的所有质数，对于每个质数 pp，把 [L,R][L,R] 中能被 pp 整除的数标记为合数。最终没有被标记的数就是质数，对相邻的质数两两比较，找出差值最小和最大的即可。

时间复杂度：O(\sum_{p\le \sqrt R}\frac{R-L}{p})=O(\sqrt R\log^2\sqrt R+(R-L)\log^2 R)O(∑p≤R​​pR−L​)=O(R​log2R​+(R−L)log2R)（该复杂度摘自：李煜东《算法进阶指南》）

Time	188ms
Memory	11956kB

 

代码

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = 1 << 16;
int prime[maxn], mm, p[1000005];
long long a[1000005];
bool isprime[maxn];
void primes()
{
    memset(isprime,true,sizeof(isprime));
    isprime[1] = false;
    mm = 0;
    for(int i = 2; i <= maxn; i++){
        if(isprime[i]){
            prime[++mm] = i;
        }
        for(int j = 1; j <= mm && prime[j]*i <= maxn; j++){
            isprime[prime[j]*i] = false;
            if(i % prime[j] == 0) break;
        }
    }
}
int main()
{
    long long n, m, cnt;
    primes();
    while(scanf("%lld%lld",&n,&m) != EOF){
        cnt = 0;
        memset(a,0,sizeof(a));
        memset(p,0,sizeof(p));
        if(n == 1) n = 2;
        for(int i = 1; i <= mm; i++){
            if(prime[i] > m) break;
            for(long long j = ceil(n*1.0/prime[i]); j <= floor(m*1.0/prime[i]); j++){
                while(j <= 1) j++;
                if(prime[i]*j < n || prime[i]*j > m) continue;
                p[prime[i]*j-n] = 1;
            }
        }
        for(long long i = 0; i <= m-n; i++){
            if(!p[i]) a[++cnt] = i+n;
        }
        if(cnt == 0 || cnt == 1){
            printf("There are no adjacent primes.\n");
        }
        else{
            long long l1, r1, l2, r2, minn = 0x3f3f3f3f, maxx = -1;
            for(long long i = 1; i < cnt; i++){
                if(a[i+1]-a[i] < minn){
                    minn = a[i+1] - a[i];
                    l1 = a[i];
                    r1 = a[i+1];
                }
                if(a[i+1]-a[i] > maxx){
                    maxx = a[i+1] - a[i];
                    l2 = a[i];
                    r2 = a[i+1];
                }
            }
            printf("%lld,%lld are closest, %lld,%lld are most distant.\n",l1,r1,l2,r2);
        }
    }
    return 0;
}