【代码超详解】CometOJ C0275 [USACO]特殊的质数肋骨（质数的判定 + 枚举 + 剪枝）

一、题目描述

1000 ms，256 MB

二、算法分析说明与代码编写指导

N 的范围是 1 到 8，当 N = 8 时，需要验证并列举 10000000 到 99999999 的全部“特殊质数”，直接暴力枚举会超时。
一种方法是：当 N 比较大的时候，只对高位（这里选高 5 位）枚举质数，然后进行验证。当然验证过程不能直接模拟题目所说的“从右边开始切下肋骨”的过程，因为判定较大数的素性是比较费时间的。如果对一个数，截去末尾若干位后判定仍然是质数，但差不多截完了以后又发现剩余的位不是质数，显然是不经济的。所以，应该先截取高位开始更快速地判断，一旦发现不是“特殊质数”，立刻终止判断并返回 false。
另一种方法来自：https://blog.csdn.net/shuati2000/article/details/97816699
从 1 位数开始枚举到 N 位数，枚举过程中发现正在枚举的数已经不是“特殊质数”则跳过该数。这种剪枝方法可以省去非常大的枚举量。最高位自然只有 2，3，5，7。以 N = 8，且正在枚举第二位为例。如果枚举到 21，验证 21 不是质数，那么剩余 6 位都不用再进行枚举，而直接开始枚举高位为 23 的数。

三、AC 代码

法一

#include<cstdio>
#include<algorithm>
#pragma warning(disable:4996)
using namespace std;
unsigned prime[9592] = {
    2,3 }, _PTy, MaxPrime, * prime_end = prime + sizeof(prime) / sizeof(prime[0]);
inline void genprime() {
   
	decltype(_PTy) a = 4, t; bool flag = true;
	for (auto i = prime + 2; i != prime_end;) {
   
		t = sqrt(a); flag = true;
		for (auto j = prime; *j <= t; ++j)if (a % *j == 0) {
    flag = false; break; }
		if (flag) {
    *i = a, ++i; }
		++a;
	}
	MaxPrime =