求解最大子序列【动态规划】

最新推荐文章于 2021-05-22 22:29:28 发布

王亿亿

最新推荐文章于 2021-05-22 22:29:28 发布

阅读量180

点赞数

分类专栏：算法与数据结构

本文链接：https://blog.csdn.net/weixin_43931465/article/details/109089609

版权

动态规划最长公共子序列字符串处理算法 C++

关键词由CSDN通过智能技术生成

算法与数据结构专栏收录该内容

13 篇文章 1 订阅

订阅专栏

问题描述

字符序列的子序列是指从给定字符序列中随意地（不一定连续）去掉若干字符（可以不去）后所形成的字符序列。ac是abcd子序列。求两序列的最长公共子序列（LCS）。

动态规划

对于数组A[m],B[n]：

边界dp=0
Am-1=Bn-1,退一步，dp[m-1][n-1]=dp[m-2][n-2]
Am-1≠Bn-1,分别退一步，dp[m-1][n-1]=max(dp[m-2][n-1],dp[m-1][n-2])

# include<bits/stdc++.h>
using namespace std;
#define MAX 50
int dp[MAX][MAX];
vector<char> subs; //存放LCS
//动态规划dp数组 
void LCSLength(string a, string b, int m, int n)
{
	for(int i=0; i<=m; i++)
		dp[i][0]=0;
	for(int j=0; j<=n; j++) 
		dp[0][j]=0;
	for(int i=1; i<=m; i++)
		for(int j=1; j<=n; j++)
		{
			if(a[i]==b[j])
				dp[i][j]=dp[i-1][j-1]+1;
			else
				dp[i][j]=max(dp[i][j-1], dp[i-1][j]);
		}
} 

//从dp查找LCS
void BuildLCS(string a, string b, int m, int n)
{
	int k=dp[m][n];
	int i=m, j=n;
	while(k>0)
	{
		if(dp[i][j]==dp[i-1][j])
			i--;
		if(dp[i][j]==dp[i][j-1])
			j--;
		else
		{
			subs.push_back(a[i-1]);
			i--, j--, k--;
		}
	} 
}

int main()
{
	string a="abcdef";
	string b="hahabhabce";
	int m=a.length();
	int n=b.length();
	cout<<"序列a:"<<a<<endl;
	cout<<"序列b:"<<b<<endl;
	LCSLength(a, b, m, n);
	BuildLCS(a, b, m, n);
	cout<<"a,b最长公共子序列：";
	vector<char>::reverse_iterator rit;
	for(rit=subs.rbegin(); rit!=subs.rend(); ++rit)
		cout<<(*rit);
	cout<<endl;
	cout<<"长度为"<<dp[m][n]<<endl;
	
	return 0;
}