(Relax 数论1.11)POJ 1595 Prime Cuts(欧拉筛法: 输出1~n区间中中间的2C个素数)

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
/**
 * 2 2 3 17
 2 2 3 29

 2 2 5 37
 */
using namespace std;

const int maxn = 100002;//不要开得太大,否则可能MLE
int su[maxn];
bool u[maxn];
int num = 0;

void prepare() {
	int i, j;
	memset(u, true, sizeof(u));

	for (i = 2; i <= 100001; ++i) {
		if (u[i]) {
			su[++num] = i;
		}

		for (j = 1; j <= num; ++j) {
			if (i * su[j] > 100001) {
				break;
			}

			u[i * su[j]] = false;

			if (i % su[j] == 0) {
				break;
			}
		}
	}
}

int main() {
	prepare();

	int n, c;
	int pri[maxn];
	while (scanf("%d%d", &n, &c) != EOF) {
		int i;
		int j = 1;
		memset(pri, 0, sizeof(pri));

		for (i = 1; i <= n; ++i) {
			if (u[i]) {
				pri[j++] = i;
			}
		}

        printf("%d %d:",n,c);

		int nm = 0;
		if (c * 2 > j) {//如果要求输出的素数的个数>[1,n]内所拥有的素数的个数
			for (i = 1; i < j; ++i) {//全部输出
				printf(" %d", pri[i]);
			}
			printf("\n\n");
		} else {//****否则根据要求输出中间的规定个数的素数
			if (j % 2 == 0) {

				int pos = (j - 2 * c) / 2;//打印中间的2C-1个元素

				for (i = pos + 1;; ++i) {
					nm++;
					if (nm > (2 * c - 1)) {
						break;
					}

					if (pri[i] != 0) {//**这个判断其实可以没有....
						printf(" %d", pri[i]);
					}
				}
				printf("\n\n");
			} else {
				int pos = (j - 2 * c) / 2;
				for (i = pos + 1;; ++i) {
					nm++;
					if (nm > (2 * c)) {
						break;
					}

					if (pri[i] != 0) {
						printf(" %d", pri[i]);
					}
				}

				printf("\n\n");
			}
		}
	}

	return 0;
}