[ DP CDQ分治 ] BZOJ4553

设 $l_i$ 表示 $a_i$ 可以变成的最小值， $r_i$ 表示最大值， $f_i$ 表示以 $i$ 为结尾的最长子序列长度。

$$f_i=\max\{f_j+1~|~j<i,r_j\le a_i,a_j\le l_i\}$$
然后cdq分治就好了。

#include<bits/stdc++.h>
using namespace std;
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void Read(int& x){
    char c=nc();
    for(;c<'0'||c>'9';c=nc());
    for(x=0;c>='0'&&c<='9';x=(x<<3)+(x<<1)+c-48,c=nc());
}
const int N=100010;
struct Node{
    int x,y,id,z;
    bool b;
    Node(int x=0,int y=0,int z=0,int id=0,bool b=0):x(x),y(y),z(z),id(id),b(b){}
}b[N<<1],t[N<<1];
int k,n,m,Ans,mx,x,y;
int f[N];
int a[N],L[N],R[N];
int cnt,v[N],c[N];
inline void Update(int x,int y){
    for(;x<=mx;x+=x&-x)
    if(v[x]!=cnt)v[x]=cnt,c[x]=y;else c[x]=max(c[x],y);
}
inline int Query(int x){
    int Ans=0;
    for(;x;x-=x&-x)
    if(v[x]==cnt)Ans=max(Ans,c[x]);
    return Ans;
}
void Solve(int l,int r){
    if(l==r)return;
    int Mid=l+r>>1;
    int p1=l,p2=Mid+1;
    for(int i=l;i<=r;i++)if(b[i].z<=Mid)t[p1++]=b[i];else t[p2++]=b[i];
    for(int i=l;i<=r;i++)b[i]=t[i];
    Solve(l,Mid);
    cnt++;
    int j=l;
    for(int i=Mid+1;i<=r;i++){
        for(;j<=Mid&&b[j].x<=b[i].x;j++)
        if(!b[j].b)Update(b[j].y,f[b[j].id]);
        if(b[i].b)f[b[i].id]=max(f[b[i].id],Query(b[i].y)+1);
    }
    Solve(Mid+1,r);
    int p=l;p1=l;p2=Mid+1;
    while(p1<=Mid||p2<=r)
    if((p2>r)||(p1<=Mid&&b[p1].x<b[p2].x))t[p++]=b[p1++];else t[p++]=b[p2++];
    for(int i=l;i<=r;i++)b[i]=t[i];
}
inline bool Cmp(Node a,Node b){
    return a.x<b.x;
}
int main(){
    Read(n);Read(m);
    for(int i=1;i<=n;i++)Read(a[i]),L[i]=R[i]=a[i],f[i]=1;
    while(m--)Read(x),Read(y),L[x]=min(L[x],y),R[x]=max(R[x],y);
    for(int i=1;i<=n;i++)mx=max(mx,R[i]);
    for(int i=1;i<=n;i++){
        m++;b[m]=Node(a[i],L[i],m,i,1);
        m++;b[m]=Node(R[i],a[i],m,i,0);
    }
    sort(b+1,b+m+1,Cmp);
    Solve(1,m);
    for(int i=1;i<=n;i++)Ans=max(Ans,f[i]);
    cout<<Ans<<endl;
    return 0;
}