這個題n^2暴力顯然，然後考慮優化，每次找前面的最大值，有點像是三維偏序，樹套樹和cdq都能做，這裏用cdq，sort的cdq好像比較簡單。。。

題幹：

題目描述

佳媛姐姐過生日的時候，她的小夥伴從某寶上買了一個有趣的玩具送給他。玩具上有一個數列，數列中某些項的值可能會變化，但同一個時刻最多只有一個值發生變化。現在佳媛姐姐已經研究出了所有變化的可能性，她想請教你，能否選出一個子序列，使得在任意一種變化中，這個子序列都是不降的？請你告訴她這個子序列的最長長度即可 。

註意：每種變化最多只有一個值發生變化。在樣例輸入1中，所有的變化是：

1 2 3 
2 2 3 
1
 
 3 3 
1 1 3
1 2 4 

選擇子序列為原序列，即在任意一種變化中均為不降子序列在樣例輸入2中，所有的變化是:

3 3 3
3 2 3

選擇子序列為第一個元素和第三個元素，或者第二個元素和第三個元素，均可滿足要
輸入輸出格式
輸入格式：

輸入的第一行有兩個正整數n, m，分別表示序列的長度和變化的個數。接下來一行有n個數，表示這個數列原始的狀態。接下來m行，每行有2個數x, y，表示數列的第x項可以變化成y這個值。1 <= x <= n。

輸出格式：

輸出一個整數，表示對應的答案

代碼：

#include<iostream>
#include
 
<cstdio>
#include<cmath>
#include<ctime>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(register int i = a;i <= n;++i)
#define lv(i,a,n) for(register int i = a;i >= n;--i)
#define
 
 clean(a) memset(a,0,sizeof(a))
const int INF = 1 << 30;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = 0;
    while(c = getchar(), c < ‘0‘ || c > ‘9‘)
        if(c == ‘-‘) op = 1;
    x = c - ‘0‘;
    while(c = getchar(), c >= ‘0‘ && c <= ‘9‘)
        x = x * 10 + c - ‘0‘;
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < 0) putchar(‘-‘), x = -x;
    if(x >= 10) write(x / 10);
    putchar(‘0‘ + x % 10);
}
#define N 100010
int n,m;
int a[N];
int minc[N],maxc[N];
int id[N],dp[N],tree[N];
int lowbit(int x)
{
    return x & -x;
}
void insert(int pos,int x)
{
    for(;pos <= 100000;pos += lowbit(pos))
    tree[pos] = max(tree[pos],x);
}
int query(int pos)
{
    int ans = 0;
    for(;pos;pos -= lowbit(pos))
    ans = max(tree[pos],ans);
    return ans;
}
void erase(int pos)
{
    for(;pos <= 100000;pos += lowbit(pos))
    {
        tree[pos] = 0;
    }
}
int Max(int l,int r)
{
    int maxn = 0;
    duke(i,l,r)
    {
        maxn = max(maxn,dp[i]);
    }
    return maxn;
}
void cdq(int l,int r)
{
    if(l == r)
    {
        dp[l] = max(dp[l],1);
        return;
    }
    int mid = (l + r) >> 1;
    cdq(l,mid);
    static int oid[N];
    duke(i,l,r)
    oid[i] = id[i];
    sort(id + l,id + mid + 1,[](int x,int y){return maxc[x] < maxc[y];});
    sort(id + mid + 1,id + r + 1,[](int x,int y){return a[x] < a[y];});
    static vector <int> todel;
    int p1 = l;
    for(int p2 = mid + 1;p2 <= r;p2++)
    {
        while(p1 <= mid && maxc[id[p1]] <= a[id[p2]])
        {
            insert(a[id[p1]],dp[id[p1]]);
            todel.push_back(a[id[p1]]);
            ++p1;
        }
        int tans = query(minc[id[p2]]) + 1;
        dp[id[p2]] = max(dp[id[p2]],tans);
    }
    while(!todel.empty())
    {
        erase(todel.back());
        todel.pop_back();
    }
    duke(i,l,r)
    id[i] = oid[i];
    cdq(mid + 1,r);
}
int main()
{
    read(n);read(m);
    duke(i,1,n)
    {
        read(a[i]);
        minc[i] = maxc[i] = a[i];
        id[i] = i;
    }
    duke(i,1,m)
    {
        int pos,x;
        read(pos);read(x);
        minc[pos] = min(minc[pos],x);
        maxc[pos] = max(maxc[pos],x);
    }
    cdq(1,n);
    int ans = Max(1,n);
    printf("%d\n",ans);
    return 0;
}

P4093 [HEOI2016/TJOI2016]序列