Codeforces10 D. LCIS（最长公共上升子序列LCIS）

题意：

给定长度为n,m的序列
求两个序列的最长公共子序列
输出方案

数据范围：n,m<=505

解法：

令d[i][j]表示序列1前i个数,序列2前j个数的LCIS,
发现这样就无法确定最后一个数是什么了,无法进行转移,

d[i][j]表示前i个a[i]中，以b[j]为结尾的最大公共上升子序列
则对于一个固定的i和j,d[i][j]的值有两种情况：
1.不含a[i],则d[i][j]=d[i-1][j]
2.包含a[i],那么a[i]是结尾,又因为b[j]是结尾,因此需要先满足a[i]=b[j]
这时候d[i][j]=max(d[i-1][k])+1,(k<j,且b[j]>b[k]),


枚举i一层循环，枚举j一层循环，枚举k一层循环，复杂度O(n^3),
优化:
上面的情况2转移方程:d[i][j]=max(d[i-1][k])+1,(k<j,且b[j]>b[k]),
因为此时a[i]==b[j],所以b[j]>b[k]可以换成a[i]>b[k]
d[i-1][k]用一个变量(指针)存起来,不断维护就可以优化掉枚举k了
即固定i,而j从1到m枚举,维护一个ma,满足d[i-1][k]是d[i-1][1到j-1]中最大的

再看一遍转移方程:
d[i][j]=max(d[i-1][k])+1,(k<j,且b[j]>b[k]),替换b[j]为a[i],得:
d[i][j]=max(d[i-1][k])+1,(k<j,且a[i]>b[k]),
利用了i的不变,而j递增的性质,维护1-(j-1)中最优k

代码:
for(int i=1;i<=n;i++){
    int ma=0;//记录f[i-1][k]中最大的
    for(int j=1;j<=n;j++){
        f[i][j]=f[i-1][j];//不含a[i]
        if(a[i]==b[j]){//包含a[i],则结尾为a[i],又因为结尾为b[j],则要满足a[i]=b[j]
            f[i][j]=max(f[i][j],ma+1);//f[i-1][k]中最大的+1,(k<j)
        }
        if(a[i]>b[j]){//更新ma
            ma=max(ma,f[i-1][j]);
        }
    }
}

第一维可以降掉,因为初始化的时候d[i][j]=d[i-1][j]总是i继承i-1

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=505;
int pre[maxm];
int d[maxm];//d[j]表示以a[j]结尾的LCIS
int a[maxm];
int b[maxm];
int ans[maxm],num;
int n,m;
signed main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        scanf("%d",&b[i]);
    }
    for(int i=1;i<=n;i++){
        int ma=0;//记录满足d[i-1][k]最大的k,其中k<当前j
        for(int j=1;j<=m;j++){
            if(a[i]==b[j]){
                if(d[ma]+1>d[j]){
                    d[j]=d[ma]+1;
                    pre[j]=ma;
                }
            }
            if(a[i]>b[j]){//不包含a[i]
                if(d[j]>d[ma]){//维护1到j满足a[i]>b[j]的d[j]最大的j,用于d[j+1]更新
                    ma=j;
                }
            }
        }
    }
    int ma=0;
    for(int i=1;i<=m;i++){
        if(d[i]>d[ma])ma=i;
    }
    printf("%d\n",d[ma]);
    for(int i=ma;i;i=pre[i]){
        ans[++num]=b[i];
    }
    for(int i=num;i>=1;i--){
        printf("%d ",ans[i]);
    }
    return 0;
}

---