合工大程序设计艺术与方法实验四动态规划.doc

. . . .......... 《程序设计艺术与方法》课程实验报告 实验名称 实验四 动态规划 姓 名 班 级 学 号 实验日期 2018.6.12 指导教师 徐本柱 成 绩 一、实验目的和要求 (1) 理解动态规划的基本思想、动态规划算法的基本步骤。 (2) 掌握动态规划算法实际步骤。 二、实验预习内容 动态规划 三、实验项目摘要 (1) 求两个字符串的最长公共子序列。 X 的一个子序列是相应于X 下标序列{1, 2, …, m}的一个子序列，求解两个序列的所有 子序列中长度最大的，例如输入：pear, peach 输出：pea。 (2) 给定两个字符串a 和b，现将串a 通过变换变为串b，可用的操作为，删除串a 中的一 个字符；在串a 的某个位置插入一个元素；将串a 中的某个字母换为另一个字母。对于 任意的串a 和串b，输出最少多少次能够将串变为串b。 思考：输出变换的步骤。 (3) 输入一个矩阵，计算所有的子矩阵中和的最大值。 例如，输入 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2 输出为：15 四、实验结果与分析（源程序及相关说明） 求两个字符串的最长公共子序列 算法思想：通过动态规划求解：seat二维数组用于表示选择的最长公共子序列以及输出时选择的方向，其中0表示可选的子序列点，1表示递归选择str1（0，i-1），str2（0，j)的公共子序列,-1选择str1（0，i），str2（0，j-1)的公共子序列，Long表示选择str1（0，i），str2（0，j)的公共子序列的长度，其中当递归结束时，Long[m][n]就存储着最大的长度。 先将Long数组第一行第一列初始化为0，方便计算，接着如果str[i] == str[j],则seat[i][j]为可选点，seat[i][j]=0，同时Long[i][j] = Long[i - 1][j - 1] + 1，如果str[i] ！= str[j]，则Long[i][j]=max（Long[i - 1][j] Long[i][j - 1]），如果选择前者seat[i][j]=1，否则seat[i][j]=-1；当循环遍历了两个字符串后，就得出结论，在根据seat中存储的元素值，从m，n开始（m,n分别为两字符串的长度），先递归，后输出对应位置在字符串中的字符，递归结束，就可以输出字符串。 #include stdafx.h #includestdio.h #includestring #includeiostream using namespace std; #define MAXLEN 100 int seat[MAXLEN][MAXLEN];//其中0表示可选的子序列点，1表示选择str1（0，i-1），str2（0，j)的公共子序列,-1选择str1（0，i），str2（0，j-1)的公共子序列 int Long[MAXLEN][MAXLEN];//表示选择str1（0，i），str2（0，j)的公共子序列的长度 void LCSLength(string str1, string str2, int m, int n) { int i, j; for (i = 0; i = m; i++) Long[i][0] = 0; for (j = 1; j = n; j++) Long[0][j] = 0; //第一行第一列不使用，方便计算 for (i = 1; i = m; i++) { for (j = 1; j = n; j++) { if (str1[i - 1] == str2[j - 1]) { Long[i][j] = Long[i - 1][j - 1] + 1; seat[i][j] = 0; } else if (Long[i - 1][j] = Long[i][j - 1]) { Long[i][j] = Long[i - 1][j]; seat[i][j] = 1; } else { Long[i][j] = Long[i][j - 1]; seat[i][j] = -1; } } } } //以str1为标准输出 bool Print(string str1, int i, int j,int m,int n) { if (i == 0 || j == 0) return true; if (seat[i][j] == 0) {