最长公共子序列

Problem Description

 给定两个序列X=

Input

输入数据有多组，每组有两行 ，每行为一个长度不超过500的字符串（输入全是大写英文字母（A,Z）），表示序列X和Y。

Output

每组输出一行，表示所求得的最长公共子序列的长度，若不存在公共子序列，则输出0。

Example Input

ABCBDAB
BDCABA

Example Output

4

分析：最长公共子序列的结构 ------ 最优子结构性质
           设序列 X=<x1，…，xm> 和 Y=<y1，…，yn>的一个最长公共子序列为 Z=<z1，…，zk>。

           则下述结论成立：
         1.若xm=yn，则zk=xm=yn 且Zk-1=<z1，…，zk-1>是 Xm-1和 Yn-1 的最长公共子序列。
          2. A、若xm！= yn 且zkxm，则 Z 是Xm-1 和 Y 的最长公共子序列。
              B、若xm！= yn 且zkyn，则 Z 是X 和 Yn-1 的最长公共子序列。
         结论：2个序列X、Y的最长公共子序列Z包含了这2个序列前缀的最长公共子序列 ---- 最长公共子序列问题具有最优子结构性质。

状态转移方程：c[i][j] = 0  (i=0或j=0)

                         c[i][j] = c[i-1][j-1] +1   (x[i] = y[i])

                         c[i][j] =  max(c[i-1][j] , c[i][j-1])    (x[i]!=y[i])

代码如下

#include<stdio.h>
#include<string.h>
int main()
{
    char x[505],y[505];
    int len1,len2,i,j;
    int c[505][505];
    while((gets(x)&&gets(y)))
    {
        memset(c,0,sizeof(c));
        len1 = strlen(x);
        len2 = strlen(y);
        for(i = 0;i<=len1;i++)
            c[i][0] = 0;//初始化
        for(j = 0;j<=len2;j++)
            c[0][j] = 0;//初始化首行首列为零 对应状态转移方程的第一种
        for(i = 1;i<=len1;i++)
        {
            for(j = 1;j<=len2;j++)
            {
                if(x[i-1]==y[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];
                }
            }
        }
        printf("%d\n",c[len1][len2]);

    }
    return 0;
}