BJ模拟 传送门【KD树+拓扑排序】

题目描述：

解题思路：

首先想到的肯定是拓扑排序，但直接建边是$O(n^2)$的，考虑优化。
注意到这个图的性质：只能往右上方向走，所以我们可以从右上角到左下角倒着拓扑，用KD树去模拟拓扑排序的过程即可。

#include<bits/stdc++.h>
#define ll unsigned long long
using namespace std;
int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')c=getchar(),f=-1;
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}
const int N=200005;
struct node{int x,y,v,mnx,mny,mxx,mxy,mxv,lc,rc;}tr[N];
struct point{int x,y,a,b;}p[N];
inline bool cmpx(const point &a,const point &b){return a.x<b.x||a.x==b.x&&a.y<b.y;}
inline bool cmpy(const point &a,const point &b){return a.y<b.y||a.y==b.y&&a.x<b.x;}
int n,rt,ans,nowv;
node newnode(int x,int y)
{
    node t;
    t.x=t.mnx=t.mxx=x;
    t.y=t.mny=t.mxy=y;
    t.v=t.mxv=t.lc=t.rc=0;
    return t;
}
void update(int u)
{
    int lc=tr[u].lc,rc=tr[u].rc;
    tr[u].mnx=min(tr[u].x,min(tr[lc].mnx,tr[rc].mnx));
    tr[u].mxx=max(tr[u].x,max(tr[lc].mxx,tr[rc].mxx));
    tr[u].mny=min(tr[u].y,min(tr[lc].mny,tr[rc].mny));
    tr[u].mxy=max(tr[u].y,max(tr[lc].mxy,tr[rc].mxy));
    tr[u].mxv=max(tr[u].v,max(tr[lc].mxv,tr[rc].mxv));
}
void build(int &u,int l,int r,int d)
{
    int mid=l+r>>1;
    d?nth_element(p+l,p+mid,p+r+1,cmpy):nth_element(p+l,p+mid,p+r+1,cmpx);
    tr[u=mid]=newnode(p[mid].x,p[mid].y);
    if(l<mid)build(tr[u].lc,l,mid-1,d^1);
    if(r>mid)build(tr[u].rc,mid+1,r,d^1);
    update(u);
}
void query(int u,int mnx,int mny,int mxx,int mxy)
{
    if(tr[u].mxv<=nowv)return;
    if(tr[u].mxx<mnx||tr[u].mnx>mxx||tr[u].mxy<mny||tr[u].mny>mxy)return;
    if(tr[u].mnx>=mnx&&tr[u].mxx<=mxx&&tr[u].mny>=mny&&tr[u].mxy<=mxy){nowv=max(nowv,tr[u].mxv);return;}
    if(tr[u].x>=mnx&&tr[u].x<=mxx&&tr[u].y>=mny&&tr[u].y<=mxy)nowv=max(nowv,tr[u].v);
    query(tr[u].lc,mnx,mny,mxx,mxy),query(tr[u].rc,mnx,mny,mxx,mxy);
}
void modify(int u,int x,int y,int v,int d)
{
    if(tr[u].x==x&&tr[u].y==y){tr[u].v=max(tr[u].v,v);update(u);return;}
    int t;
    if(!d)t=(x<tr[u].x||x==tr[u].x&&y<tr[u].y);
    else t=(y<tr[u].y||y==tr[u].y&&x<tr[u].x);
    t?modify(tr[u].lc,x,y,v,d^1):modify(tr[u].rc,x,y,v,d^1);
    update(u);
}
int main()
{
    //freopen("lx.in","r",stdin);
    n=getint();
    for(int i=1;i<=n;i++)
        p[i].x=getint(),p[i].y=getint(),p[i].a=getint(),p[i].b=getint();
    build(rt,1,n,0);
    sort(p+1,p+n+1,cmpx);
    for(int i=n;i;i--)
    {
        nowv=0;
        query(rt,p[i].x,p[i].y,p[i].a,p[i].b);
        ans=max(ans,nowv+1);
        modify(rt,p[i].x,p[i].y,nowv+1,0);
    }
    cout<<ans<<'\n';
    return 0;
}