51Nod LCS 最长公共子序列 打印路径

题目

给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。
比如两个串为：

abcicba
abdkscab

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。
Input
第1行：字符串A
第2行：字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列，如果有多个，随意输出1个。
Sample Input
abcicba
abdkscab
Sample Output
abca

题目大意：
两个字符串中最长的公共子字符串，此串不一定连续。
设str1,str2, dp[i][j] 表示str1前i个字符与str2前j个字符的最长公共子字符串。
那么转移的方程是
if (str[i]==str[j]) dp[i][j] = dp[i][j]+1;
else dp[i][j] = max(dp[i-1][j],dp[i][j-1]);
这个题需要打印子字符串，那么在状态转移的时候，记录一下上一个状态。

易错点
第一次记录路径的时候。

 for ( int i=1; i<=la; i++ ) {
        for( int j=1; j<=lb; j++ ) {
            if ( A[i-1]==B[j-1] ) {
                dp[i][j] = dp[i-1][j-1] + 1 ;
            } else {
                dp[i][j] = max(dp[i][j-1],dp[i-1][j]);
            }
        }
        if ( cnt+1== dp[i][lb] ) {
            ans[cnt] = A[i-1];
            cnt++;
        }
    }

这种记录字符串的方式，错误的以为当前的状态会延续到最后的状态。
例如 str1 = “abcc” str2 = “accb”;
在i==2的时候，str1[1]==str2[3],这种情况下是dp[2][4]的最优解，但并不是全局的最优解。
看了网上的题解之后写出了正确的记录前一个状态的方法。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>

using namespace std;
const int N = 1000+10;
char A[N];
char B[N];
char ans[N];
int dp[N][N];
int path[N][N][2];

void print( int x,int y ) {
    if ( x==0 || y==0 )
        return;
    print( path[x][y][0],path[x][y][1] ) ;
    if ( A[x-1]==B[y-1] ) {
        printf("%c",A[x-1]);
    }
}
int main()
{
    cin>>A>>B;
    int la = strlen(A);
    int lb = strlen(B);
    int cnt =  0;

    // dp[i][j] 是 最长的公共字符串的长度
    for ( int i=1; i<=la; i++ ) {
        for( int j=1; j<=lb; j++ ) {
            if ( A[i-1]==B[j-1] ) {
                dp[i][j] = dp[i-1][j-1] + 1 ;
                path[i][j][0] = i-1;
                path[i][j][1] = j-1;
            } else {
                int index1 = dp[i][j-1];
                int index2 = dp[i-1][j];
                if (index1>index2) {
                    dp[i][j] = index1;
                    path[i][j][0] = i;
                    path[i][j][1] = j-1;
                } else {
                    dp[i][j] = index2;
                    path[i][j][0] = i-1;
                    path[i][j][1] = j;
                }
            }
        }

    }
    print( la,lb );
    return 0 ;
}