BZOJ 3506 机械排序臂 splay

    这是宝宝第二次敲splay了（如果LCT不算的话），还好吧！为什么自己写的TLE了呢？找个时间写下对拍吧！暂时还是要学习别人的代码的。一开始的思路是找到最小的点后，返回距离，然后删掉。在这里更深刻地体会到了自己的天真无邪，too young too simple ， （假设当前的根是我们要删除的）以为只要把左子节点的父亲修改为右子结点， 右子结点的左儿子设成左子结点，然后把右子结点更新为根就好。（有没有发现问题所在？，右子结点原来的左儿子怎么办？）那好问题出现了就要寻求解决的方法了，我是这样想的，在右子树里一直往左走，走到第一个没有左子结点的位置，返回其标号，然后我们再把它旋转到根，然后……（你懂的）。然而这个TLE了。 其实这还不是我最开始的我，最开始的我是这样的，弄了左右边界，然后发现边界都被我翻转了，（有没有感到我的愚蠢）。后来再次Orz了黄学长的代码，当然我自己是做了修改的。加油了。splay还是很重要的。（我不想做源程狗。）

#include<cstdio>
#include<iostream>
#include<algorithm>
#define lc c[k][0]
#define rc c[k][1]
#define rep(i,j,k) for(int  i = j;  i <= k; i++) 
#define maxn 100233
using namespace std;

struct node{
int v, pos;
bool operator < (const node&rhs) const{
    return pos < rhs.pos;
}
} nod[maxn];

bool cmp(node a,node b)
{
    return a.v < b.v || (a.v == b.v && a.pos < b.pos ); 
}

int read()
{
    int s = 0, t = 1; char c = getchar();
    while( !isdigit(c) ){
        if( c == '-' )  t= -1; c = getchar();
    }
    while( isdigit(c) ){
        s =s * 10 + c - '0'; c= getchar();
    }
    return s * t;
}

int s[maxn], pa[maxn], c[maxn][2], mn[maxn], v[maxn], root;
bool rev[maxn];

void pushdown(int k)
{
    if( rev[k] ){
        rev[k] ^= 1, rev[lc] ^= 1, rev[rc] ^= 1;
        swap(lc,rc);
    } 
}

void maintain(int k)
{
    mn[k] = k;
    if( v[mn[lc]] < v[mn[k]] ) mn[k] = mn[lc];
    if( v[mn[rc]] < v[mn[k]] ) mn[k] = mn[rc];
    s[k] = s[lc] + s[rc] + 1; 
}

void rorate(int k,int &root)
{
    int fa = pa[k], gfa = pa[fa];
    int l = c[fa][1] == k, r = l ^ 1;
    if( fa != root ){
        c[gfa][c[gfa][1] == fa] = k;
    } else root = k;
    pa[fa] = k, pa[c[k][r]] = fa, pa[k] = gfa;
    c[fa][l] = c[k][r], c[k][r] = fa; maintain(fa), maintain(k);
}

int q[maxn];
void splay(int k,int &root)
{
    int top = 0; q[++top] = k;
    for(int x = k; x != root; x = pa[x] ) q[++top] = pa[x];
    while( top ) pushdown(q[top--]);
    while( k != root ){
        int fa = pa[k], gfa = pa[fa];
        if( fa != root ) {
            if( c[gfa][1] == fa ^ c[fa][1] == k ) rorate(k,root);
            else rorate(fa,root);
        }
        rorate(k,root);
    }
}

void build(int l,int r,int pre)
{
    if( l > r ) return;
    int mid = (l+r)>>1;
    c[pre][pre < mid] = mid; pa[mid] = pre; v[mid] = nod[mid].v, mn[mid] = mid;
    c[mid][0] = c[mid][1] = 0;
    if( l < r ){
        build(l,mid-1,mid); build(mid+1,r,mid);
        maintain(mid);
    } else s[mid] = 1;
}

int rank(int x,int k)
{
    if( rev[k] ) pushdown(k);
    if( s[lc] >= x ) return rank(x,lc);
    else if( s[lc] + 1 < x ) return rank(x-s[lc]-1,rc);
    else return k;
}

int query_min(int x,int y)
{
    x = rank(x,root), y = rank(y+2,root);
    splay(x,root); splay(y,c[x][1]);
    return mn[c[y][0]]; 
}

void rever(int x,int y)
{
    x = rank(x,root), y = rank(y+2,root);
    splay(x,root); splay(y,c[x][1]);
    rev[c[y][0]] ^= 1;
}

void out(int k)
{
    if( lc ) out(lc);
    cout<<v[k]<<endl;
    if( rc ) out(rc);
}

int main()
{
    int n;
    while( scanf("%d", &n) && n ){
    v[0] = nod[1].v = nod[n+2].v = 0x7fffffff;
    rep(i,2,n+1) nod[i].v = read(), nod[i].pos = i;
    sort(nod+2,nod+2+n,cmp); root = (3+n)>>1;
    rep(i,2,n+1) nod[i].v = i-1;
    sort(nod+2,nod+2+n);  build(1,n+2,0);  //out(root);
    rep(i,1,n){
        int q = query_min(i,n); splay(q,root);
        printf("%d ", s[c[q][0]] );
        rever(i,s[c[q][0]]);
    }
            puts("");
    }
    return 0;
}