week10作业题_B - LIS & LCS

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

Simple Input
5 5
1 3 2 5 4
2 4 3 1 5
Simple Output
3 2

题目思路

求LIS（最长上升子序列）时，用两个for循环来遍历，一个用来遍历序列中的所有元素，一个用来遍历外层循环之前的所有元素。
状态：定义fi表示以Ai为结尾的最长上升序列的方程。
初始化：f1 = 1
转移过程：

求LCS（最长公共子序列）时
设计状态：假设 f[i][j] 为 A1, A2, …, Ai 和 B1, B2, …, Bj 的 LCS 长度
初始化：初始 f[1][0] = f[0][1] = f[0][0] = 0
转移方程：当 Ai == Bj 时，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 <cmath>
#include <cstring>
#include <cstdio>
using namespace std;
#define _for(i,a,b) for(int i = (a);i < (b); i++)
#define _rep(i,a,b) for(int i = (a);i <= (b); i++) 

const int N = 5100;

int n,m,mx;
int A[N],B[N];
int f1[N],f2[N][N];

init(){
	cin >> n >> m;
	_rep(i,1,n){
		cin >> A[i];
		f1[i] = 1;
	}
	_rep(i,1,m){
		cin >> B[i];
	}
	mx = 1;
}

void LIS(){
	_rep(i,2,n){
		_for(j,1,i){
			if(A[i] > A[j]){
//				cout << "f1[i]: " << f1[i] << endl;
				f1[i] = max(f1[j]+1,f1[i]);
			}
//			cout << "max: " << mx << endl;
			mx = max(mx,f1[i]);
		}
	}
}

void LCS(){
	memset(f2,0,sizeof(f2));
	_rep(i,1,n){
		_rep(j,1,m){
			if(A[i] == B[j])
				f2[i][j] = f2[i-1][j-1] + 1;
			else
				f2[i][j] = max(f2[i][j-1],f2[i-1][j]);
		}
	}
} 
int main(){
	init();
	LIS();
	LCS();
	cout << mx << " " << f2[n][m];
}