【DP】LIS&LCS

问题描述

（注：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;
}