【动态规划T】最长公共子序列

最新推荐文章于 2023-11-14 17:19:36 发布

原创最新推荐文章于 2023-11-14 17:19:36 发布 · 672 阅读

0 ·

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

标签

#c++ #算法 #动态规划

算法专栏收录该内容

68 篇文章

订阅专栏

该博客探讨了计算机科学中的一个经典问题——最长公共子序列（LCS）算法，给出了具体的算法分析和C++实现。通过示例输入和输出展示了如何找到两个序列的LCS，并解释了算法的时间复杂度。LCS在序列比对、文本处理等领域有广泛应用。

【问题描述】

若给定序列X={x1,x2,…,xm}，则另一序列Z={z1,z2,…,zk}，是X的子序列是指存在一个严格递增下标序列{i1,i2,…,ik}使得对于所有j=1,2,…,k有：zj=xij。
例如，序列Z={B，C，D，B}是序列X={A，B，C，B，D，A，B}的子序列。
给定2个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是序列X和Y的公共子序列。
给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}，找出X和Y的最长公共子序列。

输入
三行.第一行m和n,分别表示两个子序列的长度（符号的个数）
接下来两行，分别表示X和Y
输出
最长公共子序列的长度

**样例输入**
7 6
ABCBDAB
BDCABA
**样例输出**
4

【算法分析】

最长公共子序列（线性DP）

【算法实现】

#include<iostream>
#include<cstdio>

using namespace std;
char a[1001],b[1001];
int f[1001][1001];
int main()
{
   int n1,n2;
   cin>>n1>>n2;
   for(int i=0;i<n1;i++)
   {
       cin>>a[i];
   }
      for(int i=0;i<n2;i++)
   {
       cin>>b[i];
   }
   for(int i=0;i<=n1;i++)
   {
       f[i][0]=0;
   }
   for(int i=0;i<=n2;i++)
   {
       f[0][i]=0;
   }
   for(int i=1;i<=n1;i++)
   {
       for(int j=1;j<=n2;j++)
       {
           if(a[i]==b[j])
              f[i][j]=f[i-1][j-1]+1;
           else
               f[i][j]=max(f[i-1][j],f[i][j-1]);
       }
   }
   cout<<f[n1][n2]<<endl;
   return 0;
}