最长公共子序列——动态规划

定义：设X和Z是两个序列，其中X=<x₁,x₂,…,x_m>，Z=<z₁,z₂,…,z_k>，如果存在X的元素构成的按下标严格排序递增序列<x_i1,x_i2,…,x_ik>，使得x_ij=z_j，j=1,2,…,k，那么Z是X的子序列，Z含有的元素个数，称为子序列的长度。

定义：设X和Y是两个序列，如果Z既是X的子序列，也是Y的子序列，则称Z是X和Y的公共子序列。

问题：最长公共子序列问题，给定序列X=<x₁,x₂,…,x_m>，Y=<y₁,y₂,…,y_j>，就X和Y的最长公共子序列。

解析

Xi=<x1,x2,…,xi>
Yj=<y1,y2,…,yj>
Zk=<z1,z2,…,zk>
如果Zk是Xi和Yj的最长公共子序列，
（1）xi = yj，那么zk = xi = yj，Zk-1是Xi-1和Yj-1的最长公共子序列
（2）xi ≠ yj，那么zk ≠ xi，Zk-1是Xi-1和Yj的最长公共子序列
（3）xi ≠ yj，那么zk ≠ yi，Zk-1是Xi和Yj-1的最长公共子序列
实例：
X=<A,B,C,D,C>
Y=<B,D,C,B>
m=1-5
n=1-4
初始化：

（1）i=1
①j=1 X.A<>Y.B,C[1,1]=max{C[1,0],C[0,1]}=0,删除y
②j=2 X.A<>Y.D,C[1,2]=max{C[1,1],C[0,2]}=0,删除y
③j=3 X.A<>Y.C,C[1,3]=max{C[1,2],C[0,3]}=0,删除y
④j=4 X.A<>Y.B,C[1,4]=max{C[1,3],C[0,4]}=0,删除y
（2）i=2
①j=1 X.B等于Y.B,C[2,1]=C[1,1]+1=1,删除两个
②j=2 X.B<>Y.D,C[2,2]=max{C[2,1],C[1,2]}=1,删除y
③j=3 X.B<>Y.C,C[2,3]=max{C[2,2],C[1,3]}=1,删除y
④j=4 X.B=Y.B,C[2,4]=C[1,3]+1=1,删除两个
（3）i=3
①j=1 X.C<>Y.B,C[3,1]=max{C[3,0],C[2,1]}=1,删除x
②j=2 X.C<>Y.D,C[3,2]=max{C[3,1],C[2,2]}=1,删除y
③j=3 X.C等于Y.C,C[3,3]=C[2,2]+1=2,删除两个
④j=4 X.C<>Y.B,C[3,4]=max{C[3,3],C[2,4]}=2,删除y
（4）i=4
①j=1 X.D<>Y.B,C[4,1]=max{C[4,0],C[3,1]}=1,删除x
②j=2 X.D=Y.D,C[4,2]=C[3,1]+1=2,删除两个
③j=3 X.D<>Y.C,C[4,3]=max{C[4,2],C[3,3]}=2,删除y
④j=4 X.D<>Y.B,C[4,4]=max{C[4,3],C[3,4]}=2,删除y
（5）i=5
①j=1 X.C<>Y.B,C[5,1]=max{C[5,0],C[4,1]}=1,删除x
②j=2 X.C<>Y.D,C[5,2]=max{C[5,1],C[4,2]}=1,删除x
③j=3 X.C==Y.C,C[5,3]=C[4,2]+1=3,删除两个
④j=4 X.C<>Y.B,C[5,4]=max{C[5,3],C[4,4]}=3,删除y

（1）X=5，Y=4
查表(5,4) “3；删除y”
X=<A,B,C,D,C>
Y=<B,D,C,B>
=>
X=<A,B,C,D,C>
Y=<B,D,C>
（2）X=5，Y=3
查表(5,3) “3；删除两个”
X=<A,B,C,D,C>
Y=<B,D,C>
=>
X=<A,B,C,D>
Y=<B,D>
输出C
（3）X=4，Y=2
查表(4,2) “2；删除两个”
X=<A,B,C,D>
Y=<B,D>
=>
X=<A,B,C>
Y=< B >
输出D
（4）X=3，Y=1
查表(3,1) “1；x”
X=<A,B,C>
Y=< B >
=>
X=<A,B>
Y=< B >
（5）X=2，Y=1
查表(2,1) “1；删除两个”
X=<A,B>
Y=< B >
=>
X=< A >
Y=<>
输出B
（6）X=1，Y=0
最终输出<B,D,C>

算法

算法1：给出最长子串长度
C[0,j]=C[i,0]=0,1≤i≤m,1≤j≤n

• For i=1 to m
  • For j=1 to n
    • C[i,j]=C[i-1,j-1]+1,(i,j>0，xi = yj，B[i,j]=←（删除两个）)；或max{C[i,j-1],C[i-1,j]},(i,j>0，xi ≠ yj，C[i,j-1]，B[i,j]=↓（删除y）；C[i-1,j]，B[i,j]=↑（删除x）)

算法2：f（B,i,j）输出最长子串
If i=0 or j=0
    then return;
If B[i,j]=← //删除两个
    then 输出X[i]
Else if B[i,j]=↑ //删除x
    f(B,i-1,j)
Else f(B,i,j-1) //删除y

#include <iostream>
#include <string>
#include <vector>
using namespace std;
//算法2：f（B, i, j）输出最长子串
void OutPutLCS(vector<vector<int>> &B, string X, int X_len, int Y_len) {
	if (X_len == 0 || Y_len == 0) {
		return;
	}
	if (B[X_len][Y_len] == 3) {
		OutPutLCS(B, X, X_len - 1, Y_len - 1);
		cout << X[X_len - 1];
	}
	else if (B[X_len][Y_len] == 2) {
		OutPutLCS(B, X, X_len - 1, Y_len);
	}
	else if (B[X_len][Y_len] == 1) {
		OutPutLCS(B, X, X_len, Y_len - 1);
	}
}
int main()
{
	string X = "ABCDC";
	string Y = "BDCB";
	int m = X.size();
	int n = Y.size();
	vector<vector<int>> C(m + 1, vector<int>(n + 1));
	vector<vector<int>> B(m + 1, vector<int>(n + 1));//1：删除y；2：删除x；3：删除两个
	//算法1：给出最长子串长度
	//初始化
	for (int i = 0; i <= m; i++) {
		C[i][0] = 0;
	}
	for (int i = 0; i <= n; i++) {
		C[0][i] = 0;
	}
	for (int i = 1; i <= m; i++) {
		for (int j = 1; j <= n; j++) {
			if (X[i - 1] == Y[j - 1]) {
				C[i][j] = C[i - 1][j - 1] + 1;
				B[i][j] = 3;//删除两个
			}
			else {
				if (C[i][j - 1] >= C[i - 1][j]) {
					C[i][j] = C[i][j - 1];
					B[i][j] = 1;//删除y
				}
				else {
					C[i][j] = C[i - 1][j];
					B[i][j] = 2;//删除x
				}
			}
		}
	}
	for (int i = 1; i <= m; i++) {
		for (int j = 1; j <= n; j++) {
			cout << C[i][j] << " ";
		}
		cout << endl;
	}
	cout << "最长公共子序列长度为：" << C[m][n] << endl;
	OutPutLCS(B, X, m, n);
	cout << endl;
}

计算结果

分析

复杂度 O(mn)

github源码地址