BZOJ 4553: [Tjoi2016&Heoi2016]序列

我就不懂了怎么又是线段树啊（好吧其实可以写平衡树）

f[i]可以从f[j](a[j]<=l[i]&&r[j]<=a[i])转移过来

观察条件发现就是某个矩阵里的权最大的点

权值线段树套权值线段树就可以n(logn)^2了，但是空间也是这么多，会MLE

所以三维偏序什么的还是cdq分治吧

然后就轻易地过掉了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<algorithm>
#include<map>
#include<set>
#define rep(i,l,r) for(int i=l;i<=r;i++)
#define per(i,r,l) for(int i=r;i>=l;i--)
#define mmt(a,v) memset(a,v,sizeof(a))
#define tra(i,u) for(int i=head[u];i;i=e[i].next)
using namespace std;
const int N=100000+5;
int l
,r
,a
;
struct Node{
int l,r,mx,ti;
}tr[N<<2];
#define lc o<<1
#define rc o<<1|1
int T;
void clr(){T++;}
void refresh(int o){
if(tr[o].ti!=T)tr[o].ti=T,tr[o].mx=0;
}
void pushup(int o){
tr[o].mx=max(tr[lc].mx,tr[rc].mx);
}
void update(int o,int p,int v){
refresh(o);
int l=tr[o].l,r=tr[o].r;
if(l==r)tr[o].mx=max(tr[o].mx,v);
else{
int mid=l+r>>1;
if(p<=mid)update(lc,p,v),refresh(rc);
else update(rc,p,v),refresh(lc);
pushup(o);
}
}
int query(int o,int a,int b){
int l=tr[o].l,r=tr[o].r;
if(l==a&&b==r){refresh(o);return tr[o].mx;}
else{
int mid=l+r>>1;
if(b<=mid)return query(lc,a,b);
else if(mid<a)return query(rc,a,b);
else return max(query(lc,a,mid),query(rc,mid+1,b));
}
}
void build(int o,int l,int r){
tr[o].l=l;tr[o].r=r;tr[o].mx=0;
if(l==r)return;
int mid=l+r>>1;
build(lc,l,mid);build(rc,mid+1,r);
}
bool cmp1(int i,int j){return a[i]<a[j];}
bool cmp2(int i,int j){return l[i]<l[j];}
int p
,ans
;
void cdq(int x,int y){
if(x==y)return;
int mid=x+y>>1;
cdq(x,mid);
sort(p+x,p+mid+1,cmp1);
sort(p+mid+1,p+y+1,cmp2);
clr();
int j=mid+1;
rep(i,x,mid){
while(j<=y&&l[p[j]]<a[p[i]])
ans[p[j]]=max(ans[p[j]],query(1,1,a[p[j]])+1),j++;
update(1,r[p[i]],ans[p[i]]);
}
while(j<=y)ans[p[j]]=max(ans[p[j]],query(1,1,a[p[j]])+1),j++;
sort(p+mid+1,p+y+1);
cdq(mid+1,y);
}
int main(){
//freopen("a.in","r",stdin);
int n,m;scanf("%d%d",&n,&m);
rep(i,1,n)scanf("%d",&a[i]),l[i]=r[i]=a[i],p[i]=i;
while(m--){
int x,y;scanf("%d%d",&x,&y);
l[x]=min(l[x],y);
r[x]=max(r[x],y);
}
build(1,1,n);
rep(i,1,n)ans[i]=1;
cdq(1,n);
int ret=0;
rep(i,1,n)ret=max(ret,ans[i]);
printf("%d\n",ret);
return 0;
}