问题描述
(注: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的长度
Sample
Simple Input
5 5
1 3 2 5 4
2 4 3 1 5
Simple Output
3 2
思路
最长上升子序列是一个序列中最长的不连续的递增子序列,考虑二分优化,因为对于第i个状态来说,最优的子序列一定在每一长度中末尾元素最小的子序列中产生,这样对于所有子序列,长度增加,则末尾元素的高度增加,满足单调,可以二分出最长的末尾元素小于第i个元素的子序列。
最长公共子序列是指A和B中可以匹配到的最长不连续子序列,一维状态不好求解,则考虑二维状态,dp[n][m]代表A中前n个字符和B中前m个字符能匹配到的最大长度,可以转移的子状态有dp[n-1][m],dp[n][m-1],dp[n-1][m-1],最好的情况是若a[n]==b[m]则直接在dp[n-1][m-1]的基础上+1,否则在dp[n-1][m],dp[n][m-1]之中取最大值。
代码
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iomanip>
using namespace std;
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define mem(a,s) memset(a,s,sizeof(a))
//LIS:严格上升
int a[5010];
int b[5010];
int lis[5010];
int lcs[5010][5010];
int len=0;
int n, m;
void make_lis(){寻找第一个大于等于x的数,替换,没有则新建元素
rep(i,1,n){
int l=1,r=len,ans=len+1;
while(l<=r){
int m = (l + r)/2;
if(lis[m]>=a[i]){
ans = m;
r = m - 1;
}else
l = m + 1;
}
lis[ans] = a[i];
if(ans>len)
len = ans;
}
}
void make_lcs(){
rep(i,1,n)
rep(j,1,m)
if(a[i]==b[j])
lcs[i][j]=lcs[i-1][j-1]+1;
else
lcs[i][j] = max(lcs[i][j - 1], lcs[i - 1][j]);
}
int main(){
// freopen("in.txt","r",stdin);
scanf("%d%d", &n,&m);
rep(i, 1, n)
scanf("%d", &a[i]);
rep(i, 1, m)
scanf("%d", &b[i]);
make_lis();
make_lcs();
cout<<len <<' '<< lcs[n][m] << endl;
return 0;
}