BZOJ4825: [Hnoi2017]单旋（Splay）

题面

传送门

题解

调了好几个小时……指针太难写了……

因为只单旋最值，我们以单旋\(\min\)为例，那么\(\min\)是没有左子树的，而它旋到根之后，它的深度变为\(1\)，它的右子树里所有节点深度不变，其它所有节点都深度\(+1\)。那么这可以看做一个区间加和单点修改的事情，可以用\(Splay\)维护

然后就是插入节点，我们在\(Splay\)里找到它的前驱和后继，那么前驱的右儿子和后继的左儿子必定只有一个是空的，而且只有深度大的那个节点是空的，然后直接把节点插入就行了

还有一个问题就是我们该怎么找\(\min\)的右子树，我们可以对于每一个点维护它在题中所说\(Spaly\)中的深度，那么\(\min\)的右子树必定是从最左边开始的连续一段区间，且每一个节点深度都大于等于\(dep_{\min}\)，可以在\(Splay\)上二分找。所以还需要顺便维护一下区间最小深度

//minamoto
#include<bits/stdc++.h>
#define R register
#define inline __inline__ __attribute__((always_inline))
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
    R int res,f=1;R char ch;
    while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
    for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
    return res*f;
}
inline int getop(){R char ch;while((ch=getc())>'9'||ch<'0');return ch-'0';}
char sr[1<<21],z[20];int K=-1,Z=0;
inline void Ot(){fwrite(sr,1,K+1,stdout),K=-1;}
void print(R int x){
    if(K>1<<20)Ot();if(x<0)sr[++K]='-',x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++K]=z[Z],--Z);sr[++K]='\n';
}
const int N=1e5+5,inf=0x3f3f3f3f;
inline int min(R int x,R int y){return x<y?x:y;}
inline int max(R int x,R int y){return x>y?x:y;}
struct node;typedef node* ptr;
struct node{
    ptr lc,rc,fa;int sz,d,mn,t,s,c;
    inline node();
    inline void init(R int ss,R int dd){sz=1,s=ss,d=mn=dd,c=1;}
    inline void ppd(R int x){t+=x,d+=x,mn+=x;}
    inline void pd(){if(t)lc->ppd(t),rc->ppd(t),t=0;}
    inline ptr upd(){sz=lc->sz+rc->sz+c,mn=min(d,min(lc->mn,rc->mn));return this;}
    inline ptr son(R int x){return x<s?lc:rc;}
}e[N],*rt;int tot;
inline node::node(){lc=rc=fa=e;}
void rotate(ptr &rt,ptr p){
    ptr s=p->fa,t=s->fa;
    s!=rt?(t->lc==s?t->lc:t->rc)=p:rt=p;
    p->fa=t,s->fa=p;
    if(s->lc==p)s->lc=p->rc,p->rc->fa=s,p->rc=s->upd();
        else s->rc=p->lc,p->lc->fa=s,p->lc=s->upd();
}
void pd(ptr rt,ptr p){if(p!=rt)pd(rt,p->fa);p->pd();}
ptr splay(ptr &rt,ptr p){
    pd(rt,p);
    while(p!=rt){
        if(p->fa!=rt)rotate(rt,p->fa->lc==p^p->fa->fa->lc==p->fa?p:p->fa);
        rotate(rt,p);
    }
    return p->upd();
}
ptr push(int s,int d){
    ptr p=rt,las=e;
    while(p!=e)p->pd(),las=p,p=p->son(s);
    p=e+(++tot),p->init(s,d),p->fa=las;
    if(las!=e)(s<las->s?las->lc:las->rc)=p;
    return splay(rt,p);
}
ptr find(int s){
    ptr p=rt;
    while(p->son(s)!=e)p->pd(),p=p->son(s);
    return splay(rt,p);
}
ptr lst(int s){
    ptr p=find(s);
    if(p->s<s)return p;
    p->pd(),p=p->lc;while(p->rc!=e)p->pd(),p=p->rc;return p;
}
ptr nxt(int s){
    ptr p=find(s);
    if(p->s>s)return p;
    p->pd(),p=p->rc;while(p->lc!=e)p->pd(),p=p->lc;return p;
}
ptr Kth(int k){
    ptr p=rt;if(k>p->sz)return false;
    while(true){
        p->pd();
        if(p->lc->sz+1==k)return p;
        p=p->lc->sz>=k?p->lc:(k-=p->lc->sz+1,p->rc);
    }
}
ptr split(int l,int r){ptr s=Kth(l-1),t=Kth(r+1);splay(rt,s);return splay(s->rc,t)->lc;}
int getl(ptr p,int d){
    if(p==e)return 0;p->pd();
    int s;
    s=min(p->d,p->lc->mn)>=d?getl(p->rc,d)+p->lc->sz+1:getl(p->lc,d);
    return s;
}
int getr(ptr p,int d){
    if(p==e)return 0;p->pd();
    return min(p->d,p->rc->mn)>=d?getr(p->lc,d)+p->rc->sz+1:getr(p->rc,d);
}
int m,op,x,l,sz;ptr s,t,S,T;
int main(){
//  freopen("testdata.in","r",stdin);
    m=read(),rt=e,e->mn=e->d=inf,S=push(-inf,inf),T=push(inf,inf),sz=2;
    while(m--){
        op=getop();
        switch(op){
            case 1:{
                x=read(),++sz,s=lst(x),t=nxt(x);
                int D=max(s==S?0:s->d,t==T?0:t->d)+1;
                push(x,D),print(D);
                break;
            }
            case 2:{
                s=Kth(2),l=min(getl(rt,s->d),sz-1)-1;
                print(s->d),split(2,sz-1)->ppd(1);
                if(l>1)split(2,l+1)->ppd(-1);
                s=splay(rt,s),s->mn=s->d=1;
                break;
            }
            case 3:{
                s=Kth(sz-1),
                l=min(getr(rt,s->d),sz-1)-1;
                print(s->d),split(2,sz-1)->ppd(1);
                if(l>1)split(sz-l,sz-1)->ppd(-1);
                s=splay(rt,s),s->mn=s->d=1;
                break;
            }
            case 4:{
                s=Kth(2),l=min(getl(rt,s->d),sz-1)-1;
                print(s->d),split(2,sz-1)->ppd(1);
                if(l>1)split(2,l+1)->ppd(-1);
                s=splay(rt,s);
                s->lc->fa=e,s->rc->fa=s->lc,s->lc->rc=s->rc,rt=s->lc;
                rt->ppd(-1),--sz,rt->upd();
                break;
            }
            case 5:{
                s=Kth(sz-1),l=min(getr(rt,s->d),sz-1)-1;
                print(s->d),split(2,sz-1)->ppd(1);
                if(l>1)split(sz-l,sz-1)->ppd(-1);
                s=splay(rt,s);
                s->rc->fa=e,s->lc->fa=s->rc,s->rc->lc=s->lc,rt=s->rc;
                rt->ppd(-1),--sz,rt->upd();
                break;
            }
        }
    }
    return Ot(),0;
}

转载于:https://www.cnblogs.com/bztMinamoto/p/10735588.html