背景

　　我是DP超級大蒟蒻，這種問題上課聽一臉懵逼，

　　於是今天下午拿出來好好補。。。

　　這裏只舉最長上升子序列的栗子

O(n^2)做法

　　在序列a　

　　設dp[i]是以第i個數結束的最長上升子序列的長度

　　那麽對於每一個i，枚舉j滿足j<i&&a[j]<a[i]的最大dp[j]

　　那麽dp[i]=dp[j]+1;

　　狀態轉移方程：dp[i]=max{dp[j]+1};(j<i&&a[j]<a[i])

　　代碼：

#include<cstdio>
#include<algorithm>
using
 
 namespace std;
const int MAXN=10010;

int n,a[MAXN],dp[MAXN],ans=0;

int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    dp[1]=1;
    for(int i=1;i<=n;i++){
        dp[i]=1;
        for(int j=1;j<i;j++){
            if(a[j]<a[i]&&dp[j]+1
 
>dp[i])
                dp[i]=dp[j]+1;
        }
    }
    for(int i=1;i<=n;i++){
        ans=max(ans,dp[i]);
    }
    printf("%d",ans);
    return 0;
} 

O(n*log2n)做法

　　用一個stack存儲當前的最長上升子序列，沒加上一個數qwq，判斷：

　　1.qwq數大於棧頂：直接壓入棧頂

　　2.否則：找到棧中最值得替換為qwq的數（即第一個大於等於qwq的數，可用stl的lower_bound求），將其替換為qwq。

例題

　　洛谷1020

　　第一問求最長不降子序列，第二問求最長上升子序列（因為打斷攔截的就是這個序列裏的元素）

　　代碼:

#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN=100010;

int n=1,a[MAXN],ans1=1,ans2=1;
int stack1[MAXN],stack2[MAXN];

struct cmp{
    bool operator()(int a,int b){
        return a>b;
    }
};

int main(){
    while(cin>>a[n])n++;
    n--;
    stack1[1]=stack2[1]=a[1];
    for(int i=2;i<=n;i++){
        if(a[i]<=stack1[ans1])
            stack1[++ans1]=a[i];
        else
            stack1[upper_bound(stack1+1,stack1+ans1+1,a[i],cmp())-stack1]=a[i];
        if(a[i]>stack2[ans2])
            stack2[++ans2]=a[i];
        else
            stack2[lower_bound(stack2+1,stack2+ans2+1,a[i])-stack2]=a[i];
    }
    printf("%d %d",ans1,ans2);
    return 0;
}

[OI學習筆記]最長上升子序列