最长公共子序列LCS-基础动态规划问题-【C++】

最长公共子序列LCS

最长公共子序列也称作最长公共子串(不要求连续)，英文缩写为LCS（Longest Common Subsequence）。其定义是，一个序列 S ，如果分别是两个或多个已知序列的子序列，且是所有符合此条件序列中最长的，则 S 称为已知序列的最长公共子序列。

输入 
第一行给出一个整数N(0< N< 100)表示待测数据组数 
接下来每组数据两行，分别为待测的两组字符串。每个字符串长度不大于1000. 
输出 
每组测试数据输出一个整数，表示最长公共子序列长度。每组结果占一行。 

样例输入 
2 
asdf 
adfsd 
123abc 
abc123abc 
样例输出 
3 
6

解题思路：

经典的动态规划问题，需要先找出状态转移方程：
第一，如果str1[M]=str2[N]的情况下，最长公共子序列长度应该就是M-1、N-1长度的字符串组的最长公共子序列长度加上1。用公式表示就是LCS[M][N]=LCS[M-1][N-1]+1。
第二，如果str1[M]！=str2[N]的情况下，最长公共子序列长度应该在M-1、N 以及 M、N-1 这两个子问题之中最大公子序列长度的大者。用公式表示就是LCS[M][N]= MAX ( LCS[M-1][N] , LCS[M][N-1] )。

#include <iostream>
#include<bits/stdc++.h>
using namespace std; 

int LCS(string a,string b){
    int n,m;
    n=a.size();
    m=b.size();
    vector<vector<int> > memo= vector<vector<int> >(n+1,vector<int>(m+1,0));
    for(int i=1;i<n+1;i++)
        for(int j=1;j<m+1;j++){
            if(a[i]==b[j]){
                memo[i][j]=memo[i-1][j-1]+1;
            }
            else{
                memo[i][j]=memo[i-1][j]>memo[i][j-1]?memo[i-1][j]:memo[i][j-1];
            }
        } 
    return memo[n][m];
}

int main(int argc, char *argv[]) {
    string a="abcdefg";
    string b="acbebg";
    cout<<LCS(a,b);
    return 0;
}