「这是我参与2022首次更文挑战的第24天,活动详情查看:2022首次更文挑战」
前言
今天分享一道51NOD上的算法题,也是以前写的一道比较经典的算法题,下面回顾一下,再从新做一遍。
题目描述
给出两个字符串A B,求A与B的最长公共子序列(子序列不要求是连续的)。
比如两个串为:
abcicba abdkscab
ab是两个串的子序列,abc也是,abca也是,其中abca是这两个字符串最长的子序列。 输入
第1行:字符串A 第2行:字符串B (A,B的长度 <= 1000)
输出
输出最长的子序列,如果有多个,随意输出1个。
输入示例
abcicba abdkscab
输出示例
abca
解题思路
引进一个二维数组c[][],用c[i][j]记录X[i]与Y[j] 的LCS 的长度,b[i][j]记录c[i][j]是通过哪一个子问题的值求得的,以决定搜索的方向。 我们是自底向上进行递推计算,那么在计算c[i,j]之前,c[i-1][j-1],c[i-1][j]与c[i][j-1]均已计算出来。此时我们根据X[i] = Y[j]还是X[i] != Y[j],就可以计算出c[i][j]。
问题的递归式为:
回溯输出的流程图如下过程:
时间复杂度分析:
由于每次调用至少向上或向左(或向上向左同时)移动一步,故最多调用(m + n)次就会遇到i = 0或j = 0的情况,此时开始返回。返回时与递归调用时方向相反,步数相同,故算法时间复杂度为o(m + n)。
具体代码如下
#include<iostream>
#include<iomanip>
#include <string.h>
using namespace std;
long int c[1001][1001];
int main()
{
string a,b;
char d[1000];
long long int n1,n2,i,j,k;
cin>>a>>b;
//n1=strlen(a);strlen(用于数组中)
//n2=strlen(b);
n1=a.size();
n2=b.size();
k=0;
for(i=1;i<=1000;i++){c[0][i]=0;c[i][0]=0;}
for(i=1;i<=n1;i++)
for(j=1;j<=n2;j++)
c[i][j]=(a[i-1]==b[j-1])?(c[i-1][j-1]+1):max(c[i][j-1],c[i-1][j]);
for(i=n1,j=n2;i>=1&&j>=1;)
{
if(a[i-1]==b[j-1])
{
d[k]=a[i-1];
k++;
i--;
j--;
}
else
{
if(c[i][j-1]>c[i-1][j])//也可以等于,无所谓,等于的话回溯路线会有所改变
{
j--;
}
else
{
i--;
}
}
}
for(i=k-1;i>=0;i--)cout<<d[i];
cout<<endl;
return 0;
}
总结
算法题最重要的是思路,有时候即使做过的题下次遇到可能还是不会,要多多练习,多接触一些题目,俗话说的好,温故而知新,没事得多看看以前做得记录,gogogo