特立独行的幸福

对一个十进制数的各位数字做一次平方和，称作一次迭代。如果一个十进制数能通过若干次迭代得到 1，就称该数为幸福数。1 是一个幸福数。此外，例如 19 经过 1 次迭代得到 82，2 次迭代后得到 68，3 次迭代后得到 100，最后得到 1。则 19 就是幸福数。显然，在一个幸福数迭代到 1 的过程中经过的数字都是幸福数，它们的幸福是依附于初始数字的。例如 82、68、100 的幸福是依附于 19 的。而一个特立独行的幸福数，是在一个有限的区间内不依附于任何其它数字的；其独立性就是依附于它的的幸福数的个数。如果这个数还是个素数，则其独立性加倍。例如 19 在区间[1, 100] 内就是一个特立独行的幸福数，其独立性为 2×4=8。

另一方面，如果一个大于1的数字经过数次迭代后进入了死循环，那这个数就不幸福。例如 29 迭代得到 85、89、145、42、20、4、16、37、58、89、…… 可见 89 到 58 形成了死循环，所以 29 就不幸福。

本题就要求你编写程序，列出给定区间内的所有特立独行的幸福数和它的独立性。

输入格式：

输入在第一行给出闭区间的两个端点：1<A<B≤10⁴。

输出格式：

按递增顺序列出给定闭区间 [A,B] 内的所有特立独行的幸福数和它的独立性。每对数字占一行，数字间以 1 个空格分隔。

如果区间内没有幸福数，则在一行中输出 SAD。

题目隐含的条件：各位数字平方和的取值范围：

x是一位数，当x=9时平方和取到最大值9²=81

x是两位数，当x=99时平方和取到最大值9²+9²=162

x是三位数，当x=999时平方和取到最大值9²+9²+9²=243

x是四位数，当x=9999时平方和取到最大值9²+9²+9²+9²=324

因此当1 < x ≤ 10⁴ 时，无论迭代多少次，各位数字平方和总是不会超过324

#include <iostream>
using namespace std;

int sqr_sum(int a) {
	int sum = 0;
	while (a) {
		int mod = a % 10;
		sum += mod * mod;
		a /= 10;
	}
	return sum;
}

int quicksqrt(unsigned int M) {
	int N, i;
	int tmp, ttp;    // 结果、循环计数

	if (M == 0)
		return 0;

	N = 0;
	tmp = (M >> 30);           // 获取最高位：B[m-1]
	M <<= 2;

	if (tmp > 1)
		tmp -= ++N;
	for (i = 15; i > 0; i--) { // 求剩余的15位
		N <<= 1;               // 左移一位
		tmp <<= 2;
		tmp += (M >> 30);      // 假设
		ttp = N;
		ttp = (ttp << 1) + 1;
		M <<= 2;
		if (tmp >= ttp) {      // 假设成立
			tmp -= ttp;
			N ++;
		}
	}
	return N;
}

int prime(int a) {
	int val[10] = {0, 0, 1, 1, 0, 1, 0, 1, 0, 0};
	if (a < 10)
		return val[a];
	if (a % 10 == 0 || a % 10 == 5)
		return 0;
	if (a & 1) {
		int sum = 0;
		int temp = a;
		while (temp) {
			sum += temp % 10;
			temp /= 10;
		}
		if (sum % 3 == 0)
			return 0;
		int sqrt = quicksqrt(a);
		for (int i = 7; i <= sqrt; i += 2)
			if (a % i == 0)
				return 0;
		return 1;
	}
	return 0;
}

int main () {
	int a, b;
	cin >> a >> b;

	int *num = new int[325]();
	int *dlx = new int[b - a + 1](); // 独立性
	int *tldx = new int[b - a + 1](); // -1:依附 0:未知 1:特立独行
	int temp;
	for (int i = a; i <= b; i++) {
		for (int j = 0; j < 325; j++)
			if (num[j] == -1)
				num[j] = 0;

		temp = i;
		for (int cnt = 1;; ++cnt) {
			temp = sqr_sum(temp);

			if (temp == 1) {	// 若i幸福 就假设i是特立独行的
				num[i] = dlx[i - a] = cnt;
				if (tldx[i - a] == 0)
					tldx[i - a] = 1;
				break;
			}
			if (a <= temp && temp <= b) {
				tldx[temp - a] = -1;
			}
			if (num[temp] == 0) {
				num[temp] = -1;
				continue;
			} else if (num[temp] == -1) { // 死循环 不幸福
				break;
			} else { // 有幸福数依附于i
				if (tldx[i - a] == 0) {
					dlx[i - a] = cnt + num[temp];
					tldx[i - a] = 1;
					break;
				}
			}

		}
	}

	int flag = 1; // 1:没有幸福数
	for (int i = 0; i <= b - a; i++) {
		if (tldx[i] == 1) {
			flag = 0;
			if (prime(i + a))
				cout << i + a << " " << dlx[i] * 2 << endl;
			else
				cout << i + a << " " << dlx[i] << endl;
		}
	}
	if (flag)
		cout << "SAD" << endl;

	return 0;
}

参考文献：

[1]开平方的快速算法(C程序)_liht_1634的博客-CSDN博客_c语言快速开方