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的长度

样例：
input:
5 5
1 3 2 5 4
2 4 3 1 5
output:
3 2

思路

• 两个题均是动态规划的问题
• 对于最大上升子序列问题，

•   状态：Fi为以Ai为结尾的最长子序列的长度
    初始化：F[1]=1
    状态转移方程： Fi=max(j<i&&A[j]<A[i])+1；
    结果：max{F[i],i=1,2,3...}
  

• 对于最大公共子序列的问题：

•   状态：F[i][j]为A1,A2,..Ai,B1,B2,..Bj的最大公共子序列的长度
    初始化：F[1][0]=F[0][0]=F[0][1]=0;
    状态转移方程：
    当A[i]=A[j]时，F[i][j]=F[i-1][j-1]+1;否则，F[i][j]=max{F[i-1][j],F[i][j-1]]}
    结果：F[n][m];
  

总结

对于动态规划的问题，首先对问题的分析之后，要进行状态的一个定义，然后转态转移方程的定义，其次就是就是初始化状态,最后也就是对结果的一个输出；

代码

#include<iostream>
#include<algorithm>
#include<string>
using namespace std;
int a[5010], b[5010];
int f[5010][5010];
int f1[5010];
int main()
{
	int n, m;
	cin >> n >> m;
	for (int i = 1; i <= n; i++)
		scanf("%d", &a[i]);
	for (int i = 1; i <= m; i++)
		scanf("%d", &b[i]);
	//求A得最长上升子序列
	memset(f1, 0, sizeof f1);
	f1[1] = 1;
	for (int i = 2; i <= n; i++) 
	{
		int max_ = 0;
		for (int j = 1; j < i; j++)
		{
			
			if (a[j] < a[i])
			{
				if (f1[j] > max_)
					max_= f1[j];
			}
		}
		f1[i] = max_ + 1;
	}
	sort(f1+1, f1 + n+1);
	int LIS = f1[n ];
	//求AB的最大公共子序列
	memset(f, 0, sizeof f);
	f[0][1] = f[0][0] = f[1][0] = 0;
	for(int i=1;i<=n;i++)
		for (int j = 1; j <=m; j++)
		{
			if (a[i] == b[j])
				f[i][j] = f[i - 1][j - 1] + 1;
			else
				f[i][j] = max(f[i][j - 1],f[i - 1][j]);
		}
	int LCS = f[n][m];
	cout << LIS << " " << LCS << endl;
	return 0;
}