【POJ 2689】 Prime Distance

【题目链接】

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

【算法】

我们知道，一个在区间[l,r]中的合数的最小质因子必然不超过sqrt(r)

那么，先暴力筛出1-50000中的质数，对于每个询问，用筛出的质数标记[l,r]中的合数，即可

【代码】

#include <algorithm>
#include <bitset>
#include <cctype>
#include <cerrno>
#include <clocale>
#include <cmath>
#include <complex>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <limits>
#include <list>
#include <map>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <utility>
#include <vector>
#include <cwchar>
#include <cwctype>
#include <stack>
#include <limits.h>
using namespace std;
#define MAXP 50000
#define MAXD 1000010
const int INF = 2e9;

int i,j,last,mx,mn,l,r,L,U;
vector<int> P;
bool not_prime[MAXD];
pair<int,int> C,D;
bool flag;

inline void init()
{
        int i,j,tmp;
        static int f[MAXP+];
        for (i = ; i <= MAXP; i++)
        {
                if (!f[i])
                {
                        P.push_back(i);
                        f[i] = i;
                }
                for (j = ; j < P.size(); j++)
                {
                        tmp = i * P[j];
                        if (tmp > MAXP) break;
                        f[tmp] = P[j];
                        if (f[i] == P[j]) break;
                }
        }
}

int main() {

        init();
        while (scanf("%d%d",&L,&U) != EOF)
        {
                memset(not_prime,false,sizeof(not_prime));
                flag = false;
                last = -;
                mx = ; mn = INF;
                for (i = ; i < P.size(); i++)
                {
                        if (L % P[i] == ) l = L / P[i];
                        else l = L / P[i] + ;
                        r = U / P[i];
                        for (j = max(l,); j <= r; j++) not_prime[P[i]*j-L] = true;
                }
                if (L == ) not_prime[] = true;
                for (i = ; i <= U - L; i++)
                {
                        if (!not_prime[i])
                        {
                                if (last == -)
                                {
                                        last = i;
                                        continue;
                                }
                                if (i - last < mn)
                                {
                                        flag = true;
                                        mn = i - last;
                                        C = make_pair(last+L,i+L);
                                }
                                if (i - last > mx)
                                {
                                        flag = true;
                                        mx = i - last;
                                        D = make_pair(last+L,i+L);
                                }
                                last = i;
                        }
                }
                if (flag) printf("%d,%d are closest, %d,%d are most distant.\n",C.first,C.second,D.first,D.second);
                else printf("There are no adjacent primes.\n");
        }

        return ;

}