1005 继续(3n+1)猜想（PAT-B）

题目如下：

卡拉兹(Callatz)猜想已经在1001中给出了描述。在这个题目里，情况稍微有些复杂。

当我们验证卡拉兹猜想的时候，为了避免重复计算，可以记录下递推过程中遇到的每一个数。例如对 n=3 进行验证的时候，我们需要计算 3、5、8、4、2、1，则当我们对 n=5、8、4、2 进行验证的时候，就可以直接判定卡拉兹猜想的真伪，而不需要重复计算，因为这 4 个数已经在验证3的时候遇到过了，我们称 5、8、4、2 是被 3“覆盖”的数。我们称一个数列中的某个数 n 为“关键数”，如果 n 不能被数列中的其他数字所覆盖。

现在给定一系列待验证的数字，我们只需要验证其中的几个关键数，就可以不必再重复验证余下的数字。你的任务就是找出这些关键数字，并按从大到小的顺序输出它们。
输入格式：

每个测试输入包含 1 个测试用例，第 1 行给出一个正整数 K (<100)，第 2 行给出 K 个互不相同的待验证的正整数 n (1<n≤100)的值，数字间用空格隔开。
输出格式：

每个测试用例的输出占一行，按从大到小的顺序输出关键数字。数字间用 1 个空格隔开，但一行中最后一个数字后没有空格。
输入样例：

6
3 5 6 7 8 11
输出样例：

7 6

解决方案：

1.思路：此题较为复杂，先要搞懂题目是什么意思，题目要找的“关键数”是在验证别的数字的过程中没有出现过的数字。以样例输入举例：

3  5  6  7  8  11
其中验证3的时候需要判断5、8、4、2、1。
那么5就不是关键数，因为在验证3的时候验证过5了。
验证6的时候需要判断3、5、8、4、2、1.
那么3也不是关键数，因为在验证6的时候验证过3了。
验证7的时候需要判断11、17、26、13、20、10、5、8、4、2、1。
由此可以得出8和11也不是关键数。
在以上判断过程中，6和7都没有在判断其他数字的时候出现过，反而在判断6和7的过程中其他数字出现了。
这说明在这一串数字中只需要判断6和7就可以知道其他数字是否符合卡拉兹猜想，所以6和7是关键数。

可以看出，本题的关键在于找出被覆盖的数和没有被覆盖的数。
2.这里使用散列来寻找这些数。具体为：使用一个数组来存取数据，然后再创立一个数组用来标记数据；将标记数组中的元素值都设置成0。一个数在进行卡拉兹循环的时候，会产生许多中间的数据，将这些产生的中间数据作为数组的下标，在标记数组中将对应位置的元素变为1；例如：

判断3的时候出现了：5、8、4、2、1，那么就将标记数组中下标为5、8、4、2、1的元素值变为1。
这样遍历数组，找值为1或者0的下标，就可以知道哪个值被标记，哪个值没有被标记。
例如上述数组中，jud[5]、jud[8]、jud[4]、jud[2]、jud[1]值为1，说明5、8、4、2、1出现过，已经被覆盖。

3.编码：首先需要写一个卡拉兹循环的函数。
这里特别需要注意的是，判断数组Jud需要很大的内存空间，因为中间数据可能很大。

int jud[100001] = {0};
int loop(int n, int jud[]) 
{
	while(n != 1) {
		if(n % 2 == 0){
			n = n / 2;
			jud[n] = 1;  //将中间数据对应下标的元素加一
		}
		else{
			n = 3 * n + 1;
			n = n / 2;
			jud[n] = 1;  //将中间数据对应下标的元素加一
		}
	} 
}

接着使用一个循环来读取数据，然后将这些数据进行卡拉兹循环，并且标记数组。

int num[101];   //数据数组
for(int i = 0; i < K; i++) {
		cin >> n;
		num[i] = n; 
		loop(n, jud);
	}

接着再用一个循环来找出没有被覆盖过的数。

int final[101];   //关键数数组
for(int i = 0; i < K; i++) {
		if(jud[num[i]] == 0) final[j++] = num[i];   
		//如果当前数据对应的标记数组的元素值为0，说明没有被覆盖
	}

接着排序输出即可。

完整代码：

#include <iostream>
#include <algorithm>
using namespace std;
bool cmp(int a, int b)
{
	return a > b;
}
int loop(int n, int jud[]) 
{
	while(n != 1)
	{
		if(n % 2 == 0){
			n = n / 2;
			jud[n]++;
		}
		else{
			n = 3 * n + 1;
			n = n / 2;
			jud[n]++;
		}
	} 
}
int main()
{
	int jud[100001] = {0};
	int K;
	int n;
	int j = 0;
	cin >> K;
	int num[101];
	int final[101];
	for(int i = 0; i < K; i++) {
		cin >> n;
		num[i] = n; 
		loop(n, jud);
	}
	for(int i = 0; i < K; i++) {
		if(jud[num[i]] == 0) final[j++] = num[i];
	}
	sort(final, final + j, cmp);
	for(int i = 0; i < j; i++) {
		cout << final[i];
		if(i != j - 1)cout << " ";
	}
	return 0;
}