【Week 10 作业B】LIS&LCS

题目描述

东东有两个序列A和B。
他想要知道序列A的LIS和序列AB的LCS的长度。
注意，LIS为严格递增的，即a1<a2<…<ak(ai<=1,000,000,000)。

LIS

LCS

输入格式

第一行两个数n，m（1<=n<=5,000,1<=m<=5,000）

第二行n个数，表示序列A

第三行m个数，表示序列B

输出格式

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

思路

LIS

对于一个序列，要选择长度最长的上升序列，这里选择动态规划的思想，首先约定一个状态fi表示以序列第i个数Ai结尾的最长上升序列长度。
易知，f1=1，对于任意fi，fi的值为max{fj}+1，其中j<i且Aj<Ai，若fj不存在，则fi=1，这样依次计算fi，直到以序列尾为结尾的序列即fn计算完，这时答案即为max{fi},其中0<i<=n（n为序列长度）。

LCS

这里依然采用动态规划的思想，约定状态fij为序列A1,A2……Ai和B1，B2……Bj的LCS长度，易知，f0*=f*0=f10=0,对于其他状态，若Ai==Bj则fij=f(i-1)(j-1)+1,否则fij=max(f(i-1)j,fi(j-1))。
所有状态fij（0<i<=n,0<j<=m）计算完成后，fnm即为所求。

代码

#include <iostream>
#include <string.h>
using namespace std;
const int size=5000+10;
int A[size];
int B[size];
int cofa[size];
int cofab[size][size];
int n,m;
int lis()
{
	int max=1;
	cofa[1]=1;
	for(int i=2;i<=n;i++)
	{
		cofa[i]=1;
		for(int j=i-1;j>0;j--)
		{
			if(A[j]<A[i])
				cofa[i]=cofa[i]>cofa[j]+1?cofa[i]:(cofa[j]+1);
		}
		if(cofa[i]>cofa[max])
			max=i;
	}
	return max;
}
int lcs()
{
	memset(cofab,0,sizeof(cofab));
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(A[i]==B[j])
				cofab[i][j]=cofab[i-1][j-1]+1;
			else
				cofab[i][j]=max(cofab[i-1][j],cofab[i][j-1]);
		}
	}
}
int main() {
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%d",&A[i]);
	for(int i=1;i<=m;i++)
		scanf("%d",&B[i]);
	int max=lis();
	lcs();
	printf("%d %d\n",cofa[max],cofab[n][m]);
	return 0;
}