最长公共子序列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;
}