最长公共子序列实现

一个给定序列的子序列是在该序列中删去若干元素后得到的序列。给定两个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是序列X和Y的公共子序列。最长公共子序列就是求给定两个序列的一个最长公共子序列。动态规划可以有效的解决此问题。

给定两个序列
X = { x1 , x2 , ... , xm }
Y = { y1 , y2 , ... , yn }
求X和Y的一个最长公共子序列

举例
X = { a , b , c , b , d , a , b }
Y = { b , d , c , a , b , a }
最长公共子序列为
LSC = { b , c , b , a }

分析：

最长公共子序列问题具有最优子结构性质

设

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 的最长公共子序列

由最长公共子序列问题的子序列的最优子结构性质，可以建立子问题最优的递归关系。用c[i][j]记录序列Xi和Yi的最长公共子序列的长度。

当 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] }

 

 

1. /**************************************************************
2.  * File: lcs.cpp
3.  * Purpose: get the longest common subsequence of two sequences
4.  * Author: cxf
5.  * Date: 07/08/2008
6.  ***************************************************************/
9. #include <iostream>
10. using namespace std;
12. /**************************************
13.  根据b记录的子问题的类型,求最长子序列                                            
14. /**************************************/
15. void lcs(int m, int n, char *x, int **b)
16. {
17.     if (0 == m || 0 == n)
18.     {
19.         return;
20.     }
21.     if ('//' == b[m][n])
22.     {
23.         lcs(m - 1, n - 1, x, b);
24.         cout << x[m] << endl;
25.     }
26.     else if ('|' == b[m][n])
27.     {
28.         lcs(m - 1, n, x, b);
29.     }
30.     else if ('-' == b[m][n])
31.     {
32.         lcs(m, n - 1, x, b);
33.     }
34. }
36. /**********************************************
37.  求两个序列的最长子序列,x和y的元素需从下标1开始
38.  m--x的长度
39.  n--y的长度
40. ***********************************************/
41. void lcsMain(int m, int n, char *x, char *y, int **c, int **b)
42. {
43.     int i, j;
44.     /* 初始化长度数组 */
45.     for (i = 0; i < m; i++)
46.     {
47.         for (j = 0; j < n; j++)
48.         {
49.             c[i][j] = 0;
50.             b[i][j] = 0;
51.         }
52.     }
53.     /* 求X和Y各个子序列的最大长度，并记录各个长度是由那个子问题得到的,由b保存
54.      * x和y中的元素需从下标1开始
55.      */
56.     for (i = 1; i < m; i++)
57.     {
58.         for (j = 1; j < n; j++)
59.         {
60.             if (x[i] == y[j])
61.             {
62.                 c[i][j] = c[i - 1][j - 1] + 1;
63.                 b[i][j] = '//';
64.             }
65.             else if (c[i - 1][j] > c[i][j - 1])
66.             {
67.                 c[i][j] = c[i - 1][j];
68.                 b[i][j] = '|'; 
69.             }
70.             else
71.             {
72.                 c[i][j] = c[i][j - 1];
73.                 b[i][j] = '-'; 
74.             }
75.         }
76.     }   
77.     cout << "最长子序列长度：" << c[i - 1][j - 1] << endl;//打印最长子序列长度
78.     lcs(m - 1, n - 1, x, b);    //求最长子序列
79. }
82. int main()
83. {
84.     char *x = " bcdefac";
85.     char *y = " abafdfa";
86.     int xLen = strlen(x);
87.     int yLen = strlen(y);
89.     /* 申请动态二维数组 */
90.     int **c = new int*[xLen];
91.     int **b = new int*[xLen];
92.     for (int i = 0; i < xLen; i++)
93.     {
94.         c[i] = new int[yLen];
95.         b[i] = new int[yLen];
96.     }
98.     lcsMain(xLen, yLen, x, y, c, b);
100.     /* 释放内存 */
101.     for (i = 0; i < xLen; i++)
102.     {
103.         delete[] c[i];
104.         delete[] b[i];
105.     }
106.     delete[] c;
107.     delete[] b;
108.     return 0;
109. }