若给定序列X={x1,x2,…,xm},则另一序列Z={z1,z2,…,zk},是X的子序列是指存在一个严格递增下标序列{i1,i2,…,ik}使得对于所有j=1,2,…,k有:zj=xij。
例如,序列Z={B,C,D,B}是序列X={A,B,C,B,D,A,B}的子序列,相应的递增下标序列为{2,3,5,7}。
给定2个序列X和Y,当另一序列Z既是X的子序列又是Y的子序列时,称Z是序列X和Y的公共子序列。
给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn},找出X和Y的最长公共子序列。
设序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}的最长公共子序列为Z={z1,z2,…,zk} ,则
1.若xm=yn,则集合Z的特点:
zk=xm=yn,且Zk-1是Xm-1和Yn-1的最长公共子序列。
2.若xm≠yn且,则Z的特点:
是Xm-1和Y的最长公共子序列。
X和Yn-1的最长公共子序列。
最优子结构的证明:
设序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}的最长公共子序列为Z={z1,z2,…,zk}
1.若xm=yn, 则zk=xm=yn,且Zk-1是Xm-1和Yn-1的最长公共子序列。
2.若xm≠yn且zk≠xm, 则Z是Xm-1和Y的最长公共子序列。
3.若xm≠yn且zk≠yn, 则Z是X和Yn-1的最长公共子序列。
计算最优值
#define NUM 100
int c[NUM][NUM]; //x[1-i] 和 y[1-j]之间的公共子序列的长度。
int b[NUM][NUM]; // 辅助完成最优解得计算
void LCSLength (int m, int n, const char x[],char y[]){
int i,j;
//数组c的第0行、第0列置0
for (i = 1; i <= m; i++) c[i][0] = 0;
for (i = 1; i <= n; i++) c[0][i] = 0;
//根据递推公式构造数组c
for (i = 1; i <= m; i++)
for (j = 1; j <= n; j++){
if (x[i]==y[j])
{c[i][j]=c[i-1][j-1]+1; b[i][j]=1; } //↖
else if (c[i-1][j]>=c[i][j-1])
{c[i][j]=c[i-1][j]; b[i][j]=2; } //↑
else { c[i][j]=c[i][j-1]; b[i][j]=3; } //←
}
}
构造最长公共子序列
void LCS(int i,int j,char x[]){
if (i ==0 || j==0) return;
if (b[i][j]== 1){
LCS(i-1,j-1,x);
cout<<x[i];
}
else if (b[i][j]== 2)
LCS(i-1,j,x);
else
LCS(i,j-1,x);
}