Chinaunix首页 | 论坛 | 博客
  • 博客访问: 2877395
  • 博文数量: 471
  • 博客积分: 7081
  • 博客等级: 少将
  • 技术积分: 5369
  • 用 户 组: 普通用户
  • 注册时间: 2012-01-04 21:55
文章分类

全部博文(471)

文章存档

2014年(90)

2013年(69)

2012年(312)

分类: C/C++

2012-05-22 09:21:22

动态规划(Dynamic Programming,简称DP)是通过组合子问题的解来解决问题的。

注意这里的programming不是指编程,而是指一种规划

因为DP用的太广泛了,并且书上也花了很大的篇幅来讲这一章,所以我也准备把这一章分几篇来总结,并且我不按照书上的顺序来总结,也不会全部用书上的例子。

这一章首先给出一些基础的概念,然后给出一个最基础入门的DP问题,三角形求值问题。

DP适用于子问题不是独立的情况,这样如果用分治法,也会作许多重复的工作,而DP只对子问题求解一次,将其结果保存在一张表中,从而避免了每次遇到各个子问题时重新计算的情况。

比较分治法于动态规划的区别:

分治法:将问题划分为一些独立的子问题,递归的求解各子问题,然后合并子问题的解而得到原问题的解。

eg.

MERGE-SORT(A, p, r) 1 if p < r 2 then q ← |(p + r)/2| 3 MERGE-SORT(A, p, q) 4 MERGE-SORT(A, q + 1, r) 5 MERGE(A, p, q, r)动态规划适用于子问题不是独立的情况,也就是各子问题包含公共的子子问题,鉴于会重复的求解各子问题,DP对每个问只求解一遍,将其保存在一张表中,从而避免重复计算。

DP算法的设计可以分为四个步骤

①.描述最优解的结构。

②.递归定义最优解的值。

③.按自底而上的方式计算最优解的值。

④.由计算出的结果创造一个最优解。

下面来看看三角形求值这个问题:

将一个由N行数字组成的三角形,如图所以,设计一个算法,计算出三角形的由顶至底的一条路径,使该路径经过的数字总和最大。

sanjiaoxing

这是在我见过的DP题目中,算是最简单的了,我相信就算没有学过DP的也会。

将上图转化一下:

sanjiaoxing2

假设上图用map[][]数组保存。

令f[i][j]表示从顶点(1, 1)到顶点(i, j)的最大值。

则可以得到状态转移方程:

f[i][j] = max(f[i+1][j], f[i+1][j+1]) + map[i][j]

此题既适合自顶而下的方法做,也适合自底而上的方法,

当用自顶而下的方法做时,最后要在在最后一列数中找出最大值,

以下代码采用自底而上的方法做:

而用自底而上的方法做时,f[1][1]即为最大值。

所以我们将图2根据状态转移方程可以得到图3:

sanjiaoxing3

最大值是30.

点击(此处)折叠或打开

  1. /*
  2. Author: Tanky Woo
  3. Blog: www.WuTianQi.com
  4. Description: 动态规划之三角形求值问题
  5. */
  6. #include <iostream>
  7. using namespace std;
  8.  
  9. int map[6][6] =
  10. {
  11.     {0, 0, 0, 0, 0, 0},
  12.     {0, 7, 0, 0, 0, 0},
  13.     {0, 3, 8, 0, 0, 0},
  14.     {0, 8, 1, 0, 0, 0},
  15.     {0, 2, 7, 4, 4, 0},
  16.     {0, 4, 5, 2, 6, 5}
  17. };
  18.  
  19. int f[6][6];
  20.  
  21. int _max(int a, int b)
  22. {
  23.     if(a >= b)
  24.         return a;
  25.     return b;
  26. }
  27.  
  28. int main()
  29. {
  30.     memset(f, 0, sizeof(f));
  31.     for(int i=0; i<6; ++i)
  32.         f[5][i] = map[5][i];
  33.     for(int i=4; i>=1; --i)
  34.         for(int j=1; j<=i; ++j)
  35.             f[i][j] = _max(f[i+1][j], f[i+1][j+1]) + map[i][j];
  36.     for(int i=1; i<=5; ++i)
  37.     {
  38.         for(int j=1; j<=5; ++j)
  39.         {
  40.             cout.width(2);
  41.             cout << f[i][j] << " ";
  42.         }
  43.         cout << endl;
  44.     }
  45.     cout << f[1][1] << endl;

转自:

阅读(856) | 评论(0) | 转发(0) |
0

上一篇:集合类框架

下一篇:动态规划(2)基本入门

给主人留下些什么吧!~~