1D1D动态规划优化初步.docxVIP

PAGE PAGE # 1D/1D动态规划优化初步 所谓1D/1D动态规划，指的是状态数为 0(n)，每一个状态决策量为 0(n)的动态规划方程。 直接求解的时间复杂度为 0( n2)，但是，绝大多数这样的方程通过合理的组织与优化都是可 以优化到0(nlogn)乃至0(n)的时间复杂度的。这里就想讲一讲我对一些比较初步的经典的 优化方法的认识。 本文中不想进行过多的证明与推导，主要想说明经典模型的建立、转化与求解方法。 由于本人认识与水平相当有限，如果出现什么错误与疏漏， 还请大牛多多指正。 另外，也 希望大牛们更多地向我们介绍一下有关动态规划优化的更深入的东西。 本文中使用两种方式表示一个函数： f(x)与f[x]，用方括号表示的函数值可以在规划之前 全部算出(常量)，而用圆括号表示的函数值必须在规划过程中计算得到 (变量)。无论是什 么函数值一经确定，在以后的计算中就不会更改。 经典模型一： x J f (x)二 miin{ f(i) w[i,x]} i丄 相信这个方程大家一定是不陌生的。另外，肯定也知道一个关于决策单调性的性质: 假如用k(x)表示状态x取到最优值时的决策，则决策单调性表述为： 切兰j k( i)兰k( j)当且仅当： 切兰j, w i, j]+ w十1, j W[十1,j+] w [j对于这个性质的证明读者可以在 任意一篇讲述四边形不等式的文章中找到，所以这里不再重复。而且，从实战的角度来看， 我们甚至都不需要验证 w函数的这个性质，最经济也是最可靠的方法是写一个朴素算法打 出决策表来观察(反正你总还是要对拍) 。当然，有的时候题目要求你做一点准备工作，去 掉一些明显不可能的决策，然后在应用决策单调性。这是上述性质也许会有点用处。 正如前文中所述，我们关注的重点是怎样实现决策单调性。 有了决策单调性，怎样高效地 实现它呢？很容易想到在枚举决策的时候，不需要从 1开始，只要从k(x-1)开始就可以了， 但这只能降低常数，不可能起到实质性的优化。 另一种想法是从 k(x-1)开始枚举决策更新f(x)，一旦发现决策u不如决策u+1来得好，就 停止决策过程，选取决策 u作为f(x)的最终决策。这样时间是很大提高了，但可惜是不正确 的。决策单调性并没有保证 f(j)+w[j,x]有什么好的性质，所以这样做肯定是不对的。 刚才我们总是沿着“ f(x)的最优决策是什么”这个思路进行思考，下面我们换一个角度， 思考对于一个已经计算出来的状态 f(j) j “ f(j)能够更新的状态有哪些”。这样，每一步过程中 某些状态的决策可能不是最优的，但是当算法结束的时候所有状态对应的决策一定是最优 的。 一开始，只有f(1)的函数值被计算出来，于是所有状态的当前最优决策都是 1。 111 现在，显然f(2)的值已经确定了：它的最有决策只能是 1。我们用决策2来更新这个决策 表。由于决策单调性，我们知道新的决策表只能有这样的形式： 222222 这意味着我们可以使用二分法来查找“转折点” ，因为如果在一个点 x上，如果决策2更 好，则所有比x大的状态都是决策 2更好；如果x上决策1更好，则所有比x小的状态都是 决策1更好。 现在决策1和决策2都已经更新完毕，则 f(3)业已确定，现在用决策 3来更新所有状态。 根据决策单调性，现在的决策表只能有以下 2种类型： 333333 而这样的决策表示绝对不会出现的： 22222，不可能。 那么，我们的更新算法就是： 1、 考察决策2的区间[b,e]的b点上是否决策3更优，如果是，则全部抛弃决策 2,将此 区间划归决策3;如果否，则在决策 2的区间[b,e]中二分查找转折点。 2、 如果第1问的回答是“是”，则用同样的方法考察决策 1。 推演到这一步，相信决策单调性的实现算法已经明了了： 使用一个栈来维护数据，占中的 每一个元素保存一个决策的起始位置与终了位置， 显然这些位置相互连接且依次递增。 当插 入一个新的决策时，从后到前扫描栈，对于每一个老决策来说，做这样两件事： 1、 如果在老决策的起点处还是新决策更好，则退栈，全额抛弃老决策，将其区间合并至 新决策中，继续扫描下一个决策。 2、 如果在老决策的起点处是老决策好，则转折点必然在这个老决策的区间中； 二分查找 之，然后新决策进栈，结束。 由于一个决策出栈之后再也不会进入，所以均摊时间为 0(1)，但是由于二分查找的存在， 所以整个算法的时间复杂度为 0(nlogn)。 下面我们来看两个例题。 例题1:玩具装箱。 题目来源：湖南省选 2008。 题目大意：有n个玩具需要装箱，每个玩具的长度为 c[i]，规定在装箱的时候， 必须严格按 照给出的顺序进行，并且同一个箱子中任意两个玩具之间必须且只能间隔一个单位长度， 换 句话说，如