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

1、1D1D动态规划优化初步1D/1D动态规划优化初步所谓1D/1D动态规划，指的是状态数为 0(n)，每一个状态决策量为 0(n)的动态规划方程。 直接求解的时间复杂度为 0( n2)，但是，绝大多数这样的方程通过合理的组织与优化都是可 以优化到0(nlogn)乃至0(n)的时间复杂度的。这里就想讲一讲我对一些比较初步的经典的 优化方法的认识。本文中不想进行过多的证明与推导，主要想说明经典模型的建立、转化与求解方法。由于本人认识与水平相当有限，如果出现什么错误与疏漏， 还请大牛多多指正。 另外，也希望大牛们更多地向我们介绍一下有关动态规划优化的更深入的东西。本文中使用两种方式表示一个函数： f(

2、x)与fx，用方括号表示的函数值可以在规划之前全部算出(常量)，而用圆括号表示的函数值必须在规划过程中计算得到 (变量)。无论是什么函数值一经确定，在以后的计算中就不会更改。X_J经典模型一： f(X)二 min f (i) wi,x相信这个方程大家一定是不陌生的。另外，肯定也知道一个关于决策单调性的性质：假如用k(x)表示状态x取到最优值时的决策，则决策单调性表述为：-i空j, k()乞k( j)当且仅当：切兰j w i, j+ w十1, j 作 w十l,j+ w ，对于这个性质的证明读者可以在任意一篇讲述四边形不等式的文章中找到，所以这里不再重复。而且，从实战的角度来看， 我们甚至都不需要

3、验证 w函数的这个性质，最经济也是最可靠的方法是写一个朴素算法打出决策表来观察(反正你总还是要对拍) 。当然，有的时候题目要求你做一点准备工作，去掉一些明显不可能的决策，然后在应用决策单调性。这是上述性质也许会有点用处。正如前文中所述，我们关注的重点是怎样实现决策单调性。 有了决策单调性，怎样高效地实现它呢？很容易想到在枚举决策的时候，不需要从 1开始，只要从k(x-1)开始就可以了，但这只能降低常数，不可能起到实质性的优化。另一种想法是从 k(x-1)开始枚举决策更新f(x)，一旦发现决策u不如决策u+1来得好，就 停止决策过程，选取决策 u作为f(x)的最终决策。这样时间是很大提高了，但可

4、惜是不正确 的。决策单调性并没有保证 f(j)+wj,x有什么好的性质，所以这样做肯定是不对的。刚才我们总是沿着“ f(x)的最优决策是什么”这个思路进行思考，下面我们换一个角度，思考对于一个已经计算出来的状态 f(j) j “ f(j)能够更新的状态有哪些”。这样，每一步过程中某些状态的决策可能不是最优的，但是当算法结束的时候所有状态对应的决策一定是最优 的。一开始，只有f(1)的函数值被计算出来，于是所有状态的当前最优决策都是 1。111111*111111111111111111111111111111111111111111111111现在，显然f(2)的值已经确定了：它的最有决策只能

5、是 1。我们用决策2来更新这个决策表。由于决策单调性，我们知道新的决策表只能有这样的形式：111111*1111111111111222222222222222222222222222222这意味着我们可以使用二分法来查找“转折点” ，因为如果在一个点 x上，如果决策2更好，则所有比x大的状态都是决策 2更好；如果x上决策1更好，则所有比x小的状态都是 决策1更好。现在决策1和决策2都已经更新完毕，则 f(3)业已确定，现在用决策 3来更新所有状态。 根据决策单调性，现在的决策表只能有以下 2种类型：111111*1111111111111111222222222222222222333333

6、33333 111111*1111111113133333333333333333333333333333333333而这样的决策表示绝对不会出现的：111111*333333333333333322222222222222222222222222222，不可能。那么，我们的更新算法就是：1、 考察决策2的区间b,e的b点上是否决策3更优，如果是，则全部抛弃决策 2,将此 区间划归决策3;如果否，则在决策 2的区间b,e中二分查找转折点。2、 如果第1问的回答是“是”，则用同样的方法考察决策 1。推演到这一步，相信决策单调性的实现算法已经明了了： 使用一个栈来维护数据，占中的每一个元素保存一个

