LCS最长公共序列是很普遍的动态规划问题,描述如下:
给定二个字符串s1,s2,求出最长的公共子序列,需要注意的是这里的公共子序列不要求是连续的
例如:s1="abcad",s2="daba",这里的"aba"就是公共子序列,所以长度为3,还有有可能会有很多个长度为3的子序列。
解题思路:这里需要用到二维整型数组dp[i][j]表示的是从s1中前i个字符串和s2中前j个字符串的LCS,我们可以这样递推:
1,当知道dp[i][j]的时候,如果s1[i+1]==s2[j+1],dp[i+1][j+1]=s[i][j]+1; 否则的话,dp[i+1][j+1]=max(dp[i+1][j],dp[i][j+1])
2,这里需要注意的是dp[][]的初始化问题,这里的dp[i][j]中的i和j表示的是s1和s2中字符串的长度分别为i,j。所以,不管任何一个子序列长度为0的时候,dp[][]的值都为0,即dp[0][i]=0;dp[i][0]=0;然后从长度为1开始计算dp[1][1]的值,以此类推,最后dp[n][n]的值就是答案。
代码:
#include<iostream>
#include<string>
using namespace std;
string s1, s2;
int n, m;
int c[1000][1000];
int max(int a, int b)
{
if (a >= b) return a;
else return b;
}
int main()
{
cin >> s1 >> s2;
n = s1.length();
m = s2.length();
for (int i = 0; i <= n; i++)//dp[][]数组的初始化
c[i][0] = 0;
for (int j = 0; j <= m; j++)//dp[][]数组的初始化
c[0][j] = 0;
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= m; j++)
{
if (s1[i - 1] == s2[j - 1])//因为dp[i][j]是指长度为i和j,但是字符串是从0开始的,所以用i-1,j-1.
{
c[i][j] = c[i - 1][j - 1] + 1;
}
else
{
int ma = max(c[i][j - 1], c[i - 1][j]);
c[i][j] = ma;
}
}
}
cout << c[n][m] << endl;
return 0;
}
不知道写的容不容易被理解,有问题欢迎评论哦~ 么么哒