greedy算法计算金字塔分数问题

         金字塔分数问题：将一个真分数拆分成数个分子为1的分数之和。

         解决该问题的greedy方法是每次从剩余的部分选择一个分子为1的最大的分数，拆分成的金字塔所有数的特点是分母越来越大，形成一个金字塔的形状。

     （1）如何每次从剩余的部分都能选择一个最大的分数，这个可以进行论证如下：

证明：原真分数为a/b

若接下来选择的一个分数为1/c,且1/c是最接近a/b的分数

则根据题意需要满足以下条件：

a/b>=1/c；

则 c>=a/b;

又因为需要1/c最大，所以分母c必须最小，所以c必须为大于a/b的最小整数，得证。

其中每次剩余部分为。

     （2）公式计算如下式（1）所示

其中， a0=a, b0=b

则                             ................................    式（1）

 

最后结果为

 

   （3）为什么最终的结果a/b一定能够等于数个1/ci的和，是否会出现无限不循环的情况，即中i等于正无穷。但是只要满足某次循环中的分子为ac-b=1，就可终止循环，那么如何证明一定存在ac-b=1呢？有待证明？

  （4）算法的实现如下：

根据公式（1）可同过一个while循环进行循环迭代计算，算法如下：

   

std::vector<int> GreedyStrategy::egyptianFraction(map<int, int>num)
{
	map<int, int>::iterator iter = num.begin();                   //本质上是赋值迭代器以头地址
	vector<int> result;
	int b = iter->first;                                          //分子
	int a = iter->second;                                         //分母
	float c;
	int d;
	while (a > b)                                //继续递归
	{
		c = a / (b*1.0);
		d = a / b;
		if (c > d)
		{
			c = d + 1;
			result.push_back(c);
			b = b*c - a;
			a = a*c;
		}
		else
		{
			result.push_back(d);
			break;
		}
	}
	return result;
}

     注：最后返回的是所有分母的一个vector序列。