动态规划（四）最长公共子序列

###题目描述
最长公共子序列问题（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;
}