1306：最长公共子上升序列

时间限制: 1000 ms         内存限制: 65536 KB
【题目描述】
给定两个整数序列，写一个程序求它们的最长上升公共子序列。
当以下条件满足的时候，我们将长度N的序列S1,S2,...,SN 称为长度为M的序列A1,A2,...,AM的上升子序列：
存在1≤i1<i2<...<iN≤M，使得对所有1≤j≤N，均有Sj=Aij，且对于所有的1≤j<N，均有Sj<Sj+1。
【输入】
每个序列用两行表示，第一行是长度M(1≤M≤500)，第二行是该序列的M个整数Ai(−2^31<=Ai<2^31)
【输出】
在第一行，输出两个序列的最长上升公共子序列的长度L。在第二行，输出该子序列。如果有不止一个符合条件的子序列，则输出任何一个即可。
【输入样例】
5
1 4 2 5 -12
4
-12 1 2 4
【输出样例】
2
1 4
【提示】
经典算法Baidu搜索，深刻体会。

#include<iostream>
#include<algorithm>
#include<string.h>
#include<string>
#define N 505
using namespace std;
long long a[N],b[N],k[N];
int dp[N][N]={0},p[N][N];
int main()
{
	int n,m,i,j,ans,mx,y;
	cin >> m;
	for (i = 1;i <= m;i++)cin >> b[i];
	cin >> n;
	for (i = 1;i <= n;i++)cin >> a[i];
	ans = 0;
	for (i = 1;i <= n;i++)
	{
		mx = 0;y = 0;
		for (j = 1;j <= m;j++)
		{
			dp[i][j] = dp[i-1][j];//保证到i、j为止的最大长度不可能小于i-1，j的
			p[i][j] = j;//只要以一个序列为准记住p就可以
			//i不变，遍历j，此时要保证j序列值要小才可以是上升序列
			if (a[i] > b[j] && mx < dp[i-1][j])
			{	
				mx = dp[i-1][j];//记录dp[i-1][0~j]的最大值
	            y = j;//记录此时的j值
			}
			else if (a[i] == b[j]) 
			{
				dp[i][j] = mx+1;
				p[i][j] = y;
			}
			ans = max(ans,dp[i][j]);
		}
	}
	cout << ans << endl;
	int T = 1;
	for (j = 2;j <= m;j++)//看作一位数组，找最大的序列的末尾为止，也可以在上面的ans更新
		if (dp[n][j] > dp[n][T])
			T = j;
	int top = 0;
	i = n;j = T;
	while (dp[i][j])
    {
		 while(a[i] != b[j] && i) i--;//找到一个相等的值为止
		 k[top++] = b[j];
		 j = p[i][j];
    }
	for (i = top-1;i >= 0;i--)
		cout << k[i] << " ";
	cout << endl;
	return 0;
}