最长上升子序列

简介

        最长上升子序列（Longest Increasing Subsequence，简称 LIS）是指在一个序列中找到一个最长的子序列，使得子序列中的元素按照从小到大的顺序排列，且子序列中的元素在原序列中的相对顺序保持不变。

算法

        朴素算法

         约定数组的长度为n,可以很方便地得到,对于任意 0<=j,i<n,如果a[i]<a[j],dp[i] = max(dp[i],dp[j]+1) 

         得到如下代码：       

#include<iostream>

using namespace std;
const int N = 1e5+10;
int dp[N];
int a[N];
int n;
#define rep(i,a,b) for(int i = a;i<b;i++)
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
int main(){
    IOS;
    cin>>n;
    rep(i,0,n)cin>>a[i];
    rep(i,0,n)
        rep(j,0,i)
            if(a[i]>a[j])dp[i] = max(dp[i],dp[j]+1);
    int maxn = 0;
    rep(i,0,n)
        maxn = max(maxn,dp[i]+1);
    cout<<maxn<<endl;
}

算法时间复杂度为O（n^2），可以看出，代码的状态数和单词转移都是O(n)，可以对于转移可以借助二分查找进行优化；优化后的算法时间复杂度为O（nlogn）

        优化

考虑维护一个数组dp及当前长度len，dp[i]表示长度为i的上升子序列的最后一个元素的最小值。 然后考虑如何转移：

        dp[len]<a[i] 就让显然可以直接把a[i]直接加入这个上升子序列的后面，len+1

        如果dp[len]>=a[i]就找到第一个大于a[i]的位置,由于只想要最长的子序列的长度，可以直接用a[i]替换，不会影响最终长度。

lower_bound()和upper_bound()

        介绍一下lower_bound()和upper_bound()函数，

        upper_bound 和 lower_bound 是 C++ STL 中的两个函数，用于在已排序的容器（如 vector、set、map 等）中查找元素的位置。

        lower_bound 函数返回的是第一个大于等于指定值的元素位置，而 upper_bound 函数返回的是第一个大于指定值的元素位置。

        例如，对于一个已排序的数组 {1, 2, 3, 4, 4, 5, 6}，lower_bound(4) 的结果是指向第一个 4 的位置，upper_bound(4) 的结果是指向第一个大于 4 的元素位置，即指向 5 的位置。

如果在给定的容器中不存在指定的值，lower_bound 和 upper_bound 函数将返回指向容器中第一个大于指定值的元素位置。

lower_bound只能用于不减序列，而最长下降子序列维护的是不增序列，可以借助greater<int>，更改成

*lower_bound(dp + 1, dp + len + 1, a[i], greater<int>()) = a[i];

 代码实现

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1e6+10;
int a[N],dp[N];
int len;
#define rep(i,a,b) for(int i = a;i<b;++i)
int main()
{
    int n;
    scanf("%d",&n);
    rep(i,0,n)scanf("%d",&a[i]);
    rep(i,0,n)
    {
        if(dp[len]<a[i])
            dp[++len] = a[i];
        else
            *lower_bound(dp+1,dp+len+1,a[i]) = a[i];
    }
    cout<<len<<endl;
}

B3637 最长上升子序列

P1020 [NOIP1999 普及组] 导弹拦截