动态规划法解决编程之美1.16节：24点问题

三，分支限界法求解 原代码出自 http://blog.csdn.net/tianshuai11/article/details/7713640，作者在main()中调用Fork((1<<N)-1);//开始1111 表示四个数，不易看出动态规划的过程（记录并利用前面已经求解过的子问题），我稍作变动，改为for (i = 1; i <= ((1 << N) - 1); i++)  Fork(i); 即可很直观的看出每次调用都是利用了前面已求解且已经记录的子问题

#include "stdafx.h"

#include <iostream>
#include <set>
#include <string>
#include <cmath>
using namespace std;

#define N	4	// 4张牌，可变
#define RES	24	// 运算结果为24，可变
#define EPS 1e-6

struct Elem
{
	Elem(double r, string& i):res(r),info(i){}
	Elem(double r, char* i):res(r),info(i){}
	double res;	// 运算出的数据
	string info; // 运算的过程
	bool operator<(const Elem& e) const
	{
		return res < e.res; // 在set的红黑树插入操作中需要用到比较操作
	}
};

int A[N];	// 记录N个数据
// 用二进制数来表示集合和子集的概念，0110表示集合包含第2,3个数
set<Elem> vset[1<<N];	// 包含4个元素的集合共有16个子集0-15

set<Elem>& Fork(int m)
{
	// memo递归
	
	
	if (vset[m].size())
	{
		
		return vset[m];
	}
	for (int i=1; i<=m/2; i++)
		if ((i&m) == i)
		{
			set<Elem>& s1 = Fork(i);
			set<Elem>& s2 = Fork(m-i);
			set<Elem>::iterator cit1;
			set<Elem>::iterator cit2;
			// 得到两个子集合的笛卡尔积，并对结果集合的元素对进行6种运算
			for (cit1=s1.begin(); cit1!=s1.end(); cit1++)
				for (cit2=s2.begin(); cit2!=s2.end(); cit2++)
				{
					string str;
					str = "("+  cit1->info + "+" + cit2->info + ")";
					vset[m].insert(Elem(cit1->res + cit2->res, str));

					str = "(" + cit1->info + "-" + cit2->info + ")";
					vset[m].insert(Elem(cit1->res - cit2->res, str));

					str = "("+ cit2->info + "-"+cit1->info+")";;
					vset[m].insert(Elem(cit2->res - cit1->res, str));

					str = "(" + cit1->info + "*" + cit2->info + ")";
					vset[m].insert(Elem(cit1->res * cit2->res, str));

					if (abs(cit2->res) > EPS) 
					{
						str = "(" + cit1->info + "/" + cit2->info + ")";
						vset[m].insert(Elem(cit1->res / cit2->res, str));
					}

					if (abs(cit1->res) > EPS)
					{
						str = "(" + cit2->info + "/" + cit1->info + ")";
						vset[m].insert(Elem(cit2->res / cit1->res, str));
					}
				}
		}
		return vset[m];
}

int main()
{
	int i;
	for (i=0; i<N; i++)
		cin >> A[i];
	// 递归的结束条件
	for (i=0; i<N; i++)
	{
		char str[10];
		sprintf(str,"%d",A[i]);
		vset[1<<i].insert(Elem(A[i],str));
	}
	for (i = 1; i <= ((1 << N) - 1); i++)
	{
		printf("i= %d\n", i);
		Fork(i);
	}
	
	//Fork((1<<N)-1);//开始1111 表示四个数 
	// 显示算出24点的运算过程
	set<Elem>::iterator it;
	for (it=vset[(1<<N)-1].begin(); 
		it!=vset[(1<<N)-1].end(); it++)
	{
		if (abs(it->res-RES) < EPS)
			cout << it->info << endl;
	}

	system("pause");
}