六、动态规划 -- 最长上升子序列 II

题目描述

给定一个长度为 $N$ 的数列，求数值严格单调递增的子序列的长度最长是多少。

输入格式

第一行包含整数 $N$。

第二行包含 $N$ 个整数，表示完整序列。

输出格式

输出一个整数，表示最大长度。

数据范围

$1\le N\le 100000，$
$-10^9\le 数列中的数\le 10^9$

样例

输入样例

7
3 1 2 1 8 5 6

输出样例

4

样例解释

在以上的输入中，最长的子序列为 1 2 5 6，长度为 4

---

思路

本题的数组范围加大，所以不能使用最长上升子序列Ⅰ的做法。需要进行优化。

(贪心) $O(nlogn)$

1. 使用数组a来存储输入序列，开辟一个新的数组q来存储结果，数组q[r]表示：长度为r的序列长度结尾的最小值。最终的数组q的长度就是最长上升子序列的长度。
2. 假设数组q中存在数，那么当循环到数组a的第i个位置时，需要判断当前的数a[i]可以插在数组q的哪个位置，如果找不到插入的位置，那么就在最后把长度 + 1，并且把这个a[i]插入到数组q的尾端。
3. 使用二分来查找大于a[i]的最小的数，然后插入到找到的数字的后面，整个数组q是一个单调递增的数组。

---

---

举个栗子

坐标轴4表示的是序列长度为4时,序列结尾的最小值
例如 1 2 4 7 和1 2 4 9
这时候选哪个比较好?如果不清楚,那么插入8
第一个 1 2 4 7 8
第二个 1 2 4 9 8
明显第一个方案要更好一点,代码中用数组q[r]来表示长度为r的序列长度结尾的最小值
所以此时q[4]=7更妥当一些。

具体模拟过程如下：

n = 10
数组a 3 1 2 0 8 5 6 4 7 9
start:

// q为空， 第一个数字3加入到q
q: 3
len = 1

// 判断第二个数字
第一步：1比3要小, r = 0, q[r + 1] = q[1] = 1;
q: 1
len = 1

第二步：2比1大, r = 1,所以len = r + 1, q[r + 1] = q[2] = 2;
q: 1 2
len = 2

第三步：没有比0小的, 看0是不是可以更新数组q, 寻找下标：r = 0, q[r + 1] = q[1] = 0;
q:0 2

第四步：8比2大,加入到q数组，r = 2,len = max(len, r + 1) = 3, q[r + 1] = q[3] = 8
q:0 2 8

第五步：5比2大,5比8小,可以更新q数组, r = 2,len = max(len,r + 1) = 3,q[r + 1] = q[3] = 5;
q:0 2 5

第六步：6比5大,加入到q数组,r = 3,len = max(len, r + 1) = 4,q[r + 1] = q[4] = 6;
q: 0 2 5 6

第七步：4比2大,4比5小,寻找5的下标并更新数组q,r = 2,len = max(len, r + 1) = 4,q[r + 1] = q[3] = 4
q: 0 2 4 6

第八步：7比6大, r = 4, len = max(len, r + 1) = 5,q[r + 1] = q[5] = 7
q: 0 2 4 6 7

第九步：9比7大,r = 5,len = max(len, r + 1) = 6, q[r + 1] = q[6] = 9
q: 0 2 4 6 7 9

最终的最长上升子序列为6，序列为 1 2 5 6 7 9,以数字9为结尾。

---

C++代码

#include <iostream>

using namespace std;

const int N = 100010;

int a[N], q[N];
int n;

int main()
{
    scanf("%d", &n);
    
    for (int i = 0; i < n; i ++) scanf("%d", &a[i]);
    
    int len = 0;
    for (int i = 0; i < n; i ++)
    {
        int l = 0, r = len;
        while (l < r)
        {
            int mid = l + r + 1 >> 1;
            if (q[mid] < a[i]) l = mid;
            else r = mid - 1;
        }
        
        len = max(len, r + 1);
        q[r + 1] = a[i];
    }
    
    cout << len << endl;
    return 0;
}

参考链接

最长上升子序列