半数集问题

问题描述：

给定一个自然数n给定一个自然数n，由n 开始可以依次产生半数集set(n)中的数如下。
(1) n∈set(n)；
(2) 在n 的左边加上一个自然数，但该自然数不能超过最近添加的数的一半；
(3) 按此规则进行处理，直到不能再添加自然数为止。

初步分析：

        以set(8)为例，分析半数集的计算过程。首先set(8)中含有8本身，即set(8) = {8}，在此基础上可以在8的左边添加1,2,3,4这四个数，即set(8) = {8, 18, 28, 38, 48}，对于18，没有可以添加的数了；对于28，可以在左边添加1，即128；对于38，可以在左边添加1，即138；对于48，可以在左边添加1,2,即148,248，当前set(8) = {8, 18, 28, 38, 48, 128, 138, 148, 248}。现在仅有248的左边可以添加一个数了，即1248，最后的半数集为set(8) = {8, 18, 28, 38, 48, 128, 138, 148, 248, 1248}。该半数集的生成过程如下图所示：

递归关系：

        从上面set(8)的生成过程不难发现，set(8)其实由其本身,set(1),set(2),set(3),set(4)组成，即set(8) = 1 + set(1) + set(2) + set(3) + set(4)。这其实是把set(8)分治成了n/2个规模更小的问题，我们可以用递归函数实现这一算法，递归关系如下：

set(n) = 1 + set(1) + set(2) + set(3) + … + set(n/2)

        根据上述递归关系，我们也可以发现下面的关系：

set(n) = set(n + 1 - n%2)

实现代码：

#include <bits/stdc++.h>
using namespace std;

//使用cache数组存储已计算的set(i) 
int cache[100];
int num_of_set(int n){
	if(n == 1) return 1;
	int tmp = 1;
	for(int i = 1; i < n/2 + 1; i++){
		if(cache[i] == 0)
			cache[i] = cache[i + 1 - i%2] = num_of_set(i);
		tmp += cache[i];
	}		
	return tmp;
}

int main(){
	int num;
	memset(cache, 0, 400);
	cache[1] = 1;
	
	//打开文件流 
	ifstream input("input.txt", ios::in);
	if(!input){
		cout << "file open error!\n";
		exit(1);
	}
	
	ofstream output("output.txt", ios::out | ios::trunc);
	if(!output){
		cout << "file open error!\n";
		exit(1);
	}
	
	input >> num;
	output << num_of_set(num);
	
	//关闭文件流 
	input.close();
	output.close();
	return 0;
}

        由于直接递归会造成很多的重复计算，造成算法时间和空间上的性能降低，所以在此处使用cache数组储存已经计算过的set(i)的值以提高算法性能。

正确性检验：