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

题目（示例）

int x[]={1,2,3,2,4,1,2};

int y[]={2,4};

求两个数组的最长公共子序列。

求解过程

思路（根据课本）

第一部分

其实看书还是非常好懂的，不过主要是

 if(x[i-1]==y[j-1]){
               
               c[i][j]=c[i-1][j-1]+1;
               b[i][j]=1;
            }

这块里面的x[i-1]和y[i-1]，因为刚开始我写的是x[i]==y[i],在运行结束后总是会发现少掉第一个，所以这里就改成x[i-1]了，其他就没什么问题，照着书上写就行

void LCS_LENGTH(int x[],int y[],int m,int n)

void LCS_LENGTH(int x[],int y[],int m,int n){
    int i,j;
    
    for(i=0;i<=m;i++) {
     b[i][0]=0;
     c[i][0]=0;
    }
    for(j=0;j<=n;j++) {
     b[0][j]=0;
     c[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]=1;
            }
            else if(c[i-1][j]>=c[i][j-1]){
               
                c[i][j]=c[i-1][j];
                b[i][j]=2;
             } 
            else{
                
                c[i][j]=c[i][j-1];
                b[i][j]=3;
             }
        }
    }
}

第二部分

void LCS(int i,int j,int x[])

void LCS(int i,int j,int x[]){
   
    if(i==0||j==0) return ;
    if(b[i][j]==1){
        LCS(i-1,j-1,x);
        // cout<<x[i];
        result[count++]=x[i-1];

        } else if(b[i][j]==2){
            LCS(i-1,j,x);
        } else {
            LCS(i,j-1,x);
        }

}

感觉也没啥好说的，因为c[i][j]最后一个肯定是该结果，也就是最长子序列的长度，我在主函数里面调用为

 LCS(m,n,x);

也就是从矩阵最下方开始：

• 如果是第一种情况，b[i][j]==1,那么就i-1,j-1，向斜对角线回移动;并输出（也可以把结果保存在一个数组里面）
• 第二种，就i-1,j,向上移动一行；
• 第三张，就i,j-1，向左移动一列;
• 一直到结束；

• 这个过程从书上那个矩阵图可以很直观的看出来；

运行结果如下

总结如下

其实感觉实现起来还是蛮简单的，但是当时一直结果不对，这种情况其实一般就是边界这块有问题，到底是i还是i-1，是>还是>=这些往往是需要我们自己注意的