浅谈LIS

danxiangfa

已于 2023-04-01 10:21:17 修改

阅读量67

点赞数

文章标签：算法动态规划

于 2023-01-15 11:09:08 首次发布

本文链接：https://blog.csdn.net/danxiangfa/article/details/128692699

版权

LIS（最长上升子序列）

不难想出 $O(n^2)$ 的dp。

记 $f [i]$ 为以 $a [i]$ 结尾的LIS长度。

$f[i]=max^{i-1}_{j=1,a[i]>a[j]} f[j]+1$

显然，有 $f [1] = 1$ 。

$O(nlog_n)$ 的做法

显然对于一个LIS序列，他的结尾越小，就越优。我们维护一个 $f$ 数组，令 $f [i]$ 表示LIS长度为 $i$ 时，结尾数字的最小值（这个 $f$ 数组是单调递增的）。对于每一个 $a [i]$ ，如果大于 $f$ 数组末尾的数字，就把他加入 $f$ 数组的末尾，否则就用 $a [i]$ 替换掉 $f$ 数组中第一个大于等于 $a [i]$ 的数。

#include <bits/stdc++.h>
using namespace std;
int n, ans, a[100001], f[100001];
int main() {
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin >> n;
    for(int i = 1; i <= n; i++)
        cin >> a[i];
    ans = 1; f[1] = a[1];
    for(int i = 2; i <= n; i++)
        if(a[i] > f[ans])
            f[++ans] = a[i];
        else *lower_bound(f + 1, f + ans + 1, a[i]) = a[i];
    cout << ans << endl;
    return 0;
}

最长不降子序列

#include <bits/stdc++.h>
using namespace std;
int n, ans, a[100001], f[100001];
int main() {
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin >> n;
    for(int i = 1; i <= n; i++)
        cin >> a[i];
    ans = 1; f[1] = a[1];
    for(int i = 2; i <= n; i++)
        if(a[i] >= f[ans])
            f[++ans] = a[i];
        else *upper_bound(f + 1, f + ans + 1, a[i]) = a[i];
    cout << ans << endl;
    return 0;
}