7、决策的起始位置与终了位置， 显然这些位置相互连接且依次递增。 当插入一个新的决策时，从后到前扫描栈，对于每一个老决策来说，做这样两件事：1、 如果在老决策的起点处还是新决策更好，则退栈，全额抛弃老决策，将其区间合并至 新决策中，继续扫描下一个决策。2、 如果在老决策的起点处是老决策好，则转折点必然在这个老决策的区间中； 二分查找 之，然后新决策进栈，结束。由于一个决策出栈之后再也不会进入，所以均摊时间为 0(1)，但是由于二分查找的存在，所以整个算法的时间复杂度为 0(nlogn)。下面我们来看两个例题。例题1:玩具装箱。题目来源：湖南省选 2008。题目大意：有n个玩具需要装箱，每个玩具的长

8、度为 ci，规定在装箱的时候， 必须严格按照给出的顺序进行，并且同一个箱子中任意两个玩具之间必须且只能间隔一个单位长度， 换句话说，如果要在一个箱子中装编号为 ij的玩具，则箱子的长度必须且只能是j丨=j -i ck，规定每一个长度为I的箱子的费用是P =(1 - L)2，其中L是给定的 k=t一个常数。现在要求你使用最少的代价将所有玩具装箱，箱子的个数无关紧要。分析：本题可以很轻松地列出一个 1D1D的动态规划方程：X jf (x)二 mi n i() w i，其中 wi, j =( j v ck - L)2。i k甘不难验证这个方程式满足决策单调性的，于是我们可以直接套用上文中的方法进行优

9、化， 时间复杂度为O(nlogn)。例题2 :土地购买题目来源：USACO Mo nthly, March, 2008, Gold题目大意：有 N块土地需要购买，每块土地都是长方形的，有特定的长与宽。你可以一次性购买一组土地，价格是这组土地中长的最大值乘以宽的最大值。 比方说一块5*3的土地和一块2*9的土地在一起购买的价格就是 9*3。显然，怎样分组购买土地是一门学问，你的任务就是设计一种方案用最少的钱买下所有的土地。分析：将所有土地按照长度降序排列， 依次检索，则当前土地的长度必然在上一块土地之内，我们只需要考虑宽度就可以了。 而在宽度的问题上，当前土地的行为只能是这样： 和前面若干块土地

10、绑定；同时这些绑定的土地和他们前后的土地分离。 这样很容易得出状态转移方程：n nf(n)=吧可( maxwi)*lk +1 + f(k)这个方程还不能满足决策单调性，下面我们试图再做一下简化。如果将每一个土地的尺寸看成是一个二维坐标的话， (如下图)其中不难看出，红色点完全可以忽略，这些点 (x,y)必然满足一个性质：存在点 (x y)同时满足x =且y = y这样它就能被一个组完全覆盖。 这些被忽略的点可以通过一次线形的扫描得出。n 二f(n )=叫% n *yk 1 f(k)这个转移方程就是标准的决策单调性了，读者可以通过 w函数的性质直接证明它。然后,就用上文中的方法在 0(nlogn

11、)时间内求解。以上两个例子都是决策单调性的直接应用。 其中第二个例子稍微复杂一些，如果不忽略那些“肯定无用”的决策，不对数据进行有序化，则方程是不满足决策单调性的。这也就提醒 我们在做这一类题目的时候不能钻牛角尖死做，还得灵活一点。另外，决策单调性提供的只是 0(nlogn)的算法，事实上上面两个例题的最佳算法都是 0(n)的，在后文中我们将详细介绍另外一种经典模型， 并且试图将这两个规划方程通过数学变换转向另一个模型。F面我们来看一类特殊的 w函数： j : k, wi, j w j, k wi ,k，显然，这一类函数都是满足决策单调性的。 但是不同的是，由于这一类函数的特殊性， 他们可以用

