算法设计与分析第3章.ppt

第3章 动态规划 学习要点: 理解动态规划算法的概念。 掌握动态规划算法的基本要素 （1）最优子结构性质 （2）重叠子问题性质 掌握设计动态规划算法的步骤。 (1)找出最优解的性质，并刻划其结构特征。 (2)递归地定义最优值。 (3)以自底向上的方式计算出最优值。 (4)根据计算最优值时得到的信息，构造最优解。 通过应用范例学习动态规划算法设计策略。 （1）矩阵连乘问题； （2）最长公共子序列； （3）最大子段和 （4）流水作业调度 （5）背包问题； （6）最优二叉搜索树。 动态规划基本步骤 找出最优解的性质，并刻划其结构特征。 递归地定义最优值。 以自底向上的方式计算出最优值。 根据计算最优值时得到的信息，构造最优解。 矩阵连乘问题 二个矩阵乘积的标准算法 void matrixMultiply(int **a,int **b,int **c,int ra,int ca,int rb,int cb) {if(ca!=rb) error(“矩阵不可乘”); for(i=0;ira;i++) for(int j=0;jrb;j++){ int sum=a[i][0]*b[0][j]; for(int k=1;kca;k++); sum+=a[i][k]*b[k][j]; c[i][j]=sum; } } 按照如下次序计算 12 23 34 45 56 13 24 35 46 14 25 36 15 26 16 下面行的计算结果将会直接用到上面的结果。 main() { int n=6;//矩阵的个数 int *p=new int[n+1]; //p[0]:第一个矩阵的行数 //p[1]:第一个矩阵的列数，第二个矩阵的行数 //p[2]:第二个矩阵的列数，第三个矩阵的行数 p[0]=30; p[1]=35; p[2]=15; p[3]=5; p[4]=10; p[5]=20; p[6]=25; int **m,**s; m=new int*[n+1]; for( int i=0;in+1;i++) m[i]=new int[n+1]; s=new int*[n+1]; for(i=0;in+1;i++) s[i]=new int[n+1]; MatrixChain(p,n,m,s); Traceback(1,n,s); } void MatrixChain(int *p,int n,int **m,int **s) { for(int i=1;i=n;i++) m[i][i]=0;//单个矩阵相乘，所需数乘次数为0 for(int r=2;r=n;r++) //对角线循环 { for(int i=1;i=n-r+1;i++) //行循环 { int j=i+r-1; //列的控制 //首先在i断开，即（Ai*（Ai+1...Aj）） m[i][j]=m[i][i]+m[i+1][j]+p[i-1]*p[i]*p[j]; s[i][j]=i; for(int k=i+1;kj;k++) { //然后在k（从i+1开始遍历到j-1）断开，即（（Ai...Ak)*（Ak+1...Aj）） int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j]; if(tm[i][j])//找到更好的断开方法 { m[i][j]=t;//记录最少数乘次数 s[i][j]=k;//记录断开位置 }}}}} void Traceback(int i,int j,int **s) { //s[i][j]记录了断开的位置，即计算A[i:j]的加括号方式为： //(A[i:s[i][j]])*(A[s[i][j]+1:j]) if(i==j)return; Traceback(i,s[i][j],s);//递归打印A[i:s