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

//动态规划法
//最长公共子序列
//c数组指出子序列长度，b数组为具体路径
#include <iostream>
#define N 1000
#define UPLEFT 1
#define UP 2
#define LEFT 3

int c[N][N],b[N][N];

using namespace std;

void PRINT_LCS(int b[][N],char X[],int i,int j)  //打印路径
{
    if(i==0 || j==0) return ;
    if(b[i][j] == UPLEFT)
    {
        PRINT_LCS(b,X,i-1,j-1);
        cout << X[i] << ends;
    }
    else if(b[i][j] == UP) PRINT_LCS(b,X,i-1,j);
    else PRINT_LCS(b,X,i,j-1);
}

LCS(char X[],char Y[],int m,int n)
{

    for(int i=0;i<=m;i++) c[i][0]=0;
    for(int j=0;j<=n;j++) c[0][j]=0;
    for(int i=0;i<=m;i++)
    {
        for(int j=0;j<=n;j++) b[i][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]=UPLEFT;
            }
            else if(c[i-1][j] >= c[i][j-1])
            {
                c[i][j]=c[i-1][j];
                b[i][j]=UP;
            }
            else
            {
                c[i][j]=c[i][j-1];
                b[i][j]=LEFT;
            }
        }
    }
    PRINT_LCS(b,X,m,n);
}

int main(int argc, char *argv[])
{
    int m,n; cin >> m >> n;

    char X[N],Y[N];

    for(int i=1;i<=m;i++) cin >> X[i];
    for(int i=1;i<=n;i++) cin >> Y[i];

    LCS(X,Y,m,n);
    cout << endl;
    cout << c[m][n] <<endl;
    cout << "Hello World!" << endl;
    return 0;
}