最长上升子序列(LIS)

最长上升子序列(LIS)

---

参考博客

数组中最长的单调递增子序列，$LIS$ 有两种做法：①时间复杂度为 $O(n^2)$ 的朴素 $DP$ ，②时间复杂度为 $O(nlogn)$ 的二分写法。

做法一：

设 $d{p}_i$ 表示以 $a_i$ 为结尾的单调递增子序列的长度，状态转移方程为 $d{p}_i={\max }_{j<i,a[j]<a[i]}d{p}_j+1$

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

void DP(){
	for(int i=1;i<=n;i++) {
        dp[i] = 1;
        for(int j=i-1;j>=1;j++) {
            if(a[i] > a[j]) dp[i] = max(dp[i],dp[j]+1);
        }
    }
}

做法二：

设 $d{p}_i$ 为长度为 $i$ 的上升子序列的最后一个元素的最小值，在维护一个 $len$ 变量，表示 $dp$ 数组的长度。

状态转移方程为：

• 如果有 $a>d{p}_{len}$ ，那么把元素 $a$ 加到 $dp$ 的后一位，上升的子序列长度加 $1$， $d{p}_{len+1}=a$
• 如果 $a<d{p}_{len}$ ，则找到 $dp$ 数组中第一个大于等于 $a$ 的值进行替换。替换后数组 $dp$ 仍然为上升子序列，且使得同一长度的序列的最后一个元素更小。

例如：序列 $1001101201234$

int a[N],dp[N],n,len;

void DP(){
	for(int i=1;i<=n;i++) {
        if(dp[len] < a[i]) {
            dp[++ len] = a[i];
        }
        else {
            int k = lower_bound(dp+1,dp+1+len,a[i]) - dp;
            dp[k] = a[i];
        }
    }
}

题目：P1020 NOIP1999 普及组 导弹拦截