迷之最长公共子序列（LCS）

最长公共子序列：

两个序列按顺序匹配相同的元素，所构成的新序列。(元素间可以不连续)
举例：

 a----b----e----d----c----f----g
 丨                  丨   丨   丨
ca-------------------c----fh---ge

最长公共子序列为 ：acfg

---

解决最长公共子序列：

1. 暴力搜索：假设长度为N的序列，其子序列个数（排列组合）：
   Cn(0)+Cn(1)+Cn(2)+Cn(3)+Cn(4)+……+Cn(n) = 2^n ；
   同理，另一长度为M的序列，其子序列个数为 2^m ；
   复杂度为O（2^(N+M)）；
   基本GG。
2. 动态规划：把求解问题分成若干个子问题。
   求前N-1元素与前M-1的最长公共字串；
   求前N-2元素与前M-2的最长公共字串；
   ……

---

如何使用动态规划解决最长公共子序列：

我们把两个序列从尾部一个一个拆分比较，会发现两种情况：
（设序列X,Y，X(a)为序列第a个元素）

1. 末尾元素相同，即X(N)=Y(M)。
   我们可以知道此时最长公共子序列的末尾元素一定也为X(N)(或Y（M）)。
   可继续把序列拆分成 比较X(N-1) 与 Y(M-1)。
2. 末尾元素不同，即X(N)≠Y(M)。
   此情况下又分两种情况：
   ——①最长子序列的末尾元素可能是X(N),非Y(M)。踢出Y(M), 比较X(N) 与 Y(M-1)
   ——②最长子序列的末尾元素可能是Y(M),非X(N)。踢出X(N), 比较X(N-1) 与 Y(M)

---

动态规划的递推结构：

放个模板：

void lcs(){
    la=strlen(a);
    lb=strlen(b);
    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-1][j],dp[i][j-1]);
            }
        }
    }
}