【模型】LIS

最新推荐文章于 2024-09-06 17:52:44 发布

Nekroz_

最新推荐文章于 2024-09-06 17:52:44 发布

阅读量529

点赞数

分类专栏： DP LIS 文章标签： c++ DP LIS

本文链接：https://blog.csdn.net/Diogenes_/article/details/79927633

版权

DP 同时被 2 个专栏收录

12 篇文章 0 订阅

订阅专栏

LIS

1 篇文章 0 订阅

订阅专栏

本文介绍了最长上升子序列（LIS）问题，包括两种解决方案：O(n^2)的时间复杂度的朴素动态规划方法和O(n log n)的优化算法。文章通过代码示例详细阐述了这两种方法的实现过程，并指出虽然定义有所调整，但不影响最终结果。最终，通过求解f数组确定最长上升子序列的长度。

摘要由CSDN通过智能技术生成

LIS

最长上升子序列，会有各种变种，比如说最长不下降子序列，最长不上升子序列等等，代码差不多，这里以最长上升子序列为例。最长上升子序列，就是在一个序列A中找到一个一个序列中最长的单调递增的子序列（可以不连续）。
LIS有两种解决方法法， $O(n^2)$ & $O(nlog_2n)$ ，在这里都给出。

$\S$ LIS $O(n^2)$

朴素的DP，方程可以直接写出：
$f_i = max{f_j + 1}, 其中1\leq j \leq i, a[j] < a[i]$
$Ans = maxelement f$

$code:$

int LIS(std::vector< int > &a) {
    int ans = 1; 
    std::vector< int > f(a.size(), 1);
    for (unsigned int i = 2; i <= a.size() - 1; i++) 
        for (unsigned int j = i - 1; j >= 1; j--) 
            if (a[i] > a[j] && f[i] < f[j] + 1) {
                f[i] = f[j] + 1;
                if (f[i] > ans)
                    ans = f[i];
            }
    return ans;
}

$\S$ LIS $O(nlog_2n)$

我们对 $f$ 数组重新定义， $f_i$ 定义为长度为i的上升子序列的最后一个元素（本人在程序中使用了vector，f中所有元素都往前平移一位，所以定义发生改变，但不影响结果），初始化 $f[1] = a[1]$ ，然后每插入一个 $a[i]$ ,判断它和 $f$ 中最后一个元素的大小，如果大于，就插入，如果不大于的话，就从f中找到第一个大于 $a[i]$ 的，插在后面（可以用upper_bound），
$Ans = f.size()$

$code:$

int LIS(std::vector< int > &a) {
    std::vector< int > f;
    f.push_back(a[1]);
    for (unsigned int i = 2; i < a.size(); i++)
        if (a[i] < f.back()) 
            f[std::upper_bound(f.begin(), f.end(), a[i]) - f.begin()] = a[i]; 
        else f.push_back(a[i]); 
    return f.size();
}