POJ-1458(Common Subsequence )

最新推荐文章于 2024-07-13 16:30:23 发布

原创最新推荐文章于 2024-07-13 16:30:23 发布 · 332 阅读

2 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#动态规划 #算法

打卡专栏收录该内容

13 篇文章

订阅专栏

本文详细解析了求解两个字符串最长公共子序列的问题，通过动态规划算法实现，阐述了状态转移方程，并提供了完整的代码示例。

题目大意

求两个字符串的最长公共子序列。如"abcd"和“abfceb”的最长公共子序列为"abc"。

解题思路

动态规划的经典问题：最长公共子序列问题。
假设str1表示输入的第一个序列，str2表示输入的第二个序列，i表示指向str1的指针,j表示指向str2的指针。
1.状态集合表示：dp[i][j]表示第一个序列的前i个字母和第2个序列的前j个字母的最长公共子序列长度
2.状态转移
–>a.选str1[i],不选str2[j], dp[i][j] = dp[i][j-1]。需要注意的是dp[i][j-1]并不等价于选str1[i],不选str2[j]，因为dp[i][j-1]表示的是str1的前i个字符和str2的前j-1个字符中的最长公共子序列长度，所有有可能包含str1[i]，也有可能不包含。即dp[i][j-1]表示的范围包含了该种情况。情况b也是如此。
–>b.不选str1[i],选str2[j],dp[i][j] = dp[i-1][j]
–>c.选str1[i]且选str2[j],dp[i][j] = dp[i-1][j-1] + 1。两者都选的前提是str1[i] == str2[j]
–>d.不选str1[i]且不选str2[j],dp[i][j] = dp[i-1][j-1]。该种情况被a,b所包含了，可以不用求。
dp[i][j] = max{a, b, c, d}

#include <iostream>
#include <string>
#include <algorithm>
#include <cstring>
using namespace std;
const int MAXN = 1005;
int dp[MAXN][MAXN];     //dp[i[j]表示第一个序列的前i个字母和第2个序列的前j个字母的最长公共子序列长度
int main() {
    string str1, str2;
    while(cin >> str1 >> str2) {
        memset(dp, 0, sizeof(dp));
        int n = str1.length();
        int m = str2.length();
        for(int i = 1; i <= n; i++) {   //注意:下标从1开始(避免状态转移时数组越界),取字符时直接-1即可
            for(int j = 1; j <= m; j++) {
               dp[i][j] = max(dp[i-1][j], dp[i][j-1]);
               if(str1[i-1] == str2[j-1]) {
                   dp[i][j] = max(dp[i][j], dp[i-1][j-1] + 1);
               }
            }
        }
        cout << dp[n][m] << endl;;
    }
    return 0;
}