[CH弱省胡策R2]TATT

description

洛谷

data range

\[ n\le 5\times 10^4\]

solution

这就是四维偏序了...

好象时间复杂度是\(O(n^{\frac{5}{3}})\)...

因为可以剪枝所以速度快了一些...

code

#include<bits/stdc++.h>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cstring>
#include<complex>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<ctime>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define FILE "a"
#define Cpy(x,y) (memcpy(x,y,sizeof(x)))
#define Set(x,y) (memset(x,y,sizeof(x)))
#define mp make_pair
#define pb push_back
#define RG register
#define il inline
using namespace std;
typedef unsigned long long ull;
typedef pair<int,int> PI;
typedef vector<int>VI;
typedef long long ll;
typedef double dd;
const dd eps=1e-10;
const int mod=1e9+7;
const int N=200010;
const int M=205;
const dd pi=acos(-1);
const int inf=2147483647;
const ll INF=1e18+1;
const ll P=100000;
il int read(){
    RG int data=0,w=1;RG char ch=getchar();
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(ch<='9'&&ch>='0')data=(data<<3)+(data<<1)+(ch^48),ch=getchar();
    return data*w;
}

il void file(){
    srand(time(NULL)+rand());
    freopen(FILE".in","r",stdin);
    freopen(FILE".out","w",stdout);
}

int n,rt,tot,cd,nxt[]={1,2,0};
struct node{int d[4];}p[N];
bool operator <(node x,node y){
    if(x.d[0]!=y.d[0])return x.d[0]<y.d[0];
    if(x.d[1]!=y.d[1])return x.d[1]<y.d[1];
    if(x.d[2]!=y.d[2])return x.d[2]<y.d[2];
    if(x.d[3]!=y.d[3])return x.d[3]<y.d[3];
    return 0;
}
struct kdnode{int d[3],l[3],r[3],s[2],v,mx;}t[N];
bool operator <(kdnode x,kdnode y){return x.d[cd]<y.d[cd];}

#define ls (t[i].s[0])
#define rs (t[i].s[1])
#define mid ((l+r)>>1)
int l[3],r[3],ans;
il void update(int i){
    for(RG int k=0;k<3;k++){
        t[i].l[k]=t[i].r[k]=t[i].d[k];
        if(ls){
            t[i].l[k]=min(t[i].l[k],t[ls].l[k]);
            t[i].r[k]=max(t[i].r[k],t[ls].r[k]);            
        }
        if(rs){
            t[i].l[k]=min(t[i].l[k],t[rs].l[k]);
            t[i].r[k]=max(t[i].r[k],t[rs].r[k]);            
        }
    }
    t[i].mx=max(t[i].v,max(t[ls].mx,t[rs].mx));
}

int rebuild(int l,int r,int D){
    if(l>r)return 0;
    cd=D;nth_element(t+l,t+mid,t+r+1);
    t[mid].s[0]=rebuild(l,mid-1,nxt[D]);
    t[mid].s[1]=rebuild(mid+1,r,nxt[D]);
    update(mid);return mid;
}

void insert(int &i,int D){
    if(!i){
        i=++tot;Cpy(t[i].d,r);Cpy(t[i].l,r);Cpy(t[i].r,r);
    }
    if(t[i].d[0]==r[0]&&t[i].d[1]==r[1]&&t[i].d[2]==r[2]){
        t[i].v=max(t[i].v,ans);t[i].mx=max(t[i].mx,t[i].v);return;
    }
    insert(t[i].s[r[D]>t[i].d[D]],nxt[D]);update(i);
}

#define ckm(dd) (l[dd]<=t[i].l[dd]&&t[i].r[dd]<=r[dd])
#define ckd(dd) (l[dd]<=t[i].d[dd]&&t[i].d[dd]<=r[dd])
#define cke(dd) (r[dd]<t[i].l[dd]||t[i].r[dd]<l[dd])
void query(int i,int D){
    if(!i||cke(2)||cke(1)||cke(0)||t[i].mx<=ans)return;
    if(ckm(2)&&ckm(1)&&ckm(0)){ans=max(ans,t[i].mx);return;}
    if(ckd(2)&&ckd(1)&&ckd(0))ans=max(ans,t[i].v);
    query(ls,nxt[D]);query(rs,nxt[D]);
}

int main()
{
    n=read();
    for(RG int i=1;i<=n;i++)for(RG int k=0;k<4;k++)p[i].d[k]=read();
    sort(p+1,p+n+1);
    for(RG int i=1;i<=n;i++){
        for(RG int k=0;k<3;k++)r[k]=p[i].d[k+1];
        ans=0;query(rt,0);ans++;
        insert(rt,0);
        if(i%10000==0)rt=rebuild(1,tot,0);
    }
    printf("%d\n",t[rt].mx);
    return 0;
}

转载于:https://www.cnblogs.com/cjfdf/p/9656384.html