最长公共子序列问题 LCS

问题分析

这里分析了三种情况，主要是为了降低问题规模

递推公式

例子

创建一个二维数组，用来存储递推公式中的计算值。填充完数组。c[ i ] [ j ]=记录序列中最长公共子序列的长度。

/*
这里所有的想法都是针对一个LCS时候的情况。
包括下面的代码
*/

#include <stdio.h>
#include <string.h>
#define MAXLEN 100

/*
    *x,y:用来存储字符串的数组
    * m,n:分别是x,y字符串的长度
    * c[][MAXKEN]:存储LCS的二维数组
    * b[][MAXLEN]:存储c[][MAXlEN]数组中LCS的箭头指向
      (其中指向左上为0，指向上为1，指向左为-1)
    * 

*/
void LCSLength(char *x,char *y,int m,int n,int c[][MAXLEN],int b[][MAXLEN])
{
    int i,j;
    for(i=0;i<=m;i++)
        c[i][0]=0;
    for(j=1;j<=n;j++)
        b[0][j]=0;
    for(i=1;i<=m;i++)
    {
        for(j=1;j<=n;j++)
        {
            if(x[i-1]==y[j-1])
            {
                //递归公式中的第一种情况
                //箭头指向左上角
                c[i][j]=c[i-1][j-1]+1;
                b[i][j]=0;
            }
            else if(c[i-1][j]>=c[i][j-1])
            {
                //较大者为LCS，箭头向上
                //c[i-1][j]==c[i][j-1],箭头默认向上
                c[i][j]=c[i-1][j];
                b[i][j]=1;
            }
            else
            {
                //最后一种情况也就是C[i][j-1]>c[i-1][j]
                //这时候箭头向左指
                c[i][j]=c[i][j-1];
                b[i][j]=-1;
            }
        }
    }
}

void PrintLCS(int b[][MAXLEN],char *x,int i,int j)
{
    //用标记好的c和b数组
    //输出LCS子序列
    if(i==0 || j==0)//递归终止条件
        return;
    if(b[i][j]==0)//箭头指向左上角
    {
        PrintLCS(b,x,i-1,j-1);
        printf("%c", x[i-1]);
    }
    else if(b[i][j]==1)//箭头指向上
        PrintLCS(b,x,i-1,j);
    else//b[i][j]==-1,箭头指向左边
        PrintLCS(b,x,i,j-1);
}


int main()
{
    char x[MAXLEN];
    char y[MAXLEN];
    int m,n;
    scanf("%d %d", &m,&n);
    scanf("%s", x);
    scanf("%s", y);
    int b[MAXLEN][MAXLEN];
    int c[MAXLEN][MAXLEN];
    LCSLength(x,y,m,n,c,b);
    PrintLCS(b,x,m,n);
    return 0;
}