LIS and 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

解题思路
LIS: 最长上升子序列
例：C= 1 2 3 4 5 0 9
最长上升序列为 1 2 3 4 5 9

LCS: 最长公共子序列
例： A = 1 5 2 6 8 7
B = 2 3 5 6 9 8 4
LCS 为5 6 8或 2 6 8

求LIS的长度：
两重循环，外层循环遍历C,从0到n
内层循环从0到 i
cnt数组记录[ 0, i ] 的最长上升子序列
也就是说，c[0,4] 记录了 到4 的 LIS的长度
cnt[i] = max ( cnt[j] + 1, cnt[i] )
遍历到C[j] ,如果从 c[j] 直接跳到了c[i] ,那么长度为多少 , 然后和 cnt[i]比较，取最大值

LCS的长度:
一维数组无法记录全部信息，所以有二维数组记录 f [] []

同样2重循环，更新f[i][j] ,表示 A [ 1 , i ] , B [ 1 , j ] 的最长公共子序列的长度

代码实现

#include<iostream>
using namespace std;
#include<algorithm>
#define ll long long

ll a[5005];
ll b[5005];

int cnt[5005];
int f[5005][5005];

int n,m;
int lis();
int lcs();
int main()
{
	ios::sync_with_stdio(0);
	cin>>n>>m;
	for(int i=0;i<n;i++)
	{cin>>a[i];cnt[i]=1;}
	for(int i=0;i<m;i++)
	cin>>b[i];
	
	int res1,res2;
	res1=lis();
	res2=lcs();
	cout<<res1<<" "<<res2;
	return 0;
}

int lis()
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<i;j++)
		{
			if(a[j]<a[i])
			{
				cnt[i] = max(cnt[i],cnt[j]+1);
			}
		}
	}
	int res=-1;
	for(int i=0;i<n;i++)
	{
		if(res < cnt[i])
		res = cnt[i];
	}
	return res;
}

int lcs() //从1到n，m的 
{
	for(int i=0;i<n;i++) cnt[i]=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(a[i-1]==b[j-1]) f[i][j] = f[i-1][j-1]+1;
			else f[i][j] = max(f[i-1][j],f[i][j-1]);
		}	
	}
	return f[n][m];
}

小结
动态规划的代码不长，但是确实非常考验思维，只要找到状态转移方程，就不难了。