(Relax 数论1.14)POJ 26889 Prime Distance(区间筛素数:求一个区间L\U中距离最近的两个素数和距离最远的两个素数)

这道题是学习素数筛法的经典，应用到了区间筛素数。具体思路是先筛出1到sqrt(2147483647)之间的所有素数，然后再通过已经晒好素数筛出给定区间的素数，关于筛素数的问题，我转载了一篇网上一个人总结的素数总结，里面就有关于筛大数区间的素数的算法。虽然说一切都已经具备了，但是这道题我第一次做的时候RUNTIME ERROR了，原因是声明变量的时候用了long，当数据比较大的时候会超，第二次我改成long long，这才A掉。。总之，还有很多要学的东西，继续努力吧！！

#include
#include 

using namespace std;

/**
 * prime[]:存储小范围的素数表
 * prime1[]:存储某一区间的素数表
 * pcount:小范围素数表中素数的个数
 * pcount1:某一区间的素数表中素数的个数
 * in[]、in1[]用来标记某一个数字是否为素数..
 */
long long prime[46500], prime1[1000010];
long long pcount, p1count;
long long in[46500], in1[1000010];

void getsprime() {//快速筛小范围素数
	for (int i = 2; i < 46500; ++i){
		in[i] = 1;
	}
	for (int i = 2; i < 46500; ++i) {
		if (!in[i]){
			continue;
		}
		for (int j = i * 2; j < 46500; j += i){
			in[j] = 0;
		}
	}
	pcount = 0;
	for (int i = 2; i < 46500; ++i){
		if (in[i]){
			prime[pcount++] = i;
		}
	}
}
void getlprime(long L, long U) {
	if (U < 46500) { //假如区间的最大值就小于46500的话，那么之际用我们刚才求完的就好了
		p1count = 0;
		for (int i = L; i <= U; ++i){
			if (in[i]){
				prime1[p1count++] = i;
			}
		}
	} else { //筛大数区间内的素数的核心部分
		long long k, size = U - L;
		for (long i = 0; i <= size; ++i){
			in1[i] = 1;
		}
		for (long i = 0; i <= pcount && prime[i] * prime[i] <= U; ++i) {
			k = L / prime[i];
			if (k * prime[i] < L){
				++k;
			}
			if (k <= 1){
				++k;
			}
			while (k * prime[i] <= U) {
				in1[k * prime[i] - L] = 0;
				++k;
			}
		}
		p1count = 0;
		for (long i = 0; i <= size; ++i){
			if (in1[i]){
				prime1[p1count++] = i + L;
			}
		}
	}
}

int main() {

	getsprime();
	long long L, U;
	while (scanf("%lld%lld", &L, &U) != EOF) {
		getlprime(L, U);

		long long distance, xmin, ymin, xmax, ymax;
//		long long min=2147483647,max=-1;
		int max = INT_MIN;
		int min = INT_MAX;
		if (p1count < 2) {
			printf("There are no adjacent primes.\n");
		} else {
			long i;
			for (i = 0; i < p1count - 1; ++i) {
				distance = prime1[i + 1] - prime1[i];
				if (distance < min) {
					min = distance;
					xmin = i;
					ymin = i + 1;
				}

				if (distance > max) {
					max = distance;
					xmax = i;
					ymax = i + 1;
				}
			}

			printf("%lld,%lld are closest, %lld,%lld are most distant.\n",
					prime1[xmin], prime1[ymin], prime1[xmax], prime1[ymax]);
		}

	}

	return 0;
}