LIS和LCS-动态规划问题

题目描述

东东有两个序列A和B。

他想要知道序列A的LIS和序列AB的LCS的长度。

注意，LIS为严格递增的，即a1<a2<…<ak(ai<=1,000,000,000)。

Input

第一行两个数n，m（1<=n<=5,000,1<=m<=5,000）
第二行n个数，表示序列A
第三行m个数，表示序列B

Output

输出一行数据ans1和ans2，分别代表序列A的LIS和序列AB的LCS的长度

Simple Input

5 5
1 3 2 5 4
2 4 3 1 5

Simple Output

3 2

解题思路

本题可分为两部分，第一部分是计算序列A的LIS，即最长上升子序列。这是动态规划的典型题目，可以定义f1[i]是以A[i]结尾的最长上升序列的方程，那么初始化f1[1]=1，转移方程为f[i]=max(f[j] | j< i && A[j] < A[i]) +1。最后结果是max{f[i] , i=1,2,…,n}。
第二部分是LCS，即最长公共子序列，但此时并不要求是上升的。
可以定义f2[i][j]是以A[i]和B[j]结尾的最长公共子序列的方程。那么初始化为f2[i][0]=0,i=1,2,…n.f2[0][j]=0,j=1,2,…m.转移方程为当A[i]==B[j]时，f2[i][j]=f2[i-1][j-1]+1.否则f2[i][j]=max(f2[i][j-1],f2[i-1][j]).据此写出递推公式即可。最后的结果是f2[n][m].

代码

#include<iostream>
//#include<queue>
//#include<vector>
//#include<utility>
//#include<cstring>
using namespace std;
const int maxm=5e3+1;
int n,m; 
int f1[maxm];//序列A的LIS 
int A[maxm]; //A序列的值
int B[maxm];
int f2[maxm][maxm];//A，B序列的LCS
int ans1,ans2;
void lis()
{
	int sum,mmax=1;
	f1[1]=1;
	for(int i=2;i<=n;i++)
	{
		sum=0;
		for(int j=1;j<i;j++)
			if(A[j]<A[i] && f1[j]>sum)
			    sum=f1[j];
		f1[i]=sum+1;
		if(f1[i]>mmax) mmax=f1[i];
	}
	cout<<mmax<<" ";
}
void lcs()
{
	for(int i=1;i<=n;i++)
	    f2[i][0]=0;
   	for(int i=1;i<=m;i++)
	    f2[0][m]=0;
    for(int i=1;i<=n;i++)
    	for(int j=1;j<=m;j++)
	    	if(A[i]==B[j]) f2[i][j]=f2[i-1][j-1]+1;
	    	else f2[i][j]=f2[i][j-1]>f2[i-1][j]?f2[i][j-1]:f2[i-1][j];

	cout<<f2[n][m]<<endl;
}
int main()
{
	std::ios::sync_with_stdio(false);
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	    cin>>A[i];
	for(int i=1;i<=m;i++)
	    cin>>B[i];
    lis();
    lcs();
    return 0; 
}