BZOJ 3223 文艺平衡树 （splay）

题目链接：BZOJ 3223

因为这几天被数据结构虐的各种不忍直视，所以我计划在这一段时间专攻数据结构。但是今天做这道很水的splay还是debug了很久QAQ。感觉这段时间状态好差= =。

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define maxn (100000+10)
int N,M,root,siz=0;
int c[maxn][2],tag[maxn],size[maxn],fa[maxn],id[maxn];

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;
}

void pushup(int x){
	int l=c[x][0],r=c[x][1];
	size[x]=size[l]+size[r]+1;
}

void pushdown(int x){
	int l=c[x][0],r=c[x][1];
	if(tag[x]){
		swap(c[x][0],c[x][1]);
		tag[x]=0; tag[r]^=1; tag[l]^=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[c[x][r]]=y; fa[x]=z; fa[y]=x;
	c[y][l]=c[x][r]; c[x][r]=y;
	pushup(y); pushup(x);
}

void splay(int x,int &k){
	while(x!=k){
		int y=fa[x],z=fa[y];
		if(y!=k){
			if((c[y][0]==x)==(c[z][0]==y))rotate(y,k);
			else rotate(x,k);
		}
		rotate(x,k);
	}
}

void built(int l,int r,int f){
	if(l>r)return ;
	int now=id[l],last=id[f];
	if(l==r){
		fa[now]=last; size[now]=1;
		if(l<f)c[last][0]=now; else c[last][1]=now;
		return ;
	}
	int mid=(l+r)>>1; now=id[mid];
	built(l,mid-1,mid); built(mid+1,r,mid);
	fa[now]=last; pushup(mid);
	if(mid<f)c[last][0]=now; else c[last][1]=now;
}

int find(int x,int rk){
	pushdown(x);
	int l=c[x][0],r=c[x][1];
	if(size[l]+1==rk)return x;
	else if(size[l]>=rk)return find(l,rk);
	else return find(r,rk-size[l]-1);
}

void rever(int l,int r){
	int x,y,z;
	x=find(root,l); y=find(root,r+2);
	splay(x,root); splay(y,c[x][1]);
	z=c[y][0]; tag[z]^=1;
}

int main(){	
	N=read(); M=read();
	for(int i=1;i<=N+2;i++)id[i]=++siz;
	built(1,N+2,0); root=(N+3)>>1;
	for(int i=1;i<=M;i++){
		int l=read(),r=read();
		rever(l,r);
	}
	for(int i=2;i<=N+1;i++){
		printf("%d ",find(root,i)-1);
	}
	return 0;
}