紫书第九章-----动态规划初步（最长公共子序列LCS）

最长公共子序列LCS

【题目 最长公共子序列Lcs 51Nod - 1006 】
给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。
比如两个串为：abcicba、abdkscab。ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。

解决方案
下面代码参考博客地址，特此说明！注意本题的打印路径还是非常巧妙的。我们可以这么思考，通过动态规划，不断地进行状态的转移，不断递推，不断更新dp数组，同时我们更新下面的pre数组，最后得到的dp[len1][len2]就是最大长度，更新pre为1的从len1,len2进行dfs回去不断搜索，遇到可以增长公共自序列的不断打印输出必定就是我们要求最长公共子序列（这一点和我们对状态的定义很有关系，因为我们定义的状态就是A字符串的前i个和B字符串的前j个的最长公共子序列）。

/*
    状态：dp[i][j]表示A字符串的前i个与B字符串的前j个的最长子序列
    状态转移：
        if(str[i]==str[j]) dp[i][j]=dp[i-1][j-1]+1
        else dp[i][j]=max(dp[i-1][j],dp[i][j-1])
    边界处理：dp[i][0]=0;dp[0][j]=0;
*/
#include<iostream>
#include<cstring>
using namespace std;

const int maxn = 1000+5;
int dp[maxn][maxn];
int pre[maxn][maxn];
char str1[maxn],str2[maxn];

void dfs(int i,int j){
    if(pre[i][j]==1){   //有状态转移才会让公共子序列的长度增长
        dfs(i-1,j-1);
        cout<<str1[i];
    }
    else if(pre[i][j]==2){
        dfs(i-1,j);
    }
    else if(pre[i][j]==3){
        dfs(i,j-1);
    }
    else{
        return ;
    }
}

int main()
{
    while(cin>>(str1+1)>>(str2+1)){
        int len1=strlen(str1+1);
        int len2=strlen(str2+1);
        memset(pre,0,sizeof(pre));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=len1;i++){
            for(int j=1;j<=len2;j++){
                if(str1[i]==str2[j]){
                    dp[i][j]=dp[i-1][j-1]+1;
                    pre[i][j]=1;//表示[i][j]可以往[i-1][j-1]转移
                }else{
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
                    if(dp[i-1][j]>dp[i][j-1]){
                        pre[i][j]=2;//表示[i][j]可以往[i-1][j]转移
                    }else{
                        pre[i][j]=3;//表示[i][j]可以往[i][j-1]转移
                    }
                }
            }
        }
        //注意dp[len1][len2]就是公共子序列的最大长度，可以如下打印输出
        //cout<<dp[len1][len2]<<endl;//本题不要求打印最大长度
        dfs(len1,len2);
    }
    return 0;
}