最长递增子序列

最长递增子序列

题目描述：

给定数组arr，设长度为n，输出arr的最长递增子序列。（如果有多个答案，请输出其中字典序最小的）

输入描述：

输出两行，第一行包括一个正整数n（n<=100000），代表数组长度。第二行包括n个整数，代表数组$arr\left(1\le ar{r}_i\le 1e9\right)$ 。

输出描述：

输出一行。代表你求出的最长的递增子序列。

示例1

输入

9
2 1 5 3 6 4 8 9 7

输出

1 3 4 8 9

示例2

输入

5
1 2 8 6 4

输出

1 2 4

说明

其最长递增子序列有3个，（1，2，8）、（1，2，6）、（1，2，4）其中第三个字典序最小，故答案为（1，2，4）

备注：

时间复杂度 $O(n\log n)$ ，空间复杂度 $O(n)$ 。

---

题解：

动态转移方程：dp[i] = max{dp[j] + 1 (0<=j<i, a[j]<a[i])}。dp[i] 表示以 a[i] 结尾的情况下，a[0…i] 中的最大递增子序列长度。普通的解法时间复杂度为 $O(n^2)$ ，因为我们需要遍历 i 前面的元素寻找最大的 dp[j] 。

那么有没有优化的方法呢？

考虑一下情况：**对于一个上升子序列，显然其结尾元素越小，越有利于在后面接其它的元素，也就越可能变得更长。**因为，我们可以维护一个 low 数组。low[i] 的含义是：所有长度为 i 的最长上升子序列中，末尾最小的数字。那么 low 一定是单调的，于是我们可以这样更新 low 数组：

• 若 a[i] > low[len] ，则 low[++len] = a[i]；
• 否则的话，二分搜索 low 中大于等于 a[i] 的位置 p ，low[p] = a[i]。

这样，最大的长度就是 low 数组的长度，并且还有一个好处：用更小的元素替换，可以帮助我们找到字典序最小的子序列。

注意：low 数组中最终存储的不一定是一个合法的最长上升子序列，它是帮助我们计算 dp[i] 的。在输出时，我们需要从 最大长度 位置往前迭代，如果 a[j] < a[i] && dp[j] == dp[i] - 1，则 a[j] 是 a[i] 的前一个元素。

代码：

#include <cstdio>
#include <vector>
#include <algorithm>

using namespace std;

const int N = 100010;

int n;
int a[N];
int f[N];
int low[N];
int ret[N];

int main(void) {
    scanf("%d", &n);
    scanf("%d", a + 0);
    int len = 0, idx = 0;
    f[0] = 0, low[0] = a[0];
    for (int i = 1; i < n; ++i) {
        scanf("%d", a + i);
        if (a[i] > low[len]) {
            low[++len] = a[i];
            f[i] = len;
            idx = i; // 最大长度的位置
        } else {
            int p = lower_bound(low, low + len + 1, a[i]) - low;
            low[p] = a[i];
            f[i] = p;
            // 这里进行更新的原因是：
            // 若更新 low 数组中的元素后，以 a[i] 为尾元素的长度与最大长度相等，则说明在长度一样时，
            // 以 a[i] 结尾的子序列字典序更小，这也是我上面说的可以帮助我们找 最小字典序
            if (f[i] == len) idx = i;
        }
    }
    int t = len;
    ret[len--] = a[idx];
    for (int i = idx; i >= 0; --i) {
        if (a[i] < a[idx] && f[i] == f[idx] - 1) {
            ret[len--] = a[i];
            idx = i;
        }
    }
    for (int i = 0; i <= t; ++i) printf("%d%c", ret[i], " \n"[i == t]);
    return 0;
}