12、一种更加简洁也更加有借鉴意义的方法解决。由于w函数满足-i冬j ：k,wi,j wj,k =wi,k，我们总是可以找到一个特定的一元函数wx，使得-i岂j, wi, j = w j 一 wx，这样，假设状态f(x)的某一个决策是k，有:，其中 g(k)二 f (k) w1,k。f(x)二 f(k) wk,x = f (k) wxwk二 g(k) wx -w1,这样我们发现：一旦 f(k)被确定，相应地 g(k)也被确定，更加关键的是， 无论k值如何，wx-w 1总是一个常数。换句话说，我们可以把方程写成下述形式：x-1 xf (x ming(k) w1,x。不难发现这个方程是无聊的，因为 mi

13、n g(k)我们可以用一个变量“打擂台”直接存储；但是，如果在 k的下界上加上一个限制，那这个方程就不是很 无聊了。于是，我们就得到了另一个经典模型。经典模型二：x Jf (x)二 min g(k) wx,其中，bx随 x 不降。k生x这个方程怎样求解呢？我们注意到这样一个性质：如果存在两个数 j, k ,使得j = k，而且g(k) = g(j)，则决策j是毫无用处的。因为根据 bx单调的特性，如果j可以作为合法决 策，那么k 一定可以作为合法决策，又因为 k比j要优，(注意：在这个经典模型中， “优”是绝对的，是与当前正在计算的状态无关的) ，所以说，如果把待决策表中的决策按照 k排序的话

14、，贝U g(k)必然是不降的。这样，就引导我们使用一个 单调队列 来维护决策表。对于每一个状态 f(x)来说，计算过程分为以下几步：1、 队首元素出队，直到队首元素在给定的范围中。2、 此时，队首元素就是状态 f(x)的最优决策，3、 计算g(x)，并将其插入到单调队列的尾部，同时维持队列的单调性(不断地出队，直 到队列单调为止)。重复上述步骤直到所有的函数值均被计算出来。不难看出这样的算法均摊时间复杂度是0(1)的。下面我们来看几个例题。例题 3: The Sound of Silence题目来源：Baltic Olympiad in In formatics 2007题目大意：给出一个 N

15、项的数列，如果对于一个连续的长度为 M的片段来说，片段内所有数中最大值与最小值的差不超过一个给定的常数 C,则我们称这样的片段是一个合法的片段。编程求出所有的合法片段的起始位置。分析：本题不难看出可以分解为两个子问题： 求所有片段的最大值以及求所有片段的最小值。而这两个任务实际上是一样的， 所以我们只需要求取所有的连续 M个数的片段中的最小值。这个任务有很多很多种对数级算法，其中用堆或者用静态最优二叉树都可以做到O(nlogm)，但是这题的0(n)算法还是不那么好想的。事实上，如果用 gx表示数列中第x个数的值，用f(x)表示以x作为结尾的有M个数的 连续片段的话，显然有：xf (x. ng

16、i，直接吻合经典模型二。套用算法，就可以在0(n)的时间内解决问题。(当然，本题还有一种别致的“窗口”算法，也漂亮地在 0(n)的时间内解决了问题，详细可以看官方的解题报告。这里引入本题的主要目的是在于佐证上文中讨论到的经典模型二)例题 4： Islands题目来源：1012008题目大意：给出一个具有 N个顶点的无向加权图，同时这个图中有且恰有 N条边。现在,对于这个图中的每一个连通分量， 求出其最长路径(权值和最大，一个节点在路径上最多只 能出现一次)。分析：当然，这个问题更多的是一个图论题。 但是在最后关键问题的处理上还是可以看到经典模型二的影子。首先，用BFS找出所有连通分量。然后，对

17、于一个连通分量来说，由于点数与边数相同， 因此必然构成 基环+外向树的结构。我们可以找出基环并确定所有外向树的结构。一条最长 路径有两种可能：完整地位于某一棵外向树中； 或者位于两棵外向树中， 其间通过基环的一段连接。第一种可能可以通过树形 DP解决，问题就在于第二种可能怎样处理。如果枚举两棵外向树，那就是 0( n2)，就不可以接受了。我们考虑破环为链，然后将链整体左移，制作一个副本。比方说，如果原来的环是：4 2，以1为首破环，得：1-2-3-4，然后制作副本，得 2-3-4-1-2-3-4，制作副本3的主要目的是使得对于每一个点的方程都有统一的形式， 使得环上所有片段都可以对应链上的一个

18、片段。这种情况下，用 gx表示x点上外向树上的最长下降路的长度， f(x)表示以该点为终点的总最长路径的长度，则有：x 4f(x) = maxgi gx disi,x，其中 w 函数即 w i, j d is ,i显j然满足 i立-n十Vi j wk* (卜 i显然满足i =x_mki j ： k,wi, j w j, k二wi,k，(注意w是一个与i和j无关的常量)经过适当的变化后可以转化为经典模型二。于是有限背包问题可以在 O(NP)的时间内解决，其中P是背包可能取到的最大价值。(其实换成重量也一样)，这也就是“背包十讲”中所说的那个单调 队列法。我们注意到，如果 mk=1的话，那么每一个

19、f(x)的决策量都是 O(1)，这没什么问题；但 如果mk = +R,有意味着什么呢？仔细观察可以发现，这实际上就拿掉了方程中的循环变量的下界，对应的是 f(x) =mi ng(k) w1,x这样的一个方程，这显然是很简单的，适k 二用单变量打擂台就可以解决了(尽管我们通常并不这样做) 。所以说，借助经典模型二，我们在一个更高的高度上统一了 0-1、有限、无限三大背包问题。下面我们再次来看一下例题 2土地购买中的那个方程：n .4f(n) =minf(k) xn* yk 1,我们来仔细地观察这个方程：f(k)是变量，yk+1是常 量，但不论怎么说，这两个量在以后的计算中都不会变化。而 x n是

20、一个比例系数，它与 k无关，只随着x的变化而变化。如果我们建立平面直角坐标系，以 f(k)作为横轴，yk+1作为纵轴，则每一个状态 f(k)都可以看作是该坐标系中的一个点。在求解状态 f(n)的过程中，我要求最小化：min P -x ky，其中x, y是我建立的直角坐标系中某一个点的坐标(表示一个决策) ，k就是方程中的xn，是只与n有关，而与决策无关的一个常量。这个最小化问题是什么呢？其实就是一个平面上的线性规划。我们把式子改写成：1 P 1y x ，就演变成了这样的一个问题：在一个直线簇 y x C中，选取一条直k k k线，使得这条直线过某个给出的数据点，同时 C要最小。既然问题变成了这

21、么有意思的线性规划问题， 就有必要进一步的研究， 看看是不时有更好的解法，这就导致了我们的另一个经典模型：x J经典模型三：f (x) = mnax* f (i) bx* g(i)。注意：这个模型写的比较抽象，其实它的涵盖范围是很广的。首先， ax, bx不一定要是常量，只要他们与决策无关， 都可以接受；另外，f(i)和g(i)不管是常量还是变量都没有关系， 只要他们是一个有最优的 f(x)决定的二元组就可以了。因此，为了方便描述，我们把这个模型写成下面这个形式：n -1f (n) = miin an* x(i) bn* y(i)，i日其中，x(i)，y(i)都是可以在常数时间内通过 f(i)

22、唯一决定的二元组。这个经典模型怎样转化求解呢？前文说过，这样的模型的求解与平面上的线性规划有关，我们以x(i)为横轴，y(i)为纵轴建立平面直角坐标系，这样一个状态 f(i)所决定的二元组就可以用坐标系中的一个点表示。然后，我们的目标是：a Pm i rP 二 a x b,y其中 a = an, b = bn，化成：y x ，假设 b 0 (反之亦 b b然)，则我们的任务是使得这条直线的纵截距最小。可以想象有一组斜率相同的直线自负无 穷向上平移，所碰到的第一个数据点就是最优决策。这个时候，有一个重要的性质，那就是： 所有最优决策点都在平面点集的凸包上。 基于这个事实，我们可以开发出很多令人满

23、意的算法。这时，根据直线斜率与数据点分布的特征，可以划分为两种情况：(具体的单调性视最优化目标情况一：决策直线的斜率与二元组的横坐标同时满足单调性。 的性质而定)这样的模型是比较好处理的， 因为这个时候由于斜率变化是单调的， 所以决策点必然在凸壳上单调移动。我们只需要维护一个单调队列和一个决策指针，每次决策时作这样几件事：1、 决策指针(也就是队首)后移，直至最佳决策点。2、 进行决策。3、 将进行决策之后的新状态的二元组加入队尾， 同时作Graham-Scan式的更新操作维护 凸壳。(注意此时当前指针所在二元组有可能被抛弃)算法的时间复杂度为 0(n)情况二：没有任何限制。这时问题的解决就比

24、较困难了。 显然，决策点还是应该在凸壳上。我们不妨考虑一个单调减的凸壳，这个凸壳上点与点之间的连线必然满足这样的性质： 斜率单调减。通过细致的观察我们可以发现，对于一个给定的斜率 k来说，对应的直线簇中具有最大纵截距的直线通 过的决策点必然满足这样的性质：该点两侧的边的斜率 k1, k2满足k1 _ k _ k2。这样，我们就可以通过二分查找来确定最佳的决策点。但是，在插入数据点的过程中，我们遇到的麻烦可能更大。 首先，肯定是二分查找确定横坐标的插入点，然后对两侧分别进行 Graham维护凸性。但接下来的问题就严重了：在维护凸形的过程中我们肯定删掉了一些点，怎样重新得到一个完整的凸壳决策表呢？

25、使用 move 是一个折中的办法，但是这与理论的时间复杂度分析根本无益。完美的解决方法是应用平衡二叉树。 我们以横坐标为关键字建立平衡二叉树， 这样查找和插入的过程都可以在 O(logn)时间内完成。当我们做 Graham维护时，首先将新数据点 Splay到根节点，此时剩下的节点必然分居左子树和右子树。然后， 我们以左子树为例， 后序遍历依次查找节点，直至查找到一个满足凸形的节点。将这个节点 Splay到根节点的左孩子，然后删掉这个节点的右孩子。这样的算法的时间复杂度是 O(nlogn)，但是实现起来非常复杂。下面我们来看几个例题。例题6:玩具装箱的线性算法。例题1中玩具装箱的动态规划方程为：

26、X xf (x)= mi rfi ( X -(- ck - L 2,下面，我们试图通过数学变换将其变成经- kM典模型三。X为了简化计算，设 sumx -7 cx，则：4x Jf (x)二 min f (i) (xi -1L sumx -sumi)2= mn i(SUmx|Mx L 1)sUm#i)不妨设ax二sumx x-丨-1,bi二sumi i，显然这两个量都是常量，则:X-1 2f(x)二 mp f (i)(ax -bi) X二 m f (i) a x b i -2axbix -4 ，二 mp f (i) b2i2axbi a2x然后问题就明朗了，设平面直角坐标系中 x(ibi, y(

27、i f(i) b2i，则问题变成：min P =2ax y，其对应的线性规划的目标直线为 y = -2ax P。回顾定义不难看出，ax随着x的增大而增大，x(i)也随着i的增大而增大。因此，问题 中直线斜率单调减，数据点横坐标单调增，符合经典模型三种的情形 1。使用单调队列维护凸壳可以在O(n)的时间内解决本题。例题七：货币兑换。题目来源：NOI2007题目大意：有3种货币体系：人民币， A券，B券，其中A券与B券的价格在每一天都是不同的。在某一天 D，你可以做3件事情：1、 如果你的手头上有 A券或B券，你可以将它们都按照当天的价格换成人民币。2、 如果你的手头上有人民币， 你可以将它们按照

28、一个特定的比例 Rate并以照当天的价格换成A券和B券(就是说你兑换的 A券和B券的价值比是 Rate)3、什么也不做。一开始你有一些人民币，请你通过最佳的操作方式在最后一天结束的时候手头上握有最多的 人民币。分析：试题中的 Hint已经告诉我们，如果我们想买进人民币，就必须全额买进；如果我们 想卖出人民币，就必须全额卖出。 由于不管是在哪一天，人民币总是越多越好；我们用 f(n)表示到了第n天最多可能持有的人民币数量，用 x(i)和y(i)表示在第i天，用最多的钱能够换成的A券和B券。(注意：由于 Rate确定，兑换金额确定，所以 A券和B券的数量是唯 一确定的)，我们有：n Jf (n) =maxan* x(i) bn* y(i)，其中an和bn代表a券和b券在第n天的牌价。 这个方程式符合经典模型三中的情形二。 所以，我们应该使用一个平衡树来维护凸形。 时间复杂度是0(nlogn)。由于数据弱，所以这道题目用 move也是可以搞定的。这篇文章中，我们着重讨论了这样三类经典模型的