最长上升子序列 （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的长度

模板题 没有侧式用例！！！！！！！！

思路:

最长上升子序列（LIS，Longest Increasing Subsequence）

定义状态i为以Ai为结尾得最长上升子序列方程，
第一个dp[1]=1;因为只有一个数，无关升降；
在dp[i]时要取[0…i）此区间内的已有的做大上升子序列，再+1，且此子序列满足此序列最大值小于A[i]。
状态转移方程：
dp[i]=max(dp[j] | j<i且Aj<Ai) + 1；
输出max(dp[i],i=1…n)

最长公共子序列（LCS，longest common subsequence）

假设A序列长度为n，B序列长度为m；
我们创建一个二维数组dp[n][m]为动态转移方程。
其中dp[i][j]代表A1,A2...Ai和B1,B2...Bj得LCS长度。
初始化:dp[0][0]=dp[1][0]=dp[0][1]=0(任意序列长度为0，则LCS为0)；
状态转移方程:
如果当前A，B序列得数字一样，则可以直接加入LCS:
当Ai==Bj时，dp[i][j]=dp[i-1][j-1]+1;
如果不一样，就取上一个状态得最优:
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);

代码:

#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn=5000+5;
const int maxm=5000+5;
int A[maxn],B[maxm];
int d[maxn];     //LIS
int dp[maxn][maxm];       //LCS
int n,m;
int ans1=0;

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]);
//LIS
	for(int i=1; i<=n; i++)
	{
		d[i]=1;
		for(int j=1; j<i; j++)
		{
			if(A[j]<A[i])
			{
				d[i]=max(d[i],d[j]+1);
			}
		}
		ans1=max(ans1,d[i]);
	}
//LCS  
	memset(dp,0,sizeof(dp));
	for(int i=1; i<=n; i++)
	{
		for(int j=1; j<=m; j++)
		{
			if(A[i]==B[j])  dp[i][j]=dp[i-1][j-1]+1;
			else
			{
				dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
			}
		}
	}
	printf("%d %d",ans1,dp[n][m]);
	return 0;
}

总结:

dp数组一定要注意其实和结束点，否则可能数组越界，或者状态转移不完整；