动态规划-最长公共子序列

#include <iostream>
#define N 6
#define M 7
using namespace std;


int Lcs_Length(char x[],char y[],int b[M][N],int c[M+1][N+1])
{
    int i,j;
    for(i=0;i<=M;i++)
    {
        c[i][0]=0;
    }
    for(i=0;i<=N;i++)
    {
        c[0][i]=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]=3;
            }
            else if(c[i][j-1]>c[i-1][j])
            {
                c[i][j]=c[i][j-1];
                b[i][j]=2;
            }
            else
            {
                c[i][j]=c[i-1][j];
                b[i][j]=1;
            }
        }
    }
    return 0;
}
int Print_Lcs(int b[M][N],char x[],int i,int j)
{
    if(i==0||j==0)
        return 0;
    if(b[i][j]==3)
    {
        Print_Lcs(b,x,i-1,j-1);
        cout<<x[i-1];
    }
    else if(b[i][j]==2)
        Print_Lcs(b,x,i-1,j);
    else Print_Lcs(b,x,i,j-1);
    return 0;
}
int main()
{
    cout << "Hello world!" << endl;
    char x[M+1]="ABCBDAB";
    char y[N+1]="BDCABA";
    int b[M][N]={{0}};
    int c[M+1][N+1]={{0}};
    Lcs_Length(x,y,b,c);
    Print_Lcs(b,x,7,6);
    /*for (int i=0;i<N+1;i++)
    {
        for(int j=0;j<N+1;j++)
            cout<<m[i][j]<<' ';
        cout<<endl;
    }
    cout<<endl;
    for (int i=0;i<N+1;i++)
    {
        for(int j=0;j<N+1;j++)
            cout<<s[i][j]<<' ';
        cout<<endl;
    }*/
    return 0;
}