1 最长公共子序列.docxVIP

1 最长公共子序列/blog/日积月累-算法与 2009-10-19 11:08:33 阅读577 评论0 字号：大中小 订阅 一个给定序列的子序列是在该序列中删去若干元素后得到的序列。最长公共子序列问题是给定两个序列X=｛x1,x2...xm}和Y={y1,y2...yn}，找出X和Y的最长公共子序列，用动态规划算法可有效的解些问题。最长公共子序列的结构：设X = { x1 , ... , xm },Y = { y1 , ... , yn }及它们的最长子序列Z = { z1 , ... , zk }则1、若 xm = yn ， 则 zk = xm = yn，且Z[k-1] 是 X[m-1] 和 Y[n-1] 的最长公共子序列2、若 xm != yn ，且 zk != xm , 则 Z 是 X[m-1] 和 Y 的最长公共子序列3、若 xm != yn , 且 zk != yn , 则 Z 是 Y[n-1] 和 X 的最长公共子序列子问题的递归结构：当 i = 0 , j = 0 时 , c[i][j] = 0当 i , j 0 ; xi = yi 时 , c[i][j] = c[i-1][j-1] + 1当 i , j 0 ; xi != yi 时 , c[i][j] = max { c[i][j-1] , c[i-1][j] }由以上分析，先计算最优值，再构造最长公共子序列，C++源程序如下：#include stdio.h #include string.h#includeiostreamusing namespace std;int length=0; //记录最长公共序列长度 //计算最优值 void lcs(string x,string y,int **b){int i,j;//获取字符串的长度 int m = x.size();int n = y.size();//创建动太二维数组 int **c;c = new int*[m+1];for(int i = 0;i m+1;i++)c[i] = new int[n+1];//初始化数组 for (i=1;i=m;i++) c[i][0] = 0;for (i=1;i=n;i++) c[0][i] = 0;c[0][0] = 0;//遍历记录最优值，并将标识存在数组b中for (i=1;i=m;i++)for (j=1;j=n;j++){if (x[i-1] == y[j-1]){ //xm = yn ， 则 zk = xm = yn，且Z[k-1] 是 X[m-1] 和 Y[n-1] 的最长公共子序列c[i][j] = c[i-1][j-1] + 1;b[i][j] = 1;}else if (c[i-1][j] c[i][j-1]){ //xm != yn ，且 zk != xm , 则 Z 是 X[m-1] 和 Y 的最长公共子序列c[i][j] = c[i-1][j];b[i][j] = 2;}else{ //xm != yn , 且 zk != yn , 则 Z 是 Y[n-1] 和 X 的最长公共子序列c[i][j] = c[i][j-1];b[i][j] = 3;}}//释放空间 for(int i = 0;i m+1;i++)delete c[i];delete c;}//构造最长公共子序列 void show(int i,int j,string x,int **b){//递归中止条件 if (i==0||j==0)return ;//从b[m][n]开始，用递归依其值在数组b中搜索if (b[i][j]==1){ //xi和yj的最长公共子序列是由x(i-1)和y(j-1)的最长公共子序列在尾部加上xi所得到的子序列show(i-1,j-1,x,b);length++;cout x[i-1];}else if (b[i][j]==2) { //xi和yi的最长公共子序列与x(i-1)和yj的最长公共子序列相同show(i-1,j,x,b);} else { //xi和yi的最长公共子序列与xi和y(j-1)的最长公共子序列相同show(i,j-1,x,b); }} int main() { string x;string y;int **b;cout 输入第一个序列; cin x;cout 输入第二个序列; cin y;int m = x.size();int n = y.size();//创建二维动态数组 b = new int*[m+1];for(int i = 0;i m+1;i++)b[i] = new int[n+1];//求解最优值 lcs(x,y,b);//构造序列并打印 cout 最长公共序列:;sh