LCS的加强输出 hdu1503

理解:

  mark数组:双重循环保证了两个字符串的每两个字符之间都进行了标记。标记了每两个字符之间不同的比较结果(相同,不相同且i这边是新增LCS(若dp相等记为新增空集),不相同且j这边是新增LCS)

  dp数组,记录了两个字符串分别到i和j的位置为止LCS多少长度。

  dp下标之所以从1开始不从0开始是因为自下而上要保证不会数组边界越出

  printLCS函数的原理在于mark已经记录了所有的比较情况,所以类似于一种状态转移。且要保证一个字符串当中的LCS字符等到另一个的LCS字符出现,一起输出。而非LCS字符总是这样的:如果字符串A先到了LCS的一个字符,它等待,而B串持续输出,直到它也出现这个字符,那就保证了在这期间输出的子串符合在B串中的顺序。

小了个结:

  第四点其实很重要,是一种acm的思想,就是一定有某种规律。

  

#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char str1[10000],str2[10000];
int len1,len2,mark[105][105],dp[105][105];
void LCS()
{
    int i,j;
    for(i=0;i<=len1;i++)
    {
        mark[i][0]=1;
    }
    for(i=0;i<=len2;i++)
    {
        mark[0][i]=2;
    }
    for(i=1;i<=len1;i++)
        for(j=1;j<=len2;j++)
    {
      if(str1[i-1]==str2[j-1])
            {
                dp[i][j]=dp[i-1][j-1]+1;
                mark[i][j]=0;
            }
      else
      {
          if(dp[i-1][j]>=dp[i][j-1])
          {
              mark[i][j]=1;//1表示i-1与j-1不同而且j-1是新增LCS
              dp[i][j]=dp[i-1][j];
          }
          else
          {
              mark[i][j]=2;//1表示i-1与j-1不同而且i-1是新增LCS
              dp[i][j]=dp[i][j-1];
          }
      }
    }
};
void PrintLCS(int i,int j)
{
    if(i==0&&j==0)
        return ;
    if(mark[i][j]==0)
    {
        PrintLCS(i-1,j-1);
        printf("%c",str1[i-1]);
    }
    else if(mark[i][j]==2)
    {
        PrintLCS(i,j-1);
        printf("%c",str2[j-1]);
    }
    else if(mark[i][j]==1)
    {
        PrintLCS(i-1,j);
        printf("%c",str1[i-1]);
    }
};
int main()
{
    while(scanf("%s %s",str1,str2)==2)
    {
        memset(dp,0,sizeof(dp));
        len1=strlen(str1);
        len2=strlen(str2);
        LCS();
        PrintLCS(len1,len2);
        putchar('\n');
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
最长公共子序列(LCS)是两个序列中都存在的最长子序列。为了修改算法 LCS 使其输出最长公共子序列,我们需要在 LCS 算法的基础上进行一些修改,具体步骤如下: 1. 构造一个二维数组 dp,其中 dp[i][j] 表示序列 A 中前 i 个元素和序列 B 中前 j 个元素的最长公共子序列长度。 2. 初始化 dp 数组的第一行和第一列,即 dp[0][j] 和 dp[i][0] 均为 0。 3. 对于序列 A 和序列 B 中的每一个元素,判断它们是否相等。 4. 如果 A[i] 和 B[j] 相等,则 dp[i][j] = dp[i-1][j-1] + 1。 5. 如果 A[i] 和 B[j] 不相等,则 dp[i][j] = max(dp[i-1][j], dp[i][j-1])。 6. 最终,最长公共子序列的长度为 dp[m][n],其中 m 和 n 分别为序列 A 和序列 B 的长度。 7. 根据 dp 数组求出最长公共子序列。从 dp[m][n] 开始,如果 A[i] 和 B[j] 相等,则将该元素加入最长公共子序列中,然后向左上方移动一格,即 i--,j--;如果不相等,则向较大的方向移动一格。 代码如下: ```python def LCS(A, B): m, n = len(A), len(B) dp = [[0] * (n+1) for _ in range(m+1)] # 构造二维数组 for i in range(1, m+1): for j in range(1, n+1): if A[i-1] == B[j-1]: dp[i][j] = dp[i-1][j-1] + 1 else: dp[i][j] = max(dp[i-1][j], dp[i][j-1]) # 求最长公共子序列 lcs = '' i, j = m, n while i > 0 and j > 0: if A[i-1] == B[j-1]: lcs = A[i-1] + lcs i -= 1 j -= 1 elif dp[i-1][j] > dp[i][j-1]: i -= 1 else: j -= 1 return lcs ``` 该算法的时间复杂度为 $O(mn)$,其中 m 和 n 分别为序列 A 和序列 B 的长度。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值