程序设计实习(II)19_动态规划教学教案.pptVIP

最长公共子序列 Sample Input abcfbc abfcab programming contest abcd mnp Sample Output 4 2 0 输入两个子串s1,s2, 设MaxLen(i,j)表示: s1的左边i个字符形成的子串，与s2左边的j个字符形成的子串的最长公共子序列的长度 MaxLen(i,j) 就是本题的“状态” 假定 len1 = strlen(s1),len2 = strlen(s2） 那么题目就是要求 MaxLen(len1,len2) 最长公共子序列 显然：MaxLen(n,0) = 0 ( n= 0…len1） MaxLen(0,n) = 0 ( n=0…len2） 递推公式： if ( s1[i-1] == s2[j-1] ) MaxLen(i,j) = MaxLen(i-1,j-1) + 1; else MaxLen(i,j) = Max(MaxLen(i,j-1), MaxLen(i-1,j) ); 最长公共子序列 S1 s1[i-1] S2 s2[j-1] S1i-1 S2j-1 S1长度为 i S2长度为 j MaxLen(S1,S2)不会比MaxLen(S1,S2j-1) 和 MaxLen(S1i-1,S2)都大，也不会比两者中任何一个小 #include iostream.h #include string.h char sz1[1000]; char sz2[1000]; int anMaxLen[1000][1000]; main() { while( cin sz1 sz2 ) { int nLength1 = strlen( sz1); int nLength2 = strlen( sz2); int nTmp; int i,j; for( i = 0;i = nLength1; i ++ ) anMaxLen[i][0] = 0; for( j = 0;j = nLength2; j ++ ) anMaxLen[0][j] = 0; for( i = 1;i = nLength1;i ++ ) { for( j = 1; j = nLength2; j ++ ) { if( sz1[i-1] == sz2[j-1] ) anMaxLen[i][j] = anMaxLen[i-1][j-1] + 1; else { int nLen1 = anMaxLen[i][j-1]; int nLen2 = anMaxLen[i-1][j]; if( nLen1 nLen2 ) anMaxLen[i][j] = nLen1; else anMaxLen[i][j] = nLen2; } } } cout anMaxLen[nLength1][nLength2] endl; } } 程序设计实习 动态规划 树形递归存在冗余计算 例1：POJ 2753 Fibonacci数列 求 Fibonacci数列的第n项 int f(int n){ if(n==0 || n==1) return n; return f(n-1)+f(n-2); } f(5) f(3) f(2) f(1) f(2) f(4) f(0) f(1) f(0) f(3) f(2) f(1) f(1) f(0) f(1) 1 1 0 0 1 0 1 0 冗余计算 计算过程中存在冗余计算，为了除去冗余计算 可以从已知条件开始计算，并记录计算过程中 的中间结果。 树形递归存在冗余计算 去除冗余： int f[n+1]; f[1]=f[2]=1; int I; for(i=3;i=n;i++) f[i] = f[i-1]+f[i-2]; cout f[n] endl; 用空间换时间 - 动态规划 例2　POJ1163 数字三角形 在上面的数字三角形中寻找一条从顶部到底边的路径，使得路径上所经过的数字之和最大。路径上的每一步都只能往左下或右下走。只需要求出这个最大和即可，不必给出具体路径。 三角形的行数大于1小于等于100，数字为 0 - 99 7 3 8 8 1 0 2 7 4 4 4 5 2 6 5 #include iostream