LCS板子(HDU-1159-Common Subsequence)
我又开始复习了~这个LCS,制作板子开始了。
dp[i][j]表示子序列Xi和Yj的最长公共子序列的长度。
当Xi == Yj时,找出Xi-1和Yj-1的最长公共子序列,然后在其尾部加上Xi即可得到X和Y的最长公共子序列;
当Xi != Yj时,求解两个子问题:1.求Xi-1和Yj的最长公共子序列 2.求Xi和Yj-1的最长公共子序列。然后取其最大值即可。
输出方案的话:和0/1背包输出的方案一样,LCS需要从后面倒推回去。
这里提到一点:因为涉及到下标-1的操作,我们需要特别注意边界情况,我们从dp[1][1]开始存储值,这样可以不会导致下标出界。
所以最后的答案存储在dp[X.length()][Y.length()]中在。
然后就是输出方案:
从后面开始判断:
1.判断dp[i][j]是否等于dp[i - 1][j],等于的话就i–
2.上面不成立,判断dp[i][j]是否等于dp[i][j - 1],等于的话就j–
3.以上都不成立,判断dp[i][j]是否大于dp[i - 1][j - 1],如果成立那么i–, j–,然后把X[i]加到我们的输出方案ans中。
最后我们要把ans给给翻转一下输出就可以啦
这里给出输出方案的板子吧,HDU-1159的题目删掉代码部分的print()就可以啦,这题是一道裸题吧。
代码部分:
#include <bits/stdc++.h>
#define mst(a, n) memset(a, n, sizeof(a))
using namespace std;
const int N = 1e3 + 10;
int dp[N][N];
string s1, s2;
string ans;
int LCS()
{
mst(dp, 0);
for (int i = 1; i <= s1.length(); i++)
{
for (int j = 1; j <= s2.length(); j++)
{
if (s1[i - 1] == s2[j - 1])
{
dp[i][j] = dp[i - 1][j - 1] + 1;
}
else
{
dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
}
}
}
return dp[s1.length()][s2.length()];
}
void print()
{
int i = s1.length(), j = s2.length();
while (dp[i][j])
{
if (dp[i][j] == dp[i - 1][j])
{
i--;
}
else if (dp[i][j] == dp[i][j - 1])
{
j--;
}
else if (dp[i][j] > dp[i - 1][j - 1])
{
i--, j--;
ans += s1[i];
}
}
reverse(ans.begin(), ans.end());
cout << ans << endl;
}
int main()
{
while (cin >> s1 >> s2)
{
ans = "";
cout << LCS() << endl;
print();
}
return 0;
}