bzoj1552 [Cerc2007]robotic sort

http://www.elijahqi.win/2018/01/22/bzoj1552-cerc2007robotic-sort/
Description

Input
输入共两行,第一行为一个整数N,N表示物品的个数,1<=N<=100000。
第二行为N个用空格隔开的正整数,表示N个物品最初排列的编号。
Output
输出共一行,N个用空格隔开的正整数P1,P2,P3…Pn,Pi表示第i次操作前第i小的物品所在的位置。
注意:如果第i次操作前,第i小的物品己经在正确的位置Pi上,我们将区间[Pi,Pi]反转(单个物品)。
Sample Input
6
3 4 5 1 6 2
Sample Output
4 6 4 5 6 6
HINT

Source

HNOI2009集训Day6
思考良久确实没有什么思路 自己可能学了假的splay
题意:每次将编号第i小的位置和第i号位置之间整体交换一下 那么不妨splay建树的时候就把编号当下标来建树 然后splay排序的顺序就是按照他们的pos 那么每次我需要splay 当前这个序号到跟 然后和我的x位置交换一下


#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 110000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x;
}
struct node{
    int v,id,pos;
}a[N];
int q[N],n,size[N],fa[N],c[N][2],rev[N],root;
inline bool cmp1(const node &a,const node &b){
    return a.v==b.v?a.pos<b.pos:a.v<b.v;
}
inline bool cmp2(const node &a,const node &b){
    return a.pos<b.pos;
}
inline void update(int x){
    size[x]=size[c[x][0]]+size[c[x][1]]+1;
}
inline void pushdown(int x){
    if (!rev[x]) return;
    rev[c[x][0]]^=1;rev[c[x][1]]^=1;
    swap(c[x][0],c[x][1]);rev[x]=0;
}
inline void build(int f,int l,int r){
    if (l>r) return;
    if (l==r){
        size[a[l].id]=1;fa[a[l].id]=a[f].id;
        c[a[f].id][l>f]=a[l].id;return;
    }
    int mid=l+r>>1;build(mid,l,mid-1);build(mid,mid+1,r);
    fa[a[mid].id]=a[f].id;c[a[f].id][mid>f]=a[mid].id;
    update(a[mid].id);
}
inline void rotate(int x,int &tar){
    int y=fa[x],z=fa[y];
    int l=(c[y][1]==x),r=l^1;
    if (y==tar) tar=x;else c[z][c[z][1]==y]=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 &tar){
    int num=0;for (int i=x;i;i=fa[i]) q[++num]=i;
    for (int i=num;i;--i) pushdown(q[i]);
    while(x!=tar){
        int y=fa[x],z=fa[y];
        if (y!=tar){
            if((c[y][0]==x)^(c[z][0]==y)) rotate(x,tar);else rotate(y,tar);
        }rotate(x,tar);
    }
}
inline int find(int x,int sz){
    pushdown(x);
    int l=c[x][0],r=c[x][1];
    if (size[l]+1==sz) return x;
    if (sz<=size[l]) return find(l,sz);else return find(r,sz-size[l]-1);
}
inline void print(int x){
    pushdown(x);
    if (c[x][0]) print(c[x][0]);
    printf("%d %d %d %d\n",x,c[x][0],c[x][1],size[x]);
    if (c[x][1]) print(c[x][1]);
}
inline void solve(int x){
    //if (x!=n+1)print(root);
    splay(x,root);printf("%d",size[c[root][0]]);if (x!=n+1) printf(" ");
    int xx=find(root,x-1),yy=find(root,size[c[root][0]]+2);
    splay(yy,root);splay(xx,c[root][0]);rev[c[xx][1]]^=1;
}
int main(){
    freopen("bzoj1552.in","r",stdin);
        n=read();
        for (int i=1;i<=n;++i) a[i+1].v=read(),a[i+1].pos=i;
        sort(a+2,a+n+2,cmp1);
        for (int i=1;i<=n+2;++i) a[i].id=i;sort(a+2,a+n+2,cmp2);
        build(0,1,n+2);root=a[(3+n)>>1].id;
        for (int i=1;i<=n+2;++i) printf("%d ",a[i].id);
        for (int i=2;i<=n+1;++i) solve(i);
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值