###题目描述
最长公共子序列问题(LCS),给定两个序列X= {x1,x2…,xm},Y= {y1,y2,…yn},求X 和Y长度最长的公共子序列。
什么叫公共子序列:Z既是X的子序列,又是Y的子序列,我们称它是X和Y的公共子序列。
子序列和子串的区别:子串是串的一个连续的部分,子序列是不改变顺序,从序列中去掉任意的元素而获得的新序列。
例如:X = {A,B,C,B,D,A,B},Y= {B,D,C,A,B,A}
输出为 4,Z= {B,C,B,A }为最长的公共子序列。
###分析
####1.刻画最长公共子序列的特征
子问题的自然分类对应两个输入序列的前缀对,前缀的定义为:给定一个序列X= {x1,x2…xm},对i= (0,m),定义X的第i前缀为Xi = {x1,x2…xi}。其实就是以xi为后缀
令X= (x1,x2,…xm)和Y= {y1,y2,…yn}为两个序列,Z= {z1,z2,…zk}为X 和Y的任意LCS。
1,如果xm = yn,则zk = xm= yn,且Zk-1 是Xm-1和Yn-1的一个LCS.
2,如果xm!= yn,zk != xm,意味着Z是Xm-1和Yn的一个LCS
3,如果xm != yn, zk != yn,意味着Z是X和Yn-1的一个LCS
两个序列的LCS包含两个序列的前缀的LCS,因此,LCS问题具有最优子结构性质。
####2.建立递归式
我们可一很容易看出LCS问题的重叠子问题性质,为了求X和Y的一个LCS,我们可能需要求X和Yn-1的一个LCS及Xm-1和Y的一个LCS,但是这几个问题都包含求解Xm-1和Yn-1的LCS的子子问题。
建立最优解的递归式
定义c[i][j]表示Xi和Yj的LCS的长度。
c[i][j] = 0 i = 0或j= 0 i从1开始到m
c[i][j] = c[i-1][j-1]+1 ;xi = yj
c[i][j] = max(c[i-1][j],c[i,j-1]) ;xi!= yj
我们通过限制条件限定了需要求解哪些子问题。
####3.动态规划方法
#include<stdio.h>
#include <stdlib.h>
#include <string.h>
int lcs(char * a,char * b)
{
int m = strlen(a);
int n = strlen(b);
int (*c)[n+1] = (int (*)[n+1])malloc((m+1) *(n+1)* sizeof(int));
int i,j;
for(i = 0;i<=m;i++)
c[i][0] = 0;
for(j = 0;j<=n;j++)
c[0][j] = 0;
for(i = 1;i<= m;i++) //对应从x1-xm,x0 = 0
{
for(j = 1;j<= n;j++)
{
if(a[i-1] == b[j-1])
c[i][j] = c[i-1][j-1]+1;
else
{
if(c[i-1][j] > c[i][j-1])
c[i][j] = c[i-1][j];
else
c[i][j]= c[i][j-1];
}
}
}
int result =c[m][n];
free(c);
return result;
}
int main()
{
char a[10] = "abcbdab";
char b[10] = "bdcaba";
printf("%d\n",lcs(a,b));
return 0;
}