动态规划01-钢条切割问题

动态规划(dynamic programming)与分治算法相似,都是通过组合子问题的解来求解原问题。区别在于,分治算法是将原问题划分为互不相交的子问题,递归求解子问题,再把它们的解组合起来,求出原问题的解;动态规划应用于子问题重叠的情况,即不同的子问题具有公共的子问题(子问题的求解是递归进行的,将其划分为更小的子子问题),在这种情况下,分治算法会反复求解那些公共的子问题,而动态规划算法对每个子子问题只求解一次,将其保存在一个表格中,从而无需每次都求解一个子问题都重新计算,查表即可。
动态规划通常用来求解最优化问题,通常按以下4个步骤来设计一个动态规划算法:
(1)刻画一个最优解的结构特征
(2)递归地定义最优解的值
(3)计算最优解的值,通常采用自低向上的方法
(4)利用计算出的信息构造一个最优解
我们应用动态规划的第一个例子是钢条切割的问题:

长度i12345678910
价格1589101717202430

钢条切割问题是这样的:
给定一段长度为n英寸的钢条和一个价格表 pi ,求钢条的切割方案,使得销售收益 rn 最大。
问题定义为:

rn=max(pn,r1+rn1...,rn1+r1))

注意到,为了求解规模为n的原问题,我们先求解形式完全一样,但规模更小的子问题,我们通过组合两个子问题的最优解,并在所有可能的两段切割方案中选取组合收益最大者,构成原问题的最优解,我们称此满足最优子结构性质:问题的最优解由相关子问题的最优解组合而成,而这些子问题可以独立求解。

直接的求解方法为:

CUT-ROD(p,n)
  if n==0
    return 0
  q=INT_MIN
  for i=1 to n
    q=max(q,p[i]+CUT-ROD(p,n-i))
  return q

这种直接的方法的问题在于对相同的子问题进行了多次求解,运行时间为 2n1
动态规划方法会仔细安排求解顺序,对每个子问题只求解一次,并将其结果保存下来,如果随后再次需要此子问题的解,只需查找保存的结果,而不必重复计算。一般,我们采用自底向上(bottom-up)的方法

vector<int> steelCut(vector<int> &p,int n){
    vector<int> ret(n+1,0);
    for(int j=1;j<=n;j++){
        int q=INT_MIN;
        for(int i=1;i<=j;i++)
            q=std::max(q,p[i]+ret[j-i]);
        ret[j]=q;
    }
    return ret;
}

动态规划最最重要的两个问题是状态和状态转移方程!!状态ret[i]表示长度为i的钢条的最大效益!
状态转移方程就是

q=std::max(q,p[i]+ret[j-i])

其中ret[j-i]中保存的就是比j规模小的问题的解,这样可以避免重复求解!

总结:
动态规划问题的两个性质是
1.为了求解规模为n的原问题,我们先求解形式完全一样,但规模更小的子问题,我们通过组合两个子问题的最优解,并在所有可能的两段切割方案中选取组合收益最大者,构成原问题的最优解,我们称此满足最优子结构性质
2.子问题存在重叠,可通过自底而上求解子问题并保存结果来避免重复计算

动态规划对应最重要的两个关键点是:
1.状态方程
2.状态转移方程

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
动态规划是解决钢条切割问题的有效方法。下面是解决该问题的步骤: 1.定义问题:假设我们有一段长度为n的钢条和一个价格表pi(i=1,2,3….),其中pi表示长度为i英寸的钢条的价格。我们需要找到一种切割方案,使得切割后的钢条总价格最大。 2.找到最优子结构:假设我们已经找到了长度为n的钢条的最优切割方案,即rn。那么rn必然由一个或多个长度小于n的钢条的最优切割方案组成。因此,问题的最优解包含了子问题的最优解,具有最优子结构性质。 3.确定状态:我们可以用r[n]表示长度为n的钢条的最大收益。因此,我们需要求解的就是r[n]。 4.确定状态转移方程:假设我们已经知道了r, r, ..., r[n-1]的值,我们需要求解r[n]。我们可以枚举钢条的第一段切割位置,假设第一段切割位置为i,则有以下状态转移方程: r[n] = max(p[i] + r[n-i]),其中1<=i<=n 5.确定边界条件:当钢条长度为0时,收益为0,即r=0。 6.计算最优解:根据状态转移方程和边界条件,我们可以使用自底向上的方法计算出r[n]的值,从而得到最优解。 下面是Python代码实现: ```python def cut_rod(p, n): r = [0] * (n + 1) for j in range(1, n + 1): q = -1 for i in range(1, j + 1): q = max(q, p[i] + r[j - i]) r[j] = q return r[n] # 示例 p = [0, 1, 5, 8, 9, 10, 17, 17, 20, 24, 30] n = 4 print("长度为{}的钢条的最大收益为{}".format(n, cut_rod(p, n))) ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值