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

问题描述

最长公共子序列和最长公共字符串的概念

最优子结构

设序列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的最长公共子序列。

子问题的递归结构

代码

#include <iostream>
using namespace std;
#define MAXLEN 100
/*用c[i][j]记录序列Xi和Yj的最长公共子序列的长度
  用b[i][j]记录c[i][j]的值是由哪个子问题得到的
  m表示字符串x的长度，n表示字符串y的长度
*/
void LCSLength(int m,int n,char*x,char*y,int b[][MAXLEN],int c[][MAXLEN])
{
    for(int i=1;i<=m;i++)               //当j等于0时，空序列时Xi和Yj的最长公共子序列
       c[i][0]=0;    
    for(int j=1;j<=n;j++)               //当i等于0时，空序列时Xi和Yj的最长公共子序列
       c[0][j]=0;   
    for(int i=1;i<=m;i++)
       for(int j=1;j<=n;j++)
       {
           if(x[i]==y[j])
           {
               c[i][j]=c[i-1][j-1]+1;
               b[i][j]=1;                //1表示箭头指向左上方
           }
           else if(c[i-1][j]>=c[i][j-1])
           {
               c[i][j]=c[i-1][j];
               b[i][j]=2;               //2表示箭头指向上方
           }
           else
           {
               c[i][j]=c[i][j-1];
               b[i][j]=3;              //3表示箭头指向箭头指向左方
           }
           
       }    
}
void LCS(int i,int j,char *x,int b[][MAXLEN])
{
    if(i==0||j==0)
       return;
    if(b[i][j]==1)
    {
        LCS(i-1,j-1,x,b);
        cout<<x[i];
    }
    if(b[i][j]==2)
        LCS(i-1,j,x,b);
    if(b[i][j]==3)
        LCS(i,j-1,x,b);
}
int main()
{
   int m,n;
   char x[MAXLEN],y[MAXLEN];
   int b[MAXLEN][MAXLEN];
   int c[MAXLEN][MAXLEN];
   cin>>m>>n;
   for(int i=1;i<=m;i++)
      cin>>x[i];
   for(int i=1;i<=n;i++)
      cin>>y[i];
   LCSLength(m,n,x,y,b,c);
   LCS(m,n,x,b);
   return 0;
}