算法合集之《动态规划算法的优化技巧》.doc
《算法合集之《动态规划算法的优化技巧》.doc》由会员分享,可在线阅读,更多相关《算法合集之《动态规划算法的优化技巧》.doc(16页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、 福州第三中学 毛子青动态规划算法的优化技巧福州第三中学 毛子青关键词 动态规划、 时间复杂度、优化、状态摘要动态规划是信息学竞赛中一种常用的程序设计方法,本文着重讨论了运用动态规划思想解题时时间效率的优化。全文分为四个部分,首先讨论了动态规划时间效率优化的可行性和必要性,接着给出了动态规划时间复杂度的决定因素,然后分别阐述了对各个决定因素的优化方法,最后总结全文正文一、引言动态规划是一种重要的程序设计方法,在信息学竞赛中具有广泛的应用。使用动态规划方法解题,对于不少问题具有空间耗费大、时间效率高的特点,因此人们在研究动态规划解题时更多的注意空间复杂度的优化,运用各种技巧将空间需求控制在软硬件
2、可以承受的范围之内。但是,也有一部分问题在使用动态规划思想解题时,时间效率并不能满足要求,而且算法仍然存在优化的余地,这时,就需要考虑时间效率的优化。本文讨论的是在确定使用动态规划思想解题的情况下,对原有的动态规划解法的优化,以求降低算法的时间复杂度,使其能够适用于更大的规模。二、动态规划时间复杂度的分析使用动态规划方法解题,对于不少问题之所以具有较高的时间效率,关键在于它减少了“冗余”。所谓“冗余”,就是指不必要的计算或重复计算部分,算法的冗余程度是决定算法效率的关键。动态规划在将问题规模不断缩小的同时,记录已经求解过的子问题的解,充分利用求解结果,避免了反复求解同一子问题的现象,从而减少了
3、冗余。但是,动态规划求解问题时,仍然存在冗余。它主要包括:求解无用的子问题,对结果无意义的引用等等。下面给出动态规划时间复杂度的决定因素:时间复杂度=状态总数*每个状态转移的状态数*每次状态转移的时间1下文就将分别讨论对这三个因素的优化。这里需要指出的是:这三者之间不是相互独立的,而是相互联系,矛盾而统一的。有时,实现了某个因素的优化,另外两个因素也随之得到了优化;有时,实现某个因素的优化却要以增大另一因素为代价。因此,这就要求我们在优化时,坚持“全局观”,实现三者的平衡。三、动态规划时间效率的优化3.1 减少状态总数我们知道,动态规划的求解过程实际上就是计算所有状态值的过程,因此状态的规模直
4、接影响到算法的时间效率。所以,减少状态总数是动态规划优化的重要部分,本节将讨论减少状态总数的一些方法。1、改进状态表示状态的规模与状态表示的方法密切相关,通过改进状态表示减小状态总数是应用较为普遍的一种方法。例一、 Raucous Rockers 演唱组(USACO96)问题描述现有n首由Raucous Rockers 演唱组录制的珍贵的歌曲,计划从中选择一些歌曲来发行m张唱片,每张唱片至多包含t分钟的音乐,唱片中的歌曲不能重叠。按下面的标准进行选择:(1) 这组唱片中的歌曲必须按照它们创作的顺序排序;(2) 包含歌曲的总数尽可能多。输入n,m,t,和n首歌曲的长度,它们按照创作顺序排序,没有
5、一首歌超出一张唱片的长度,而且不可能将所有歌曲的放在唱片中。输出所能包含的最多的歌曲数目。(1n, m, t20)算法分析本题要求唱片中的歌曲必须按照它们创作顺序排序,这就满足了动态规划的无后效性要求,启发我们采用动态规划进行解题。分析可知,该问题具有最优子结构性质,即:设最优录制方案中第i首歌录制的位置是从第j张唱片的第k分钟开始的,那么前j-1张唱片和第j张唱片的前k-1分钟是前1.i-1首歌的最优录制方案,也就是说,问题的最优解包含了子问题的最优解。设n首歌曲按照写作顺序排序后的长度为long1.n,则动态规划的状态表示描述为:gi, j, k,0in,0jm,0kt,表示前i首歌曲,用
6、j张唱片另加k分钟来录制,最多可以录制的歌曲数目,则问题的最优解为gn,m,0。由于歌曲i有发行和不发行两种情况,而且还要分另加的k分钟是否能录制歌曲i。这样我们可以得到如下的状态转移方程和边界条件:当klongi,i1时:gi, j, k=maxgi-1,j,k-longi,gi-1,j,k当klongi,i1时:gi, j, k=maxgi-1,j-1,t-longi,gi-1,j,k规划的边界条件为:当0kt-b时: a=a+1; b=longi;规划的边界条件:gi,0=(0,0) 0in这样题目所求的最大值是:ans=maxk| gn, k(m-1,t)改进后的算法,状态总数为O(n
7、2),每个状态转移的状态数为O(1),每次状态转移的时间为O(1),所以总的时间复杂度为O(n2)。值得注意的是,算法的空间复杂度也由改进前的O(m*n*t)降至优化后的O(n2)。(程序及优化前后的运行结果比较见附件)通过对本题的优化,我们认识到:应用不同的状态表示方法设计出的动态规划算法的性能也迥然不同。改进状态表示可以减少状态总数,进而降低算法的时间复杂度。在降低算法的时间复杂度的同时,也降低了算法的空间复杂度。因此,减少状态总数在动态规划的优化中占有重要的地位。2、选择适当的规划方向 动态规划方法的实现中,规划方向的选择主要有两种:顺推和逆推。在有些情况下,选取不同的规划方向,程序的时
8、间效率也有所不同。一般地,若初始状态确定,目标状态不确定,则应考虑采用顺推,反之,若目标状态确定,而初始状态不确定,就应该考虑采用逆推。那么,若是初始状态和目标状态都已确定,一般情况下顺推和逆推都可以选用,但是,能否考虑选用双向规划呢?双向搜索的方法已为大家所熟知,它的主要思想是:在状态空间十分庞大,而初始状态和目标状态又都已确定的情况下,由于扩展的状态量是指数级增长的,于是为了减少状态的规模,分别从初始状态和目标状态两个方向进行扩展,并在两者的交汇处得到问题的解。上述优化思想能否也应用到动态规划之中呢?来看下面这个例子。例二、 Divide (Merc2000)问题描述有价值分别为1.6的大
9、理石各a1.6块,现要将它们分成两部分,使得两部分价值和相等,问是否可以实现。其中大理石的总数不超过20000。(英文试题详见附件)算法分析令S=(i*ai),若S为奇数,则不可能实现,否则令Mid=S/2,则问题转化为能否从给定的大理石中选取部分大理石,使其价值和为Mid。这实际上是母函数问题,用动态规划求解也是等价的。mi, j,0i6,0jMid,表示能否从价值为1.i的大理石中选出部分大理石,使其价值和为j,若能,则用true表示,否则用false表示。则状态转移方程为:mi, j=mi, j OR mi-1,j-i*k (0kai)规划的边界条件为:mi,0=true; 0i6若mi
10、, Mid=true,0i6,则可以实现题目要求,否则不可能实现。我们来分析上述算法的时间性能,上述算法中每个状态可能转移的状态数为ai,每次状态转移的时间为O(1),而状态总数是所有值为true的状态的总数,实际上就是母函数中项的数目。算法优化实践发现:本题在i较小时,由于可选取的大理石的价值品种单一,数量也较少,因此值为true的状态也较少,但随着i的增大,大理石价值品种和数量的增多,值为true的状态也急剧增多,使得规划过程的速度减慢,影响了算法的时间效率。另一方面,我们注意到我们关心的仅是能否得到价值和为Mid的值为true的状态,那么,我们能否从两个方向分别进行规划,分别求出从价值为
11、1.3的大理石中选出部分大理石所能获得的所有价值和,和从价值为4.6的大理石中选出部分大理石所能获得的所有价值和。最后通过判断两者中是否存在和为Mid的价值和,由此,可以得出问题的解。状态转移方程改进为:当i3时:mi, j=mi, j OR mi-1,j-i*k (1kai)当i3时:mi, j=mi, j OR mi+1,j-i*k (1kai)规划的边界条件为:mi,0=true; 0i7这样,若存在k,使得m3,k=true, m4,Mid-k=true,则可以实现题目要求,否则无法实现。(程序及优化前后的运行结果比较见附件)从上图可以看出双向动态规划与单向动态规划在计算的状态总数上的
12、差异。回顾本题的优化过程可以发现:本题的实际背景与双向搜索的背景十分相似,同样有庞大的状态空间,有确定的初始状态和目标状态,状态量都迅速增长,而且可以实现交汇的判断。因此,由本题的优化过程,我们认识到,双向扩展以减少状态量的方法不仅适用于搜索,同样适用于动态规划。这种在不同解题方法中,寻找共通的属性,从而借用相同的优化思想,可以使我们不断创造出新的方法。3.2 减少每个状态转移的状态数在使用动态规划方法解题时,对当前状态的计算都是进行一些决策并引用相应的已经计算过的状态,这个过程称为“状态转移”。因此,每个状态可能做出的决策数,也就是每个状态可能转移的状态数是决定动态规划算法时间复杂度的一个重
13、要因素。本节将讨论减少每个状态可能转移的状态数的一些方法。1、四边形不等式和决策的单调性例三、石子合并问题(NOI95)问题描述 在一个操场上摆放着一排n(n20)堆石子。现要将石子有次序地合并成一堆。规定每次只能选相邻的2堆石子合并成新的一堆,并将新的一堆石子数记为该次合并的得分。试编程求出将n堆石子合并成一堆的最小得分和最大得分以及相应的合并方案。算法分析这道题是动态规划的经典应用。由于最大得分和最小得分是类似的,所以这里仅对最小得分进行讨论。设n堆石子依次编号为1,2,.,n。各堆石子数为d1.n,则动态规划的状态表示为:mi,j,1ijn,表示合并di.j所得到的最小得分,则状态转移方
14、程和边界条件为:mi,j=0 i=j ij同时令si,j=k,表示合并的断开位置,便于在计算出最优值后构造出最优解。上式中的计算,可在预处理时计算,i=1.n;t0=0, 则:上述算法的状态总数为O(n2),每个状态转移的状态数为O(n),每次状态转移的时间为O(1),所以总的时间复杂度为O(n3)。算法优化当函数wi,j满足时,称w满足四边形不等式2。当函数wi,j满足wi,jwi,j 时称w关于区间包含关系单调。在石子归并问题中,令wi,j= ,则wi,j满足四边形不等式,同时由di0,ti0可知wi,j满足单调性。mi,j=0 i=j ij 对于满足四边形不等式的单调函数w,可推知由递推
15、式定义的函数mi,j也满足四边形不等式,即。这一性质可用数学归纳法证明如下:我们对四边形不等式中“长度”l=j-i进行归纳:当i=i或j=j时,不等式显然成立。由此可知,当l1时,函数m满足四边形不等式。下面分两种情形进行归纳证明:情形1:ii=jj。下面只讨论kj,kj的情况是类似的。情形1.1:kj,此时:情形2:iijy。下面只讨论zy,zy的情况是类似的。由izyj有:综上所述,mi,j满足四边形不等式。令si,j=maxk | mi,j=mi,k-1+mk,j+wi,j 由函数mi,j满足四边形不等式可以推出函数si,j的单调性,即si,jsi,j+1si+1,j+1, ij当i=j
16、时,单调性显然成立。因此下面只讨论ij的情形。由于对称性,只要证明si,jsi,j+1。令mki,j=mi,k-1+mk,j+wi,j。要证明si,jsi,j+1,只要证明对于所有ikkj且mki,jmki,j,有:mki,j+1mki,j+1。事实上,我们可以证明一个更强的不等式mki,j-mki,jmki,j+1-mki,j+1也就是: mki,j+mki,j+1mki,j+1+mki,j利用递推定义式将其展开整理可得:mk,j+mk,j+1mk,j+mk,j+1,这正是kkjj+1时的四边形不等式。综上所述,当w满足四边形不等式时,函数si,j具有单调性。于是,我们利用si,j的单调性,
17、得到优化的状态转移方程为:mi,j=0 i=j iy,下面仅讨论zy,zy的情况是类似的。由izyj有: 接着,我们用数学归纳法证明函数m也满足四边形不等式。对四边形不等式中“长度”l=j-i进行归纳:当i=i或j=j时,不等式显然成立。由此可知,当l1时,函数m满足四边形不等式。下面分两种情形进行归纳证明:情形1:ii=jj。下面只讨论kj,kj的情况是类似的。情形2:iijy。情形2.1,当zyjj时:情形2.2,当i-1i-1yzj时:最后,我们证明决策si,j满足单调性。为讨论方便,令mki,j=mi-1,k+wk+1,j;我们先来证明si-1,jsi,j,只要证明对于所有ikkj且m
18、ki-1,jmki-1,j,有:mki,jmki,j。类似地,我们可以证明一个更强的不等式mki-1,j-mki-1,jmki,j-mki,j也就是: mki-1,j+mki,jmki,j+mki-1,j利用递推定义式展开整理的:mi-2,k+mi-1,kmi-1,k+mi-2,k,这就是i-2i-1kk时m的四边形不等式。我们再来证明si,jsi,j+1,与上文类似,设kkj,则我们只需证明一个更强的不等式: mki,j-mki,jmki,j+1-mki,j+1也就是: mki,j+mki,j+1mki,j+1+mki,j利用递推定义式展开整理的:wk+1,j+wk+1,j+1wk+1,j+
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 动态规划算法的优化技巧 算法 动态 规划 优化 技巧
限制150内