程序设计与分析——半数集问题

一、题目描述及输入输出要求

•问题描述：给定一个自然数n，由n开始可以依次产生半数集set(n)中的数如下：

   （1）n∈set(n)；

   （2）在n的左边加上一个自然数，但该自然数不能超过最近添加的数的一半；

   （3）按此规则进行处理，直到不能再添加自然数为止。例如，set(6)={6, 16, 26, 126, 36, 136}。半数集set(6)中有6个元素。注意，该半数集是多重集。

•算法设计：对于给定的自然数n，计算半数集set(n)中的元素个数。

•数据输入：输入数据由文件名为input.txt的文本文件提供。每个文件只有一行，给出整数n（0<n<1000）。

•结果输出：将计算结果输出到文件output.txt。输出文件只有一行，给出半数集set(n)中的元素个数。

                 输入文件示例                      输出文件示例

                     input.txt                                output.txt

                       6                                               6

二、题目分析与设计思路

可以先分析set(6)这个案例集合是如何得到的：先是例举数字6左边一位的值可以是1，2，3;然后再分别例举数字1，2，3可以是哪些值（只要不超过最近数的一半），不断重复上述操作直到最左边的值是1为止（因为没有比1小的数，列到1之后不能再往左列了），所有列举出来的都是set(6)的元素，用一个变量来累加得到元素的总个数。发现每次列举相当于set（6）分解出来的子问题set(1),set(2),set(3)，而每个子问题又可以分解出子问题，很容易想到要用递归树来解决这个问题。

三、代码与测试结果

#include <iostream>
#include<fstream>
#include<vector>
#include <cmath>
using namespace std;
vector<int> memo(100, 0);

int set(int n) {
	    int c = 1;
		if (memo[n]>0) {
			return memo[n];
		}
		else {
			for (int i = 1; i <= n / 2; i++) {
				c += set(i);
			}
		}
		memo[n] = c;
		return c;

}
int main() {
	ifstream in("C://Users//86133//Desktop/input.txt");
	ofstream out("C://Users//86133//Desktop/output.txt");

	int n;
	in >> n;
	out << set(n);
	in.close();
	out.close();
	return 0;

}

input.txt:

 

output.txt:

四、代码分析 

set是递归函数，返回结果就是数n对应的半数集元素个数。数组memo用来储存set函数运算得到的结果，当下次输入参数相同时可以直接从数组中索取不必再重复进行复杂的递归，减少运算次数。

时间复杂度分析：

假设输入数字是n，时间复杂度取决于计算memo[i]的次数。将memo[i]看成一个状态，计算memo[n]时，memo[1]到memo[n/2]会被计算一次，因此共有O(n/2)个状态需要计算。对于每个状态，计算memo[i]的时间复杂度为O(i/2)，因此总时间复杂度为O(n^2)。

五、总结

1. 在计算机算法设计与分析中，递归技术是十分有用的。使用递归技术往往使函数的定义和算法的描述更加便捷且易于理解。
2. 当然，递归算法存在一个很大的问题在于很容易出现子问题重复计算，增大了时间复杂度。可以考虑增设一个记忆数组用来储存函数运算得到的结果，等下次输入参数相同时直接从记忆数组中获取，从而大大减少计算量。也可以采用动态规划的算法优化半数集问题。