Bzoj1552 [cerc2007]robotic sort

阅读全文大概需要 1分钟
本文总阅读量
Posted by yjjr's blog on December 19, 2017

标签:splay

题目

本题同BZOJ3506

题目传送门

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模板(区间翻转)+维护子树内最小值

找到最小值把它翻转到根root的位置然后它左子树节点个数就是当前所处位置

code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
#define mem(x,num) memset(x,num,sizeof x)
#define inf 0x7fffffff
using namespace std;
inline ll read()
{
    ll f=1,x=0;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;
}
const int maxn=1e5+6;
struct node{int val,pos;}a[maxn];
inline bool operator < (node a,node b){return a.pos<b.pos;}
inline bool cmp (node a,node b){return (a.val<b.val)||(a.val==b.val&&a.pos<b.pos);}
 
int n,root,top,ans[maxn],s[maxn],fa[maxn],c[maxn][2],v[maxn],mn[maxn],size[maxn],pos[maxn];
bool rev[maxn];
 
 
void update(int x)
{
    int l=c[x][0],r=c[x][1];
    size[x]=size[l]+size[r]+1;
    mn[x]=v[x];pos[x]=x;
    if(mn[l]<mn[x]){mn[x]=mn[l];pos[x]=pos[l];}
    if(mn[r]<mn[x]){mn[x]=mn[r];pos[x]=pos[r];}
}
 
void pushdown(int x)
{
    int l=c[x][0],r=c[x][1];
    rev[x]^=1;rev[l]^=1;rev[r]^=1;
    swap(c[x][0],c[x][1]);
}
 
void rotate(int x,int &k)
{
    int y=fa[x],z=fa[y],l,r;
    if(c[y][0]==x)l=0;else l=1;r=l^1;
    if(y==k)k=x;
    else{if(c[z][0]==y)c[z][0]=x;else c[z][1]=x;}
    fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
    c[y][l]=c[x][r];c[x][r]=y;
    update(y);update(x);
}
 
void splay(int x,int &k)
{
    top=0;s[++top]=x;
    for(int i=x;fa[i];i=fa[i])s[++top]=fa[i];
    dep(i,top,1)
        if(rev[s[i]])pushdown(s[i]);
    while(x!=k){
        int y=fa[x],z=fa[y];
        if(y!=k){
            if(c[y][0]==x^c[z][0]==y)rotate(x,k);
            else rotate(y,k);
        }
        rotate(x,k);
    }
}
 
void build(int l,int r,int f)
{
    if(l>r)return;
    if(l==r){
        fa[l]=f;size[l]=1;
        if(l<f)c[f][0]=l;else c[f][1]=l;
        mn[l]=v[l]=a[l].val;pos[l]=l;
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid-1,mid);build(mid+1,r,mid);
    fa[mid]=f;v[mid]=a[mid].val;update(mid); 
    if(mid<f)c[f][0]=mid;else c[f][1]=mid;
}
 
int find(int x,int rank)
{
    if(rev[x])pushdown(x);
    int l=c[x][0],r=c[x][1];
    if(size[l]+1==rank)return x;
    else if(size[l]>=rank)return find(l,rank);
    else return find(r,rank-size[l]-1);
}
 
void rever(int l,int r)
{
    int x=find(root,l),y=find(root,r+2);
    splay(x,root);splay(y,c[x][1]);
    int z=c[y][0];
    rev[z]^=1;
}
 
int querymn(int l,int r)
{
    int x=find(root,l),y=find(root,r+2);
    splay(x,root);splay(y,c[x][1]);
    int z=c[y][0];
    return pos[z];
}
 
int main()
{
    n=read();root=(n+3)>>1;
    a[1].val=a[n+2].val=inf;mn[0]=inf; 
    rep(i,2,n+1)a[i].val=read(),a[i].pos=i;
    sort(a+2,a+2+n,cmp);
    rep(i,2,n+1)a[i].val=i;
    sort(a+2,a+2+n);
    build(1,n+2,0);
    rep(i,1,n){
        int x=querymn(i,n);
        splay(x,root);
        ans[i]=size[c[x][0]];
        rever(i,ans[i]);
    }
    rep(i,1,n-1)printf("%d ",ans[i]);printf("%d",ans[n]);
    return 0;
}
本文可以转载,但必须附上原文链接,否则你会终生找不到妹子!!!欢迎关注我的CSDN: ahyjjr