CF10D LCIS

CF10D LCIS - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

是远古的 CF 题一枚呀。

不难想到状态，$f(i, j)$ 表示 $A[1, i]$ 和 $B[1, j]$ 内可构成的 LCIS 长度。

然后会发现，保证上升的这一部分转移是比较难以进行的，因为我们想得到单调性，却没有一个明确的比较目标。

回顾 LIS 的做法，$f(i)$ 表示 $A[1, i]$ 且以 $A[i]$ 结尾的 LIS 长度。

那么修改状态，$f(i, j)$ 表示 $A[1, i]$ 和 $B[1, j]$ 内，且以 $B[j]$ 结尾的 LCIS 长度。

这样转移方程就很好写了，我们令 $A[0] = B[0] = -\infty$：

$$f(i, j) =\begin{cases}f(i - 1, j) & A[i] \ne B[j] \\\max\limits_{0\le k < j, B[k] < B[j]}\{f(i - 1, k)\} + 1 & A[i] = B[j]\end{cases}$$

边界：$f(i, 0) = f(0, j) = 0$。

容易看出，直接这么写是 $\mathcal{O}(nm^2)$ 的。可能已经足以通过本题？不过事实上还有一种 $\mathcal{O}(nm)$ 的做法。

发现瓶颈在于 $k$ 的循环枚举，事实上这是可以优化的。我们发现在上面转移方程的第二种情况中，$B[k] < B[j]$ 等价于 $B[k] < A[i]$，也就是这个条件的是否成立和 $j$ 的取值无关。我们记满足 $0 \le k < j, B[k] < A[i]$ 的所有 $k$ 构成的集合记为 $S(i, j)$。发现有：

$$S(i, j +1) =\begin{cases}S(i, j) & B[j] \ge A[i]\\S(i, j) \cup \{j\} & B[j] < A[i]\end{cases}$$

因此枚举 $j$ 时同时维护 $S$ 中所有 $k$ 对应的 $f(i - 1, k)$ 的最大值即可。具体来说，设 $g(i, j)$ 为 $\max \limits _{k\in S(i, j)}\{f(i-1, k)\}$，则有：

$$g(i, j +1) = \begin{cases}g(i, j) & B[j] \ge A[i] \or f(i - 1, j) \le g(i, j)\\ f(i - 1, j) & B[j] < A[i] \and f(i - 1, j) > g(i, j)\end{cases}$$

用一个变量滚动地维护 $g$ 即可。

记录方案直接 dp 时记录前驱。

/*
 * @Author: crab-in-the-northeast 
 * @Date: 2022-11-06 00:59:09 
 * @Last Modified by: crab-in-the-northeast
 * @Last Modified time: 2022-11-06 03:02:53
 */
#include <bits/stdc++.h>
inline int read() {
    int x = 0;
    bool flag = true;
    char ch = getchar();
    while (!isdigit(ch)) {
        if (ch == '-')
            flag = false;
        ch = getchar();
    }
    while (isdigit(ch)) {
        x = (x << 1) + (x << 3) + ch - '0';
        ch = getchar();
    }
    if(flag)
        return x;
    return ~(x - 1);
}

const int maxn = 505;
const int maxm = 505;

int a[maxn], b[maxm];
int f[maxn][maxm];
int pre[maxn][maxm];

void pr(int x, int n) {
    if (x == 0)
        return ;
    pr(pre[n][x], n);
    printf("%d ", b[x]);
}

int main() {
    int n = read();
    for (int i = 1; i <= n; ++i)
        a[i] = read();
    int m = read();
    for (int i = 1; i <= m; ++i)
        b[i] = read();
    
    a[0] = b[0] = -1;

    for (int i = 1; i <= n; ++i) {
        int k = 0;
        for (int j = 1; j <= m; ++j) {
            if (a[i] != b[j]) {
                f[i][j] = f[i - 1][j];
                pre[i][j] = pre[i - 1][j];
            } else {
                f[i][j] = f[i - 1][k] + 1;
                pre[i][j] = k;
            }

            if (b[j] < a[i] && f[i - 1][j] > f[i - 1][k])
                k = j;
        }
    }

    int st = 0;
    for (int i = 1; i <= m; ++i)
        if (f[n][i] > f[n][st])
            st = i;

    printf("%d\n", f[n][st]);
    pr(st, n);
    puts("");
    return 0;
}