【Week10作业 B】LIS&LCS【动态规划】

题意：

东东有两个序列A和B。
他想要知道序列A的LIS和序列AB的LCS的长度。
注意，LIS为严格递增的，即a1<a2<…<ak(ai<=1,000,000,000)。
输入：第一行两个数n，m（1<=n<=5,000,1<=m<=5,000），第二行n个数，表示序列A，第三行m个数，表示序列B。
输出一行数据ans1和ans2，分别代表序列A的LIS和序列AB的LCS的长度。

---

思路：

首先求A的LIS最长上升子序列：定义fi表示以Ai为结尾的最长上升序列的方程，初始化f1=1，然后开始求解：fi=max{fj | j<i且Aj<Ai}+1。最终答案为max{f[i]}。
然后求AB的LCS最长公共子序列：定义f[i][j]为A1,A2,…,Ai和B1,B2,…,Bj的LCS长度。初始化f[1][0]=f[0][1]=f[0][0]=0。然后开始求解：遍历j=1->m,i=1->n,若Ai==Bj，则f[i][j]=f[i-1][j-1]+1；否则f[i][j]=max(f[i-1][j],f[i][j-1])。最终答案为f[n][m]。

---

总结：

一道动态规划题，要求掌握LIS以及LCS的思想和求解方法。

---

代码：

#include <iostream>
using namespace std;

int a[5010],b[5010];
int n,m;
int ans1=0,ans2;
int f1[5010];
int f2[5010][5010];
int main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
		cin>>a[i];
	for(int i=1;i<=m;i++)
		cin>>b[i];
	//求A的LIS
	f1[1]=1;
	for(int i=2;i<=n;i++)
	{
		int maxf=0;
		for(int j=1;j<i;j++)
			if(a[j]<a[i]&&f1[j]>maxf)
				maxf=f1[j];
		f1[i]=maxf+1;
	}
	for(int i=1;i<=n;i++)
		if(ans1<f1[i])
			ans1=f1[i];
	//求AB的LCS
	f2[1][0]=0,f2[0][1]=0,f2[0][0]=0;
	for(int j=1;j<=m;j++)
		for(int i=1;i<=n;i++)
			if(a[i]==b[j])
				f2[i][j]=f2[i-1][j-1]+1;
			else
				f2[i][j]=max(f2[i-1][j],f2[i][j-1]);
	ans2=f2[n][m];
	cout<<ans1<<" "<<ans2<<endl;
}