bzoj1552 [Cerc2007]robotic sort(splay)

题目中的编号排完序其实就没嘛用了。。。主要就是一个区间翻转的操作。但是要找到数x的位置,把它转到根,就要提前把根到x的路径上的点都下放了,把x先做对,才能去转。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,a[N],fa[N],c[N][2],v[N],sz[N],rt,pos[N],sta[N];
bool rev[N];
struct Data{
    int val,id;
}aa[N];
inline bool cmp(Data a,Data b){return a.val==b.val?a.id<b.id:a.val<b.val;}
inline void update(int p){
    int l=c[p][0],r=c[p][1];
    sz[p]=sz[l]+sz[r]+1;
}
inline void build(int &p,int l,int r){
    int mid=l+r>>1;p=mid;v[p]=a[mid];
    if(l<mid) build(c[p][0],l,mid-1),fa[c[p][0]]=mid;
    if(r>mid) build(c[p][1],mid+1,r),fa[c[p][1]]=mid;
    update(p);
}
inline void pushdown(int p){
    if(!rev[p]) return;
    int l=c[p][0],r=c[p][1];
    if(l){
        rev[l]^=1;swap(c[l][0],c[l][1]);
    }if(r){
        rev[r]^=1;swap(c[r][0],c[r][1]);
    }rev[p]=0;
}
inline int find(int p,int x){
    pushdown(p);
    if(sz[c[p][0]]+1==x) return p;
    if(x<=sz[c[p][0]]) return find(c[p][0],x);
    return find(c[p][1],x-sz[c[p][0]]-1);
}
inline void rotate(int x,int &k){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(y==k) k=x;
    else c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x,int &k){
    int xx=x,top=0;
    while(xx) sta[++top]=xx,xx=fa[xx];
    while(top) pushdown(sta[top--]);//注意先把x做对。。。 
    while(x!=k){
        int y=fa[x],z=fa[y];
        if(y!=k){
            if(x==c[y][1]^y==c[z][1]) rotate(x,k);
            else rotate(y,k);
        }rotate(x,k);
    }
}
inline int split(int xx,int yy){
    int x=find(rt,xx),y=find(rt,yy);
    splay(x,rt);splay(y,c[x][1]);return c[y][0];
}
inline void rever(int l,int r){
    int x=split(l-1,r+1);
    rev[x]^=1;swap(c[x][0],c[x][1]);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i) aa[i].val=read(),aa[i].id=i;
    sort(aa+1,aa+n+1,cmp);a[1]=inf;a[n+2]=inf;
    for(int i=1;i<=n;++i) a[aa[i].id+1]=i,pos[i]=aa[i].id+1;
    build(rt,1,n+2);
    for(int i=1;i<=n;++i){
        splay(pos[i],rt);
        int ans=sz[c[rt][0]]+1;
        printf("%d",ans-1);if(i!=n) putchar(' ');
        rever(i+1,ans);
    }return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值