问题提出:
例如,若X= { A, B, C, B, D, A, B}和Y= {B, D, C, A, B, A},则序列{B,C,A}是X和Y的一个公共子序列,序列{B,C,B,A}也是X和Y的一个公共子序列。但是,后者是X和Y的一个最长公共子序列。
问题分析:
设序列X={x1,x2,……xm}和Y={y1,y2,……yn}的最长公共子序列为Z={z1,z2,……zk}。则有:
-
若xm=yn,则zk=xm=yn,且zk-1是Xm-1和Yn-1的最长公共子序列。
-
若xm!=yn且zk!=xm,则Z是Xm-1和Y的最长公共子序列。
-
若xm!=yn且zk!=yn,则Z是X和Yn-1的最长公共子序列。
其中,Xm-1={x1,x2……xm-1},Yn-1={y1,y2……yn-1},Zk-1={z1,z2……zk-1}。
所以,最长公共子序列问题具有最优子结构性质。
子问题的递归结构:
在计算X和Y的最长公共子序列时,可能要计算[X和Yn-1]及[Xm-1和Y]的最长公共子序列。这两个子问题都包含一个公共子问题,即计算Xm-1和Yn-1的最长公共子序列。
用 c[i][j] 记录序列 Xi 和 Yj 的最长公共子序列的长度。其中,Xi = {x1,x2……xi},Yj = {y1, y2……yj}。当 i = 0 或 j = 0 时,空序列是xi和yj的最长公共子序列。此时,c[i][j] = 0。其他情况如下:
代码如下:
#include<iostream>
#include<vector>
using namespace std;
const int M = 8;
const int N = 7;
void LCSLength(int m, int n, char x[], char y[], vector< vector<int> > &c){
int i, j;
for(i = 1; i < m; i++)
//y为空串
c[i][0] = 0;
for