经典算法分析与设计教学案例.ppt
《经典算法分析与设计教学案例.ppt》由会员分享,可在线阅读,更多相关《经典算法分析与设计教学案例.ppt(57页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第3章 动态规划 Those who cannot remember the past are doomed to repeat itGeorge Santayana 1 学习要点学习要点:理解动态规划算法的概念。掌握动态规划算法的基本要素(1)最优子结构性质(2)重叠子问题性质掌握设计动态规划算法的步骤。(1)找出最优解的性质,并刻划其结构特征。(2)递归地定义最优值。(3)以自底向上的方式计算出最优值。(4)根据计算最优值时得到的信息,构造最优解。2通过应用范例学习动态规划算法设计策略。(1)矩阵连乘问题;(2)最长公共子序列;(3)最大子段和(4)凸多边形最优三角剖分;(5)多边形游戏;
2、(6)图像压缩;(7)电路布线;(8)流水作业调度;(9)背包问题;(10)最优二叉搜索树。3n动态规划算法与分治法类似,其基本思想也是将待求解问题分解成若干个子问题算法总体思想nT(n/2)T(n/2)T(n/2)T(n/2)T(n/2)T(n/2)T(n/2)T(n/2)T(n)=4n但是经分解得到的子问题往往不是互相独立的。不同子问题的数目常常只有多项式量级。在用分治法求解时,有些子问题被重复计算了许多次。算法总体思想nT(n)=n/2T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)n/2T(n/4)T(n/4)T(n/4)T(n/4)T(
3、n/4)T(n/4)T(n/4)T(n/4)n/2T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)n/2T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)5n如果能够保存已解决的子问题的答案,而在需要时再找出已求得的答案,就可以避免大量重复计算,从而得到多项式时间算法。算法总体思想n=n/2T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)n/2n/2T(n/4)T(n/4)T(n/4)T(n/4)n/2T(n/4)T(n/4)T(n/4)T(n/4)T(
4、n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n/4)T(n)6求组合数求组合数每次求解可将其分为两个子问题的和。把 按递推式分解得到。7首先,将的位置上以及0的位置上元素皆填为1。填某一中间表目时,只要把它右边表目的元素与右下方表目的元素之和填入即可。这样,很快就能求出 10。8动态规划解决问题的基本特征1.动态规划一般解决最值(最优,最大,最小,最长)问题;2.动态规划解决的问题一般是离散的,可以分解(划分阶段)的;3.动态规划解决的问题必须包含最优子结构,即可以由(n1)的最优推导出n的最优9动态规划基本步骤n找出最优解的性质,并刻划其结构特征。n递归地定义最优值。n以自底
5、向上的方式计算出最优值。n根据计算最优值时得到的信息,构造最优解。10矩阵连乘问题n给定n个矩阵 ,其中 与 是可乘的,。考察这n个矩阵的连乘积 n由于矩阵乘法满足结合律,所以计算矩阵的连乘可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。n若一个矩阵连乘积的计算次序完全确定,也就是说该连乘积已完全加括号,则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积11(1)单个矩阵是完全加括号的;(2)矩阵连乘积 是完全加括号的,则 可 表示为2个完全加括号的矩阵连乘积 和 的乘积并加括号,即 16000,10500,36000,87500,34500u完全加括号的矩阵连乘积可
6、递归地定义为:u设有四个矩阵 ,它们的维数分别是:u总共有五中完全加括号的方式完全加括号的矩阵连乘积12矩阵连乘问题给定n个矩阵A1,A2,An,其中Ai与Ai+1是可乘的,i=1,2,n-1。如何确定计算矩阵连乘积的计算次序,使得依此次序计算矩阵连乘积需要的数乘次数最少。u穷举法穷举法:列举出所有可能的计算次序,并计算出每一种计算次序相应需要的数乘次数,从中找出一种数乘次数最少的计算次序。算法复杂度分析:算法复杂度分析:对于n个矩阵的连乘积,设其不同的计算次序为P(n)。由于每种加括号方式都可以分解为两个子矩阵的加括号问题:(A1.Ak)(Ak+1An)可以得到关于P(n)的递推式如下:13
7、矩阵连乘问题u穷举法穷举法u动态规划动态规划将矩阵连乘积 简记为Ai:j,这里ij 考察计算Ai:j的最优计算次序。设这个计算次序在矩阵Ak和Ak+1之间将矩阵链断开,ikj,则其相应完全加括号方式为计算量:Ai:k的计算量加上Ak+1:j的计算量,再加上Ai:k和Ak+1:j相乘的计算量14n特征:计算Ai:j的最优次序所包含的计算矩阵子链 Ai:k和Ak+1:j的次序也是最优的。n矩阵连乘计算次序问题的最优解包含着其子问题的最优解。这种性质称为最优子结构性质最优子结构性质。问题的最优子结构性质是该问题可用动态规划算法求解的显著特征。分析最优解的结构15建立递归关系n设计算Ai:j,1ijn
8、,所需要的最少数乘次数mi,j,则原问题的最优值为m1,n n当i=j时,Ai:j=Ai,因此,mi,i=0,i=1,2,nn当ij时,n可以递归地定义mi,j为:这里 的维数为 的位置只有 种可能16计算最优值n对于1ijn不同的有序对(i,j)对应于不同的子问题。因此,不同子问题的个数最多只有n由此可见,在递归计算时,许多子问题被重复计算多次许多子问题被重复计算多次。这也是该问题可用动态规划算法求解的又一显著特征。n用动态规划算法解此问题,可依据其递归式以自底向上的方式进行计算。在计算过程中,保存已解决的子问题答案。每个子问题只计算一次,而在后面需要时只要简单查一下,从而避免大量的重复计算
9、,最终得到多项式时间的算法。17用动态规划法求最优解void MatrixChain(int*p,int n,int*m,int*s)for(int i=1;i=n;i+)mii=0;for(int r=2;r=n;r+)for(int i=1;i=n-r+1;i+)int j=i+r-1;mij=mi+1j+pi-1*pi*pj;sij=i;for(int k=i+1;k j;k+)int t=mik+mk+1j+pi-1*pk*pj;if(t 0)return mij;if(i=j)return 0;int u=LookupChain(i,i)+LookupChain(i+1,j)+pi-
10、1*pi*pj;sij=i;for(int k=i+1;k j;k+)int t=LookupChain(i,k)+LookupChain(k+1,j)+pi-1*pk*pj;if(t u)u=t;sij=k;mij=u;return u;21最长公共子序列n一个给定序列的一个给定序列的子序列子序列是在该序列中删去若干元素后得到的是在该序列中删去若干元素后得到的序列。序列。n给定两个序列给定两个序列X和和Y,当另一序列当另一序列Z既是既是X的子序列又是的子序列又是Y的子的子序列时,称序列时,称Z是序列是序列X和和Y的的公共子序列。公共子序列。n最长公共子序列最长公共子序列:公共子序列中长度最长
11、的子序列。公共子序列中长度最长的子序列。n最长公共子序列问题最长公共子序列问题 给定两个序列给定两个序列X=x1,x2,xm和和Y=y1,y2,yn,找出找出X和和Y的一个最长公共子序列。的一个最长公共子序列。n例如例如:X=(A,B,C,B,D,A,B)X的子序列的子序列:q所有所有X的子集的子集(集合中元素按原来在集合中元素按原来在X中的顺序排列中的顺序排列)(A,B,D),(B,C,D,B),等等等等.22最长公共子序列X=(A,B,C,B,D,A,B)X=(A,B,C,B,D,A,B)Y=(B,D,C,A,B,A)Y=(B,D,C,A,B,A)n(B,C,B,A)和(B,D,A,B)都
12、是X和Y 的最长公共子序列(长度为4)n但是,(B,C,A)就不是X和Y的最长公共子序列例:例:23最长公共子序列n对于每一个Xm的子序列,验证它是否是Yn的子序列.nXm有2m个子序列n每个子序列需要o(n)的时间来验证它是否是Yn的子序列.q从Yn的第一个字母开始扫描下去,如果不是则从第二个开始n运行时间:o(n2m)穷举法:24最长公共子序列的结构设序列X=x1,x2,xm和Y=y1,y2,yn的最长公共子序列为Z=z1,z2,zk,则:(1)若xm=yn,则zk=xm=yn,且zk-1是xm-1和yn-1的最长公共子序列。(2)若xmyn且zkxm,则Z是xm-1和Y的最长公共子序列。
13、(3)若xmyn且zkyn,则Z是X和yn-1的最长公共子序列。由此可见,2个序列的最长公共子序列包含了这2个序列的前缀的最长公共子序列。因此,最长公共子序列问题具有最优子结最优子结构性质构性质。25子问题的递归结构由最长公共子序列问题的最优子结构性质建立子问题最优值的递归关系。用cij记录序列和的最长公共子序列的长度。其中,Xi=x1,x2,xi;Yj=y1,y2,yj。当i=0或j=0时,空序列是Xi和Yj的最长公共子序列。故此时Cij=0。其它情况下,由最优子结构性质可建立递归关系如下:26计算最优值 由于在所考虑的子问题空间中,总共有(mn)个不同的子问题,因此,用动态规划算法自底向上
14、地计算最优值能提高算法的效率。矩阵 bi,j 告诉我们要获得最优结果应该如何选择:l 如果 xi=yj bi,j=1l 如果 ci-1,j ci,j-1 bi,j=2 否则 bi,j=3void LCSLength(int m,int n,char*x,char*y,int*c,int*b)int i,j;for(i=1;i=m;i+)ci0=0;for(i=0;i=n;i+)c0i=0;for(i=1;i=m;i+)for(j=1;j=cij-1)cij=ci-1j;bij=2;else cij=cij-1;bij=3;00000000000yj:DACFABxij012nm12033CDc
15、i,j-1ci-1,ji27计算最优值X=(A,B,C,B,D,A,B)Y=(B,D,C,A,B,A)0126345yjBDACAB51203467DABxiCBAB00000000000000000 11 1 1111 2211 2222 1122 331 222331232 3 4 1223 44如果 xi=yj bi,j=“”如果 ci-1,jci,j-1bi,j=“”否则bi,j=“”例:28构造最优解构造最长公共子序列构造最长公共子序列void LCS(int i,int j,char*x,int*b)if(i=0|j=0)return;if(bij=1)LCS(i-1,j-1,x,
16、b);coutxi;else if(bij=2)LCS(i-1,j,x,b);else LCS(i,j-1,x,b);29最大子段和n问题描述:给定由n个整数(包含负整数)组成的序列a1,a2,.,an,求该序列子段和的最大值。当所有整数均为负值时定义其最大子段和为0。q依此定义,所求的最优值为:例如,当(a1,a2,a3,a4,a5,a6)=(-2,11,-4,13,-5,-2)时,最大子段和为:30最大子段和1.一个简单算法一个简单算法int MaxSum(int n,a,&besti,&bestj)int sum=0;for(i=1;i=n;i+)for(j=1;j=n;j+)int t
17、hissum=0;for(k=i;ksum)sum=thissum;besti=i;Bestj=j;return sum;n算法有3重循环,复杂性为O(n3)。n由于有:n算法可作如下改进:int MaxSum(int n,a,&besti,&bestj)int sum=0;for(i=1;i=n;i+)int thissum=0;for(j=i;jsum)sum=thissum;besti=i;bestj=j;n改进后的算法复杂性为O(n2)。31最大子段和2.分治方法求解分治方法求解n从问题的解的结构可以看出,它适合于用分治策略求解:q如果将所给的序列a1:n分为长度相等的两段a1:n/2
18、和an/2+1:n,分别求出这两段的最大子段和,则a1:n的最大子段和有三种情形:A.a1:n的最大子段和与的最大子段和与a1:n/2的最大子段和相同;的最大子段和相同;B.a1:n的最大子段和与的最大子段和与an/2+1:n的最大子段和相同;的最大子段和相同;C.a1:n的最大子段和为下面的形式。的最大子段和为下面的形式。qA、B这两种情形可递归求得。对于情形C,容易看出,an/2与an/2+1在最优子序列中。因此,我们可以在a1:n/2和an/2+1:n中分别计算出如下的s1和s2。则s1+s2即为出现情形C使得最优值。q int MaxSubSum(int a,int left,int
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 经典 算法 分析 设计 教学 案例
限制150内