最长公共子序列（LCS）

给定两个字符串s1s2s3s4…sn和t1t2t3t4…tn。求出这两个字符串最长的公共子序列的长度。
比如s = “abcd“，t=”becd”
输出 3，最长公共子序列是 (“bcd”)

好吧，首先想到的就是暴力暴力暴力。。。暴力过。

然后就是动态规划。
dp[i][j]: s1…si和t1…tj对应的LCS的长度。
由此，s¬¬1…si+1和t1…tj+1对应的公共子序列可能是：
1. 当si+1 == tj+1时，直接在其LCS后加上si+1
2. 求s1… si+1 和t1… tj 最长子序列
3. 求s1… si 和t1。。。tj+1最长子列
这三点中的某一点，所以有如下递推关系成立：
dp[i+1][j+1]: 1. Max(dp[i][j]+1, dp[i+1][j], dp[i][j+1]) si+1 = tj+1
2. Max(dp[i][j+1], dp[i+1][j])

同时，使用数组temp[i][j]，记录每一次si和tj的关系，用于构造最长子序列。

直接看代码吧：

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

using namespace std;

const int MAX_N = 1010;
const int MAX_M = 1010;
int n, m;
int dp[MAX_N+1][MAX_M+1];
char s[MAX_N], t[MAX_M];
char temp[MAX_N+1][MAX_M+1];//记录子列相等时候的状态

void solve()//计算最长子序列的长度
{
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < m; j++)
        {
            if(s[i] == t[j])
            {
                dp[i+1][j+1] = dp[i][j] + 1;
                temp[i+1][j+1] = 0;
            }else{
                if(dp[i+1][j] >= dp[i][j+1])
                {
                    dp[i+1][j+1] = dp[i+1][j];
                    temp[i+1][j+1] = 1;
                }
                else{
                    dp[i+1][j+1] = dp[i][j+1];
                    temp[i+1][j+1] = 2;
                }
            }
        }
    }
    cout << dp[n][m] << endl;
}

//构造最长子序列
void LBSString(int i, int j)
{
    if(i == 0 || j == 0)
    {
        return ;
    }
    if(temp[i][j] == 0)
    {
        LBSString(i-1, j-1);
        printf("%c", s[i-1]);
    }else if(temp[i][j] == 1)
    {
        LBSString(i, j-1);
    }else{
        LBSString(i-1, j);
    }
}

int main()
{
    memset(dp, 0, sizeof(dp));
    cin >> s >> t;
    n = strlen(s);
    m = strlen(t);
    solve();
    LBSString(n, m);
    return 0;
}