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

二、思路与算法

求LIS（最长上升子序列）时，主要使用两个for循环，一个用来遍历序列中的所有元素，一个用来遍历序列中这个元素之前的所有元素。以i为最大元素的子序列的最大长度，是所有以j为最大元素的子序列长度+1的最大值，两个循环之后，LIS就是所有最大长度的最大值。
而求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]中存储的就是LCS。

三、代码实现

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;

int a[5100]={0};
int b[5100]={0};
int f[5100]={0};
int n=0,m=0;
int ans1=1,ans2=1;
int DP[5100][5100];

int main(){
    scanf("%d %d",&n,&m);
    for(int i=0;i<n;i++){
        scanf("%d",&a[i]);
        f[i]=1;
    } 
    for(int i=0;i<m;i++){
    	scanf("%d",&b[i]);
	}
    //计算序列A的LIS 
    for(int i=0;i<n;i++){
      for(int j=0;j<i;j++){
            if(a[j]<a[i]){
            	f[i]=max(f[i],f[j]+1);
			}   
        }
        ans1=max(f[i],ans1);
    }
    //计算序列AB的LCS
	memset(DP, 0, sizeof(DP)); 
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(a[i-1]==b[j-1]){
				DP[i][j]=max(DP[i][j],DP[i-1][j-1]+1);
			}
			else{
				DP[i][j]=max(DP[i][j-1],DP[i-1][j]); 
			}		
		}
	}	
	ans2=DP[n][m];
    printf("%d %d",ans1,ans2);
	return 0;
}

四、经验与总结

1. 代码中算法的时间复杂度都为平方级别，可以使用二分法等算法，优化时间复杂度到log级别。