【题解】[CH弱省胡策R2]TATT

本蒟蒻第一道\(K-D-Tree\)维护\(dp\)

Question

题目大意：求一条路径，使得其四个维度单调不降。

先排序消掉一维再说。

对于每一个点，初始的时候绝对长度是1啊。于是，先赋值一个1，对于每一个点。

设计\(dp\)数组

\[f[i]=max_{f[j]}(a[j]<=a[i],b[j]<=b[i],c[j]<=c[i],d[j]<=d[i]) \]

那问题就转为，对于每一个点，如何求出在它之前的最大\(f[i]\)值。

发现问题类似于三维偏序，正好\(K-D-Tree\)硬刚即可：

对于每一个点，维护这个点的\(f[i]\)值，以及这个点的四个坐标；而树上只需要维护三个维度的上下界以及区间\(MAX\),动态\(query\)即可。

（有点懵的是，不能提前处理排序后的第一个点，而必须对所有点一次插入）

\(hole:\)对于\(sort\)函数的\(cmp\)函数，应该以四个维度分别为第\(1,2,3,4\)关键字排顺序，而不是只排一个。

\(Code:\)

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=5e5+10;
inline int read(){
	int s=0,w=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-')w=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		s=(s<<1)+(s<<3)+(ch^48);
		ch=getchar();
	}
	return s*w;
}
struct pt{
	int x[4],cnt;
}p[MAXN];
#define Alpha 0.75
#define inf 2147483647
struct node{
	int mx[3],mi[3],siz,maxn;
	pt c;
}tr[MAXN];
int n,m,D,rt,tot,top,ans=-inf;
int ls[MAXN],rs[MAXN],rub[MAXN];
int operator<(pt a,pt b){return a.x[D]<b.x[D];}
inline void pushup(int x){
	int l=ls[x],r=rs[x];
	tr[x].siz=tr[l].siz+tr[r].siz+1;
	for(int i=0;i<=2;++i){
		tr[x].mi[i]=tr[x].mx[i]=tr[x].c.x[i];
		tr[x].maxn=tr[x].c.cnt;
		if(l)tr[x].mi[i]=min(tr[x].mi[i],tr[l].mi[i]),
			 tr[x].mx[i]=max(tr[x].mx[i],tr[l].mx[i]),
			 tr[x].maxn=max(tr[x].maxn,tr[l].maxn);
		if(r)tr[x].mi[i]=min(tr[x].mi[i],tr[r].mi[i]),
			 tr[x].mx[i]=max(tr[x].mx[i],tr[r].mx[i]),
			 tr[x].maxn=max(tr[x].maxn,tr[r].maxn);
	}
}
inline bool cmp(pt a,pt b){
	return a.x[3]<b.x[3]||(a.x[3]==b.x[3]&&a.x[2]<b.x[2])||(a.x[3]==b.x[3]&&a.x[2]==b.x[2]&&a.x[1]<b.x[1])||(a.x[3]==b.x[3]&&a.x[2]==b.x[2]&&a.x[1]==b.x[1]&&a.x[0]<b.x[0]);
}
inline int New(){
	if(top)return rub[top--];
	else return ++tot;
}
int build(int l,int r,int d){
	if(l>r)return 0;
	int x=New(),mid=l+r>>1;
	D=d;nth_element(p+l,p+mid,p+r+1);
	tr[x].c=p[mid];ls[x]=build(l,mid-1,(d+1)%3);
	rs[x]=build(mid+1,r,(d+1)%3);pushup(x);return x;
}
void clear(int x,int pos){
	if(ls[x])clear(ls[x],pos);
	p[pos+tr[ls[x]].siz+1]=tr[x].c;rub[++top]=x;
	if(rs[x])clear(rs[x],pos+tr[ls[x]].siz+1);
}
void check(int &x,int d){
	double C=Alpha*(double)(tr[x].siz);
	if(C<(double)tr[ls[x]].siz||C<(double)tr[rs[x]].siz){clear(x,0);x=build(1,tr[x].siz,d);}
}
void Ins(int &x,pt s,int d){
	if(!x){x=New();ls[x]=rs[x]=0;tr[x].c=s;pushup(x);return;}
	if(s.x[d]<=tr[x].c.x[d])Ins(ls[x],s,(d+1)%3);
	else Ins(rs[x],s,(d+1)%3);
	pushup(x);check(x,d);
}
int f[MAXN];
bool Check(int x,pt s){
	int fg=0;
	for(int i=0;i<=2;++i){
		if(tr[x].mx[i]>s.x[i]){
			fg=1;break;
		}
	}
	if(fg)return false;
	else return true;
} 
bool CK(int x,pt s){
	int fg=0;
	for(int i=0;i<=2;++i){
		if(tr[x].mi[i]>s.x[i])return true;
	}
	return false;
}
inline bool ck(pt a,pt b){
	int fg=1;
	for(int i=0;i<=2;++i){
		if(a.x[i]>b.x[i]){
			fg=0;
			break;
		}
	}
	return fg;
} 
int query(int x,pt a){
	int res=-inf;
	if(CK(x,a))return 0;
	else if(Check(x,a))return tr[x].maxn;
	else {
		if(ck(tr[x].c,a))res=max(res,tr[x].c.cnt);
		res=max(query(ls[x],a),res);res=max(query(rs[x],a),res);
		return res;
	}
}
int main(){
	n=read();
	for(int i=1;i<=n;++i){
		p[i].x[0]=read(),p[i].x[1]=read(),p[i].x[2]=read(),p[i].x[3]=read();
	}
	for(int i=1;i<=n;++i)f[i]=1;
	sort(p+1,p+n+1,cmp);
	for(int i=1;i<=n;++i){
		f[i]+=query(rt,p[i]);
	//	cout<<"当前点f[i]:"<<i<<"<- ->"<<f[i]<<endl;
		p[i].cnt=f[i];
		Ins(rt,p[i],0);
	}
	for(int i=1;i<=n;++i)ans=max(ans,f[i]);
	printf("%d\n",ans);
	return 0;
}