[CH弱省胡策R2]TATT

一、题目

点此看题

二、解法

先排序，把$a$当作第一关键字，$b$当第二关键字，$c$当第三关键字，$d$当第四关键字排序，然后设$f[i]$前$i$个组成的最大经过点数量（必须包含$i$），转移要找到最大$f[j]$，使得$a_j\le a_i,b_j\le b_i.....$，第一维通过排序解决了，剩下的三个我们可以建一个$\text{k-d}$树直接查找，知道$f[i]$之后插入$\text{k-d}$树即可。

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int M = 100005;
const double al = 0.7;
const int inf = 1e9;
int read()
{
    int x=0,flag=1;char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,m,k,opt,res,top,ans,rt,cnt,siz[M],b[M];
int mi[M][3],mx[M][3],Max[M],ls[M],rs[M];
struct node
{
    int x[3],val;
}a[M],v[M];
struct data
{
    int x[4];
    bool operator < (const data &b) const
    {
        for(int i=0;i<4;i++)
            if(x[i]^b.x[i])
                return x[i]<b.x[i];
        return 0;
    }
}s[M];
int newnode()
{
    if(top) return b[top--];
    return ++cnt;
}
void up(int x)
{
    for(int i=0;i<3;i++)
        mi[x][i]=mx[x][i]=v[x].x[i];
    if(ls[x])
        for(int i=0;i<3;i++)
        {
            mx[x][i]=max(mx[x][i],mx[ls[x]][i]);
            mi[x][i]=min(mi[x][i],mi[ls[x]][i]);
        }
    if(rs[x])
        for(int i=0;i<3;i++)
        {
            mx[x][i]=max(mx[x][i],mx[rs[x]][i]);
            mi[x][i]=min(mi[x][i],mi[rs[x]][i]);
        }
    siz[x]=siz[ls[x]]+siz[rs[x]]+1;
    Max[x]=v[x].val;
    if(ls[x]) Max[x]=max(Max[x],Max[ls[x]]);
    if(rs[x]) Max[x]=max(Max[x],Max[rs[x]]);
}
int cmp(node a,node b)
{
    return a.x[opt]<b.x[opt];
}
void build(int &x,int l,int r,int wd)
{
    if(l>r) return ;
    x=newnode();
    int mid=(l+r)>>1;
    opt=wd;
    nth_element(a+l,a+mid,a+r+1,cmp);
    v[x]=a[mid];
    build(ls[x],l,mid-1,(wd+1)%3);
    build(rs[x],mid+1,r,(wd+1)%3);
    up(x);
}
void dfs(int &x)
{
    if(!x) return ;
    dfs(ls[x]);
    b[++top]=x;
    a[++n]=v[x];
    dfs(rs[x]);
    x=0;
}
void check(int &x,int wd)
{
    if(siz[ls[x]]>al*siz[x] || siz[rs[x]]>al*siz[x])
    {
        n=0;dfs(x);
        build(x,1,n,wd);
    }
}
void ins(int &x,int wd,node y)
{
    if(!x)
    {
        x=newnode();v[x]=y;up(x);
        return ;
    }
    if(v[x].x[wd]>=y.x[wd]) ins(ls[x],(wd+1)%3,y);
    else ins(rs[x],(wd+1)%3,y);
    up(x);check(x,wd);
}
int in(int x1,int x2,int y1,int y2,int z1,int z2,int X1,int X2,int Y1,int Y2,int Z1,int Z2)
{
    return x1<=X1 && X2<=x2 && y1<=Y1 && Y2<=y2 && z1<=Z1 && Z2<=z2;
}
int out(int x1,int x2,int y1,int y2,int z1,int z2,int X1,int X2,int Y1,int Y2,int Z1,int Z2)
{
    return x1>X2 || X1>x2 || y1>Y2 || Y1>y2 || z1>Z2 || Z1>z2;
}
void ask(int x,int a,int b,int c)
{
    if(!x || ans>=Max[x]) return ;
    if(in(-inf,a,-inf,b,-inf,c,mi[x][0],mx[x][0],mi[x][1],mx[x][1],mi[x][2],mx[x][2]))
        ans=max(ans,Max[x]);
    if(out(-inf,a,-inf,b,-inf,c,mi[x][0],mx[x][0],mi[x][1],mx[x][1],mi[x][2],mx[x][2]))
        return ;
    if(in(-inf,a,-inf,b,-inf,c,v[x].x[0],v[x].x[0],v[x].x[1],v[x].x[1],v[x].x[2],v[x].x[2]))
        ans=max(ans,v[x].val);
    ask(ls[x],a,b,c);ask(rs[x],a,b,c);
}
signed main()
{
    m=read();
    for(int i=1;i<=m;i++)
        for(int j=0;j<4;j++)
            s[i].x[j]=read();
    sort(s+1,s+1+m);
    for(int i=1;i<=m;i++)
    {
        ans=0;ask(rt,s[i].x[1],s[i].x[2],s[i].x[3]);
        ans++;
        ins(rt,0,node{s[i].x[1],s[i].x[2],s[i].x[3],ans});
        res=max(ans,res);
    }
    printf("%d\n",res);
}