动态规划求解LCS最长公共子序列问题c++

求最长公共子序列

最长公共子序列问题

（Longest-Common-Subsequence problem）给定两个序列 $X$ = <$x_1$，$x_2$，$\dots$，$x_m$> 和 $Y$= <$y_1$，$y_2$，$\dots$，$y_n$> ，求$X$和$Y$长度最长的公共子序列。

一个字符串 $s$ 被称作另一个字符串 $S$ 的子序列，说明从序列 $S$
通过去除某些元素但不破坏余下元素的相对位置（在前或在后）可得到序列 $s$ 。

例如：有两个字符串 $S$ 和 $T$ ，$S$ = “abcde”，$T$ = “acbfe”，那么$S$，$T$的最长公共子序列就是 “ace”，长度为3。

刻画最长公共子序列的特征

LCS的最优子结构：令$X$ = <$x_1$，$x_2$，$\dots$，$x_m$> 和 $Y$=
<$y_1$，$y_2$，$\dots$，$y_n$> 为两个序列，$Z$ =
<$z_1$，$z_2$，$\dots$，$z_k$> 为$X$和$Y$的任意$LCS$。

1. 如果 $x_m$ = $y_n$ ，则$z_k$ = $x_m$ = $y_n$，且$Z_{k-1}$是$X_{m-1}$和$Y_{n-1}$的一个$LCS$。
2. 如果$x_m$ $\ne$ $y_n$，那么$z_k$ $\ne$ $x_m$意味着$Z$是$X_{m-1}$ 和 $Y$ 的一个$LCS$。
3. 如果$x_m$ $\ne$ $y_n$，那么$z_k$ $\ne$ $y_n$意味着$Z$是$X$ 和 $Y_{n-1}$ 的一个$LCS$。

递归公式

根据 $LCS$ 问题的最优子结构性质，可以得到如下公式：

其中dp[ i ] [ j ]表示 $S{1}_i$ 和 $S{2}_j$ 的$LCS$长度。

C++代码

1. 求最长公共子序列长度

int dp[1010][1010];
int main()
{
	string s1,s2;
	cin>>s1>>s2
	int len1=s1.length();
	int len2=s2.length();
	memset(dp,0,sizeof(dp));
	for(int i=1;i<=len1;i++)
		for(int j=1;j<=len2;j++)
		{
			if(s1[i-1]==s2[j-1])
				dp[i][j]=dp[i-1][j-1]+1;
			else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
		}
	cout<<dp[len1][len2]<<endl;
}

2. 打印任意一组最长公共子序列

int dp[3010][3010];
int b[3010][3010];//记录路径  
void print_LCS(string s,int i,int j)
{
	if(i==0 || j==0) return;
	if(b[i][j]==0)
	{
		print_LCS(s,i-1,j-1);
		cout<<s[i-1];                                                     
	}
	else if(b[i][j]==-1)
	{
		print_LCS(s,i-1,j);
	}
	else print_LCS(s,i,j-1);
}
int main()
{
	string s,t;
	cin>>s>>t;
	int len1=s.length();
	int len2=t.length();
	memset(dp,0,sizeof(dp));
	for(int i=1;i<=len1;i++)
		for(int j=1;j<=len2;j++)
		{
			if(s[i-1]==t[j-1])
			{
				dp[i][j]=dp[i-1][j-1]+1;
				b[i][j]=0;//0表示指向左上方 
			}
			else if(dp[i-1][j]>=dp[i][j-1])
			{
				dp[i][j]=dp[i-1][j];
				b[i][j]=-1;//-1表示指向上方 
			}
			else
			{
				dp[i][j]=dp[i][j-1];
				b[i][j]=1;//1表示指向左方 
			}
		}
	print_LCS(s,len1,len2);
	return 0;
}

模板题

求最长公共子序列的长度

poj1458 Common Subsequence

打印最长公共子序列

luogu AT4527 LCS