动态规划经典案例

对于k个阶段的动态规划问题，从第k个阶段到第1阶段的求解过程称为逆序解法。状态转移方程f（s）表示：从状态s到终点E的最短路径长度，用k表示阶段，则状态方程：

#include <string.h>
#include <map>
using namespace std;
#define MAX 21
#define INF 0x3f3f3f3f
//问题表示
int n;							//顶点个数
int start;						//起点编号
int endz;						//终点编号
int c[MAX][MAX];				//存放边长度
int nextz[MAX];					//存放最短路径上当前顶点的后继顶点
map<int,char *> vname;			//存放编号对应的顶点名称

int dp[MAX];
int Count=1;					//计算步骤
void Init()						//初始化
{
	n=10;
	start=0;
	endz=9;
	memset(c,INF,sizeof(c));
	for (int i=0;i<n;i++)		//初始化为-1
		dp[i]=-1;
	for (int j=0;j<n;j++)
		c[j][j]=0;
	c[0][1]=2; c[0][2]=4; c[0][3]=3;
	c[1][4]=7; c[1][5]=4;
	c[2][4]=3; c[2][5]=2; c[2][6]=4;
	c[3][4]=6; c[3][5]=2; c[3][6]=5;
	c[4][7]=3; c[4][8]=4;
	c[5][7]=6; c[5][8]=3;
	c[6][7]=3; c[6][8]=3;
	c[7][9]=3;
	c[8][9]=4;
	vname[0]="A";
	vname[1]="B1";	vname[2]="B2";	vname[3]="B3";
	vname[4]="C1";	vname[5]="C2";	vname[6]="C3";
	vname[7]="D1";	vname[8]="D2";
	vname[9]="E";
}
int f(int s)								//动态规划问题的逆序解法
{
	if (dp[s]!=-1) return dp[s];
	if (s==endz)								//找到终点
	{
		dp[s]=0;
		printf("  (%d) f(%s)=0\n",Count++,vname[s]);
		return dp[s];
	}
	else
	{
		int cost,mincost=INF,minj;
		for (int j=0;j<n;j++)				//查找顶点s的后继顶点
		{
			if (c[s][j]!=0 && c[s][j]!=INF)	//顶点s到j有边
			{
				cost=c[s][j]+f(j);			//先求出后继顶点的f值
				if (mincost>cost)
				{
					mincost=cost;
					minj=j;
				}
			}
		}
		dp[s]=mincost;
		nextz[s]=minj;
		printf("  (%d) f(%s)=c(%s,%s)+f(%s)=%d, ",Count++,vname[s],vname[s],vname[minj],vname[minj],dp[s]); 
		printf("next(%s)=%s\n",vname[s],vname[minj]);
		return dp[s];
	}
}
void main()
{
	Init();
	printf("%s→%s求解过程\n",vname[endz],vname[start]);
	f(start);
}