C++:最长公共子序列问题(动态规划)

Time Limit: 1000 ms Memory Limit: 65536 KiB

Submit Statistic

Problem Description

给定两个序列 X={x1,x2,…,xm} 和 Y={y1,y2,…,yn},找出X和Y的最长公共子序列。

Input

输入数据有多组,每组有两行 ,每行为一个长度不超过500的字符串(输入全是大写英文字母(A,Z)),表示序列X和Y。

Output

每组输出一行,表示所求得的最长公共子序列的长度,若不存在公共子序列,则输出0。

Sample Input

ABCBDAB
BDCABA

Sample Output

4

Hint

Source

#include <bits/stdc++.h>

using namespace std;

int main()
{
    char s1[510], s2[510];
    while (cin >> s1 >> s2)
    {
        int m = strlen(s1);
        int n = strlen(s2);
        int c[m + 1][n + 1];
        for (int i = 0; i <= m; i++)
        {
            for (int j = 0; j <= n; j++)
            {
                c[i][j] = 0;
            }
        }
        for (int i = 1; i <= m; i++)
        {
            for (int j = 1; j <= n; j++)
            {
                if (s1[i - 1] == s2[j - 1])
                    c[i][j] = c[i - 1][j - 1] + 1;
                else
                    c[i][j] = max(c[i - 1][j], c[i][j - 1]);
            }
        }
        cout << c[m][n] << endl;
    }
}

 

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
公共序列问题是一种经典的动态规划问题,它可以用来求两个字符串之间的最公共序列。以下是该问题动态规划解法。 我们令 $dp[i][j]$ 表示字符串 $s1$ 的前 $i$ 个字符和字符串 $s2$ 的前 $j$ 个字符的最公共序列度,那么状态转移方程为: $$ dp[i][j] = \begin{cases} 0, & i=0 \text{ 或 } j=0 \\ dp[i-1][j-1]+1, & s1[i]=s2[j] \\ \max(dp[i][j-1], dp[i-1][j]), & s1[i] \neq s2[j] \end{cases} $$ 其中,当 $s1[i]=s2[j]$ 时,$dp[i][j]$ 等于 $s1$ 的前 $i-1$ 个字符和 $s2$ 的前 $j-1$ 个字符的最公共序列度 $dp[i-1][j-1]$ 加上当前字符相同的贡献 $1$。 当 $s1[i] \neq s2[j]$ 时,$dp[i][j]$ 取 $s1$ 的前 $i$ 个字符和 $s2$ 的前 $j-1$ 个字符的最公共序列度 $dp[i][j-1]$ 和 $s1$ 的前 $i-1$ 个字符和 $s2$ 的前 $j$ 个字符的最公共序列度 $dp[i-1][j]$ 的较大值。 最终结果即为 $dp[n][m]$,其中 $n$ 和 $m$ 分别为字符串 $s1$ 和 $s2$ 的度。 以下是该问题的 C 代码实现: ```c #include <stdio.h> #include <string.h> #define MAX_N 1000 int dp[MAX_N + 1][MAX_N + 1]; int lcs(char *s1, char *s2) { int n = strlen(s1), m = strlen(s2); for (int i = 0; i <= n; ++i) dp[i][0] = 0; for (int j = 0; j <= m; ++j) dp[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] = dp[i - 1][j - 1] + 1; } else { dp[i][j] = dp[i][j - 1] > dp[i - 1][j] ? dp[i][j - 1] : dp[i - 1][j]; } } } return dp[n][m]; } int main() { char s1[MAX_N + 1], s2[MAX_N + 1]; scanf("%s%s", s1, s2); printf("%d\n", lcs(s1, s2)); return 0; } ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值