【JZOJ 4606】 序列

Description

Analysis

记mx[i]为i位置上变化的MAX，mn[i]同理。
DP方程显然：

$$f[i]=Max_{j<i,mx[j]<a[i],a[j]<mn[i]}f[j]+1$$
把后两个约束化成二维平面上的，然后树状数组套线段树。
注意动态开节点。
时间 $O(nlog_2^2n)$，空间同时间，但实际情况远没有达到最坏值。

Code

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int N=100010,M=N*100,INF=2147483647;
int n,m,mxa,f[N],a[N],mx[N],mn[N],root[N],tr[M],lch[M],rch[M];
char ch;
void read(int &n)
{
    int t=0,p=1;
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar())
        if(ch=='-') p=-1;
    for(;'0'<=ch && ch<='9';ch=getchar()) t=t*10+ch-'0';
    n=t*p;
}
void change(int &v,int l,int r,int x,int y)
{
    if(!v) v=++m;
    tr[v]=max(tr[v],y);
    if(l==r) return;
    int mid=(l+r)>>1;
    if(x<=mid) change(lch[v],l,mid,x,y);
    else change(rch[v],mid+1,r,x,y);
    tr[v]=max(tr[lch[v]],tr[rch[v]]);
}
int query(int v,int l,int r,int x,int y)
{
    if(l==x && r==y) return tr[v];
    int mid=(l+r)>>1;
    if(y<=mid) return query(lch[v],l,mid,x,y);
    else
    if(x>mid) return query(rch[v],mid+1,r,x,y);
    else
    return max(query(lch[v],l,mid,x,mid),query(rch[v],mid+1,r,mid+1,y));
}
int lowbit(int x){return x&-x;}
void add(int x,int y,int z)
{
    for(;x<=mxa;x+=lowbit(x))
    change(root[x],1,mxa,y,z);
}
int getmax(int x,int y)
{
    int t=0;
    for(;x;x-=lowbit(x)) t=max(t,query(root[x],1,mxa,1,y));
    return t;
}
int main()
{
    int _,x,y,ans=0;
    read(n),read(_);
    fo(i,1,n) read(a[i]),mx[i]=mn[i]=a[i],mxa=max(mxa,a[i]);
    fo(i,1,_)
    {
        read(x),read(y);
        mx[x]=max(mx[x],y);
        mn[x]=min(mn[x],y);
    }
    fo(i,1,n)
    {
        f[i]=getmax(a[i],mn[i])+1;
        add(mx[i],a[i],f[i]);
        ans=max(ans,f[i]);
    }
    printf("%d",ans);
    return 0;
}