平衡树板子

Achen大佬说不要沉迷数据结构否则智商会降低的。

 

从省选考完后就开始学treap,首先是自己yy了一个打了两百多行,然后debug了2个月还是3个月记不清了。

最后弃疗,去找了网上别人的代码抄了一遍。

noip考完后补常规的一段时间,羡慕Achen能20分钟打出一个treap模板,于是自己也开始走上打板子的不归路。

到了后来可以10分钟左右打出一个结构体版的treap,看了Achen的数组版treap,觉得自己结构体版的太不优秀啦,于是就换成数组版的。

然后现在有几周没有碰过treap,感觉又不会打了。。。。不过还好,大概看一下又想起来了。

先放板子,还不会打题。

结构体版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<ctime>
using namespace std;
const int maxn=100000+10;
int n,tot=0,ans1,ans2,root;

int aa,ff;char cc;
int read() {
    aa=0;ff=1;cc=getchar();
    while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
    if(cc=='-') ff=-1,cc=getchar();
    while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
    return aa*ff;
}

struct Node{
    int num,rnd,son[2],sum,x;
}node[maxn];

void rotate(int &pos,int p) {
    int s=node[pos].son[p];
    node[s].sum=node[pos].sum;
    node[pos].son[p]=node[s].son[!p];
    node[s].son[!p]=pos;
    node[pos].sum=node[pos].x+node[node[pos].son[0]].sum+node[node[pos].son[1]].sum;
    pos=s;
}

void add(int &pos,int x) {
    if(!pos) {
        pos=++tot;node[pos].num=x;node[pos].rnd=rand();
        node[pos].sum=node[pos].x=1;
        return ;
    }
    node[pos].sum++;
    if(node[pos].num==x) {
        node[pos].x++; return;
    }
    int p=x>node[pos].num;
    add(node[pos].son[p],x);
    if(node[node[pos].son[p]].rnd<node[pos].rnd) rotate(pos,p);
}

void del(int &pos,int x) {
    if(!pos) return;
    if(node[pos].num==x) {
        if(node[pos].x>1) {
            node[pos].x--;node[pos].sum--;return;
        }
        if(node[pos].son[0]*node[pos].son[1]==0) {
            pos=node[pos].son[0]+node[pos].son[1]; return;
        }
        int p= node[node[pos].son[1]].rnd<node[node[pos].son[0]].rnd;
        rotate(pos,p);
        node[pos].sum--;
        del(node[pos].son[!p],x);
    }
    else {
        node[pos].sum--;
        if(node[pos].num>x) del(node[pos].son[0],x);
        else del(node[pos].son[1],x);
    }
}

int qrank(int pos,int x) {
    if(node[pos].num==x) return node[node[pos].son[0]].sum+1;
    if(node[pos].num>x) return qrank(node[pos].son[0],x);
    return node[node[pos].son[0]].sum+node[pos].x+qrank(node[pos].son[1],x);
}

int qnum(int pos,int x) {
    if(x>node[node[pos].son[0]].sum&&x<=node[node[pos].son[0]].sum+node[pos].x) return node[pos].num;
    if(x<=node[node[pos].son[0]].sum) return qnum(node[pos].son[0],x);
    return qnum(node[pos].son[1],x-node[pos].x-node[node[pos].son[0]].sum);
}

void q1(int pos,int x) {
    if(!pos) return;
    if(x>node[pos].num) {
        ans1=node[pos].num;
        q1(node[pos].son[1],x);
    }
    else q1(node[pos].son[0],x);
}

void q2(int pos,int x) {
    if(!pos) return;
    if(x<node[pos].num) {
        ans2=node[pos].num;
        q2(node[pos].son[0],x);
    }
    else q2(node[pos].son[1],x);
}

int main() {
    srand((unsigned)time(NULL));
    n=read();
    int opt,x,y;
    for(int i=1;i<=n;++i) {
        opt=read();x=read();
        if(opt==1) add(root,x);
        else if(opt==2) del(root,x);
        else if(opt==3) printf("%d\n",qrank(root,x));
        else if(opt==4) printf("%d\n",qnum(root,x));
        else if(opt==5) q1(root,x),printf("%d\n",ans1);
        else q2(root,x),printf("%d\n",ans2);
    }
    return 0;
}

数组版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
const int maxn=1e5+10;
int n,tot,root,ans;
int sum[maxn],num[maxn],rnd[maxn],fx[maxn],son[maxn][2];

int aa,ff;char cc;
int read() {
    aa=0;cc=getchar();ff=1;
    while(cc<'0'||cc>'9') {
        if(cc=='-') ff=-1;
        cc=getchar();
    }
    while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
    return aa*ff;
}

void rotate(int &pos,int p) {
    int s=son[pos][p];
    sum[s]=sum[pos];
    son[pos][p]=son[s][!p];
    son[s][!p]=pos;
    sum[pos]=sum[son[pos][0]]+sum[son[pos][1]]+fx[pos];
    pos=s;
}

void add(int &pos,int x) {
    if(!pos) {
        pos=++tot;
        num[pos]=x;
        rnd[pos]=rand();
    }
    sum[pos]++;
    if(num[pos]==x) {
        fx[pos]++;
        return;
    }
    int p=x>num[pos];
    add(son[pos][p],x);
    if(rnd[son[pos][p]]<rnd[pos]) rotate(pos,p);
}

void del(int &pos,int x) {
    if(!pos) return;
    if(num[pos]==x) {
        if(fx[pos]>1) {
            fx[pos]--;
            sum[pos]--;
            return;
        }
        if(son[pos][0]*son[pos][1]==0) {
            pos=son[pos][0]+son[pos][1];
            return;
        }
        int p=rnd[son[pos][1]]<rnd[son[pos][0]];
        rotate(pos,p);sum[pos]--;
        del(son[pos][!p],x);
    }
    else {
        sum[pos]--;
        del(son[pos][x>num[pos]],x);
    }
}

int qnum(int pos,int x) {
    if(x>sum[son[pos][0]]&&x<=sum[son[pos][0]]+fx[pos]) return num[pos];
    if(x<=sum[son[pos][0]]) return qnum(son[pos][0],x);
    return qnum(son[pos][1],x-sum[son[pos][0]]-fx[pos]);
}

int qrank(int pos,int x) {
    if(x==num[pos]) return sum[son[pos][0]]+1;
    if(x<num[pos]) return qrank(son[pos][0],x);
    return sum[son[pos][0]]+fx[pos]+qrank(son[pos][1],x);
}

void q1(int pos,int x) {
    if(!pos) return;
    if(num[pos]<x) ans=num[pos],q1(son[pos][1],x);
    else q1(son[pos][0],x);
}

void q2(int pos,int x) {
    if(!pos) return;
    if(num[pos]>x) ans=num[pos],q2(son[pos][0],x);
    else q2(son[pos][1],x);
}
int main() {
    srand(1031);
    n=read();
    int op,x;
    for(int i=1;i<=n;++i) {
        op=read();x=read();
        switch(op) {
            case 1:add(root,x);break;
            case 2:del(root,x);break;
            case 3:printf("%d\n",qrank(root,x));break;
            case 4:printf("%d\n",qnum(root,x));break;
            case 5:q1(root,x);printf("%d\n",ans);break;
            case 6:q2(root,x);printf("%d\n",ans);break;
        }
    }
    return 0;
}

 

而非旋treap支持split和merge很是有趣,这次直接照着Achen的板子学,感受到Achen大佬的数据结构多么优秀。

merge的时候不想用pair然后就乱搞了一搞,调了一阵子,最后强制把代码改成自己的码风。感觉不是很难打。

普通平衡树版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10,base=28;
int n,root,tot,ans;
int son[maxn][2],num[maxn],sum[maxn],rnd[maxn];

int aa,ff;char cc;
int read() {
	aa=0;ff=1;cc=getchar();
	while(cc!='-'&&(cc<'0'||cc>'9')) cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa*ff;
}

void ud(int pos){sum[pos]=sum[lc]+sum[rc]+1;}
ll pr(int x,int y) {return ((ll)x<<base)+(ll)y;}//第一个ll要套在里面 
int fi(ll x){return x>>base;}
int se(ll x){return x-((ll)fi(x)<<base);}

int merge(int x,int y) {
	if((ll)x*y==0) return x^y;
	if(rnd[x]<rnd[y]) return son[x][1]=merge(son[x][1],y),ud(x),x;
	else return son[y][0]=merge(x,son[y][0]),ud(y),y;
}

ll split(int pos,int x) {
	if(!pos) return 0;
	ll rs;
	if(sum[lc]>=x) {
		rs=split(lc,x);
		lc=se(rs);
		rs=pr(fi(rs),pos);
	}
	else {
		rs=split(rc,x-sum[lc]-1);
		rc=fi(rs);
		rs=pr(pos,se(rs));
	}
	return ud(pos),rs;
}

int qrank(int pos,int x) {
	if(!pos) return 1;
	if(num[pos]<x) return sum[lc]+1+qrank(rc,x);
	return qrank(lc,x);
}

int qnum(int pos,int x) {
	if(x==sum[lc]+1) return num[pos];
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-1);
}

void q1(int pos,int x) {
    if(!pos) return;
    if(num[pos]<x) ans=num[pos],q1(son[pos][1],x);
    else q1(son[pos][0],x);
}
 
void q2(int pos,int x) {
    if(!pos) return;
    if(num[pos]>x) ans=num[pos],q2(son[pos][0],x);
    else q2(son[pos][1],x);
}

void add(int &pos,int x) {
	int k=qrank(pos,x); ++tot;
	num[tot]=x;sum[tot]=1;rnd[tot]=rand();
	ll p=split(pos,k-1);
	pos=merge(fi(p),tot);
	pos=merge(pos,se(p));
}

void del(int &pos,int x) {
	int k=qrank(pos,x);
	ll p=split(pos,k-1);
	ll q=split(se(p),1);
	pos=merge(fi(p),se(q));
}

int main() {
	srand(1031);n=read();
	int op,x;
	while(n--) {
		op=read(); x=read();
		switch(op) {
			case 1:add(root,x);break;
			case 2:del(root,x);break;
			case 3:printf("%d\n",qrank(root,x));break;
			case 4:printf("%d\n",qnum(root,x));break;
			case 5:q1(root,x);printf("%d\n",ans);break;
            case 6:q2(root,x);printf("%d\n",ans);break;
		}
	}
	return 0;
}

文艺平衡树版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10,base=28;
int n,m,root,tot;
int son[maxn][2],num[maxn],sum[maxn],rnd[maxn],laz[maxn];

int aa;char cc;
int read() {
	aa=0;cc=getchar();
	while(cc<'0'||cc>'9') cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa;
}

void ud(int pos){sum[pos]=sum[lc]+sum[rc]+1;}
void pd(int pos){if(!laz[pos])return;swap(lc,rc);laz[lc]^=1;laz[rc]^=1;laz[pos]=0;}
ll pr(int x,int y) {return ((ll)x<<base)+(ll)y;}
int fi(ll x){return x>>base;}
int se(ll x){return x-((ll)fi(x)<<base);}

int bld(int l,int r) {
	if(l>r) return 0;
	int mid=(l+r)>>1,pos=++tot;
	num[pos]=mid;sum[pos]=1;rnd[pos]=rand();
	lc=bld(l,mid-1);rc=bld(mid+1,r);
	return ud(pos),pos;
}

int merge(int x,int y) {
	if((ll)x*y==0) return x^y; pd(x); pd(y);
	if(rnd[x]<rnd[y]) return son[x][1]=merge(son[x][1],y),ud(x),x;
	else return son[y][0]=merge(x,son[y][0]),ud(y),y;
}

ll split(int pos,int x) {
	if(!pos) return 0;
	ll rs; pd(pos);
	if(sum[lc]>=x) {
		rs=split(lc,x);
		lc=se(rs);
		rs=pr(fi(rs),pos);
	}
	else {
		rs=split(rc,x-sum[lc]-1);
		rc=fi(rs);
		rs=pr(pos,se(rs));
	}
	return ud(pos),rs;
}

void rever() {
	int l=read(),r=read();
	ll p=split(root,l-1);
	ll q=split(se(p),r-l+1);
	laz[fi(q)]^=1;
	root=merge(fi(p),fi(q));
	root=merge(root,se(q));
}

void print_ans(int pos) {
	if(!pos) return; pd(pos);
	print_ans(lc);
	printf("%d ",num[pos]);
	print_ans(rc);
}

int main() {
	srand(1031);
	n=read(); m=read();
	root=bld(1,n);
	while(m--) rever();
	print_ans(root);
	return 0;
}

 

 

Achen说我的Splay不是正常人的打法,我也觉得,但是主要原因是我打Splay的时候有打treap的感觉,然后函数都是递归的不优秀。

特别烦rotate函数感觉贼难受。还需要多多练习才能熟悉啊。

普通平衡树版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10;
int n,root,tot,ans;
int sum[maxn],fx[maxn],num[maxn],son[maxn][2],fa[maxn];

int aa,ff;char cc;
int read() {
	aa=0;ff=1;cc=getchar();
	while(cc<'0'||cc>'9') {
		if(cc=='-') ff=-1;
		cc=getchar();
	}
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa*ff;
}

void ud(int pos){sum[pos]=sum[lc]+sum[rc]+fx[pos];}
void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	fa[son[pos][!p]]=x;son[x][p]=son[pos][!p];
	fa[x]=pos;         son[pos][!p]=x;
	fa[pos]=y;         son[y][son[y][1]==x]=pos;
	ud(x);ud(pos);
	if(root==x) root=pos; fa[root]=fa[0]=0;
}

void splay(int pos,int &to) {
	int x,y;
	for(;pos!=to;rotate(pos)) {
		y=fa[x=fa[pos]];
		if(x!=to) (son[x][1]==pos ^ son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void add(int &pos,int x,int f) {
	if(!pos) {
		pos=++tot;
		fa[pos]=f;
		num[pos]=x;
	}
	++sum[pos];
	if(num[pos]==x) {
		++fx[pos];
		return;
	}
	add(son[pos][x>num[pos]],x,pos);
} 

void q1(int pos,int x) {
	if(!pos) return;
	if(num[pos]<x) ans=pos,q1(rc,x);
	else q1(lc,x);
}

void q2(int pos,int x) {
	if(!pos) return;
	if(num[pos]>x) ans=pos,q2(lc,x);
	else q2(rc,x);
}

int qnum(int pos,int x) {
	if(x>sum[lc]&&x<=sum[lc]+fx[pos]) return ans=pos,num[pos];
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-fx[pos]);
}

int qrank(int pos,int x) {
	if(x==num[pos]) return ans=pos,sum[lc]+1;
	if(x<num[pos]) return qrank(lc,x);
	return sum[lc]+fx[pos]+qrank(rc,x);
}

void del(int x) {
	int a,b; a=b=qrank(root,x);
	splay(ans,root); --sum[root];
	if((--fx[root])!=0) return;
	if(!son[root][0]) fa[root=son[root][1]]=0;
	else if(!son[root][1]) fa[root=son[root][0]]=0;
	else {
		qnum(root,a-1);a=ans;
		qnum(root,b);  b=ans;
		splay(b,root); splay(a,son[b][0]);
		fa[son[a][1]]=0;
		son[a][1]=0;
	}
}

int main() {
	 n=read(); int op,x;
    for(int i=1;i<=n;++i) {
        op=read(); x=read(); fa[ans=0]=sum[0]=fx[0]=0;
        switch(op) {
        	case 1:add(root,x,0);break;
            case 2:del(x);break;
            case 3:printf("%d\n",qrank(root,x));break;
            case 4:printf("%d\n",qnum(root,x)); break;
            case 5:q1(root,x);printf("%d\n",num[ans]);break;
            case 6:q2(root,x);printf("%d\n",num[ans]);break;
        }
        if(op!=2&&ans&&ans<=tot) splay(ans,root);
    }
	return 0;
}  

文艺平衡树版:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10;
int n,m,tot,root,ans;
int sum[maxn],num[maxn],son[maxn][2],fa[maxn],laz[maxn],zz[maxn];

int aa;char cc;
int read() {
	aa=0;cc=getchar();
	while(cc<'0'||cc>'9') cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa;
}

void ud(int pos) {sum[pos]=sum[lc]+sum[rc]+1;}
void pd(int pos) {if(!laz[pos])return;swap(lc,rc);laz[lc]^=1;laz[rc]^=1;laz[pos]=0;}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	fa[son[pos][!p]]=x;son[x][p]=son[pos][!p];
	fa[x]=pos;         son[pos][!p]=x;
	fa[pos]=y;         son[y][son[y][1]==x]=pos;
	ud(x); ud(pos);
	if(root==x) root=pos;fa[root]=fa[0]=0;
}

void splay(int pos,int &to) {
	int x,y,t=0;
	for(int p=pos;p!=to;p=fa[p]) zz[++t]=p;
	zz[++t]=to;
	while(t) pd(zz[t--]);
	for(;pos!=to;rotate(pos)) {
		y=fa[x=fa[pos]];
		if(x!=to) (son[x][1]==pos ^ son[y][1]==x)? rotate(pos):rotate(x);
	}
}

int bld(int l,int r,int f) {
	if(l>r) return 0;
	int pos=++tot,mid=(l+r)>>1;
	num[pos]=mid; sum[pos]=1; fa[pos]=f;
	lc=bld(l,mid-1,pos); rc=bld(mid+1,r,pos);
	return ud(pos),pos;
}

int qnum(int pos,int x) {
	pd(pos);if(x<1||x>sum[pos]) return ans=-1;
	if(x==sum[lc]+1) return num[ans=pos];
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-1);
}

void rever() {
	int l=read(),r=read(),a,b; ans=0;
	qnum(root,l-1); a=ans;
	qnum(root,r+1); b=ans;
	if(l==1&&r==n) laz[root]^=1;
	else if(l==1) splay(b,root),laz[son[b][0]]^=1;
	else if(r==n) splay(a,root),laz[son[a][1]]^=1;
	else {
		splay(b,root);
		splay(a,son[b][0]);
		laz[son[a][1]]^=1;
	}
}

void print_ans(int pos) {
    if(!pos) return; pd(pos);
    print_ans(lc);
    printf("%d ",num[pos]);
    print_ans(rc);
}

int main() {
	n=read(); m=read();
	root=bld(1,n,0);
	while(m--) rever();
	print_ans(root);
	return 0;
}

 

  

bzoj1503郁闷的出纳员

没有什么好说的,就是一个模板题,treap、非旋treap、splay都可以。

一个技巧:对于所有员工工资上涨下降可以直接改变工资的最低值,但是注意加入一个员工时对员工工资的处理,还有查询第k大输出-1的情况。

普通treap:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10;
int n,slry,now,tot,root,ans,delans;
char s[23];
int sum[maxn],num[maxn],fx[maxn],son[maxn][2],rnd[maxn];

int aa;char cc;
int read() {
	aa=0;cc=getchar();
	while(cc<'0'||cc>'9') cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa;
}

void rotate(int& pos,int p) {
	int s=son[pos][p];
	sum[s]=sum[pos];
	son[pos][p]=son[s][!p];
	son[s][!p]=pos;
	sum[pos]=sum[lc]+sum[rc]+fx[pos];
	pos=s;
}

void add(int& pos,int x) {
	if(!pos) {
		pos=++tot;
		num[pos]=x;
		rnd[pos]=rand();
	}
	++sum[pos];
	if(num[pos]==x) {
		++fx[pos];
		return;
	}
	int p=x>num[pos];
	add(son[pos][p],x);
	if(rnd[son[pos][p]]<rnd[pos]) rotate(pos,p); 
}

void del(int &pos,int x) {
	if(!pos) return;
	if(num[pos]==x) {
		if(fx[pos]>1) {
			--sum[pos];
			--fx[pos];
			return;
		}
		if((ll)son[pos][0]*son[pos][1]==0) {
			pos=son[pos][0]^son[pos][1];
			return;
		}
		int p=rnd[rc] < rnd[lc];
		rotate(pos,p); --sum[pos];
		del(son[pos][!p],x);
	}
	else {
		--sum[pos];
		del(son[pos][x>num[pos]],x);
	}
}

int qnum(int pos,int x) {
	if(x<=0) return -1+now-slry;
	if(x>sum[lc]&&x<=sum[lc]+fx[pos]) return num[pos];
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-fx[pos]);
}

int q(int pos) {
	if(lc) return q(lc);
	return num[pos];
}

void getadd(int x) {
	if(x<slry) return;
	add(root,x+now-slry);
}

void getdel(int x) {
	int num;
	while(root&&(num=q(root))<x) {
		++delans;
		del(root,num);
	}
}

int main() {
	srand(1031);
	n=read(); now=slry=read();
	for(int i=1;i<=n;++i) {
		scanf("%s",s+1);
		switch(s[1]) {
			case 'I':getadd(read());break;
			case 'A':now-=read(); break;
			case 'S':now+=read();getdel(now);break;
			case 'F':printf("%d\n",qnum(root,sum[root]-read()+1)-now+slry);break;
		}
	}
	printf("%d",delans);
	return 0;
}

非旋treap:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+10,base=25;
int n,slry,now,tot,root,ans,delans;
char s[23];
int sum[maxn],num[maxn],son[maxn][2],rnd[maxn];

int aa;char cc;
int read() {
	aa=0;cc=getchar();
	while(cc<'0'||cc>'9') cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	return aa;
}

void ud(int pos) {sum[pos]=sum[lc]+sum[rc]+1;}
ll pr(int x,int y) {return ((ll)x<<base)+(ll)y;}
int fi(ll x) {return x>>base;}
int se(ll x) {return x-((ll)fi(x)<<base);}

int merge(int x,int y) {
	if((ll)x*y==0) return x^y;
	if(rnd[x]<rnd[y]) return son[x][1]=merge(son[x][1],y),ud(x),x;
	return son[y][0]=merge(x,son[y][0]),ud(y),y;
}

ll split(int pos,int x) {
	if(!pos) return 0;
	ll rs;
	if(sum[lc]>=x) {
		rs=split(lc,x);
		lc=se(rs);
		rs=pr(fi(rs),pos);
	}
	else {
		rs=split(rc,x-sum[lc]-1);
		rc=fi(rs);
		rs=pr(pos,se(rs));
	}
	return ud(pos),rs;
}

int qrank(int pos,int x) {
	if(!pos) return 1;
	if(num[pos]<x) return sum[lc]+1+qrank(rc,x);
	return qrank(lc,x);
}

int qnum(int pos,int x) {
	if(x<=0) return -1+now-slry;
	if(x>sum[lc]&&x<=sum[lc]+1) return num[pos];
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-1);
}

void getadd(int &pos,int x) {
	if(x<slry) return;
	x+=now-slry;
    int k=qrank(pos,x); ++tot;
    num[tot]=x;sum[tot]=1;rnd[tot]=rand();
    ll p=split(pos,k-1);
    pos=merge(fi(p),tot);
    pos=merge(pos,se(p));
}

void getdel(int &pos,int x) {
	int k=qrank(pos,x);
	delans+=k-1;
	if(k==1) return;
	ll p=split(root,k-1);
	pos=se(p);
}

int main() {
	srand(1031);
	n=read(); now=slry=read();
	for(int i=1;i<=n;++i) {
		scanf("%s",s+1);
		switch(s[1]) {
			case 'I':getadd(root,read());break;
			case 'A':now-=read(); break;
			case 'S':now+=read();getdel(root,now);break;
			case 'F':printf("%d\n",qnum(root,sum[root]-read()+1)-now+slry);break;
		}
	}
	printf("%d",delans);
	return 0;
}

  

LCT

只要会Splay,LCT就很好打,第一次打是12月18日,因为直接是抄的板,所以1A,

昨天打SAM+LCT的时候发现有点快忘了LCT怎么打,就重新打了LCT,

然后link那里写成了newroot(y),fa[x]=y,所以死WA

然后用二分拼接代码法很快找到了bug

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=3e5+10;
int n,m,tot;
int fa[maxn],son[maxn][2],zz[maxn],sum[maxn],num[maxn],laz[maxn];

char cc;ll ff;
template<typename T>void read(T& aa) {
	aa=0;ff=1; cc=getchar();
	while(cc!='-'&&(cc<'0'||cc>'9')) cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

bool isroot(int pos) {return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}
void ud(int pos) {sum[pos]=sum[lc]^sum[rc]^num[pos];}
void pd(int pos) {
	if(!laz[pos]) return;
	laz[pos]=0; swap(lc,rc);
	if(lc) laz[lc]^=1;
	if(rc) laz[rc]^=1;
}

void rotate(int pos) {
	int x,y,p;y=fa[x=fa[pos]];
	p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos;fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

void splay(int pos) {
	int x,y,t=0;
	for(x=pos;!isroot(x);x=fa[x]) zz[++t]=x; zz[++t]=x;
	while(t) pd(zz[t--]);
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) 
			(son[x][1]==pos)^(son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos); rc=t; ud(pos);
	}
}

int findroot(int pos) {
	access(pos);
	splay(pos);
	while(lc) pos=lc;
	return pos;
}

void newroot(int pos) {
	access(pos);
	splay(pos);
	laz[pos]^=1;
}

void lk(int x,int y) {
	if(findroot(x)==findroot(y)) return ;
	newroot(x);
	fa[x]=y;
}

void del(int x,int y) {
	newroot(x);
	access(y);
	splay(y);
	if(son[y][0]==x) fa[x]=son[y][0]=0;
}

int Yth(int x,int y) {
	newroot(x);
	access(y);
	splay(y);
	return sum[y];
}

void mdf(int x,int y) {
	splay(x);
	num[x]=y;
	ud(x);
}

int main() {
	read(n); read(m);
	int op,x,y;
	for(int i=1;i<=n;++i) read(num[i]);
	for(int i=1;i<=m;++i) {
		read(op); read(x); read(y);
		switch(op) {
			case 0:printf("%d\n",Yth(x,y)); break;
			case 1:lk(x,y); break;
			case 2:del(x,y); break;
			case 3:mdf(x,y); break;
		}
	}
	return 0;
}

 

可持久化treap

就是在merge和split的时候新开一下节点就可以了

克服了心理障碍就没什么难打的

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define lc son[pos][0]
#define rc son[pos][1]
const ll maxn=1e6+7,base=25,INF=2147483647;
int n,tot,root[maxn];
int num[25*maxn],sum[25*maxn],rnd[25*maxn],son[25*maxn][2];

char cc;ll ff;
template<typename T>void read(T& aa) {
    aa=0;ff=1; cc=getchar();
    while(cc!='-'&&(cc<'0'||cc>'9')) cc=getchar();
    if(cc=='-') ff=-1,cc=getchar();
    while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
    aa*=ff;
}

ll pr(ll x,ll y) {return (x<<base)+y;}
ll fi(ll x) {return x>>base;}
ll se(ll x) {return x-(fi(x)<<base);}
void ud(int pos) {sum[pos]=sum[lc]+sum[rc]+1;}

void copy(int pos,int x) {
    num[pos]=num[x];
    sum[pos]=sum[x];
    rnd[pos]=rnd[x];
    lc=son[x][0]; rc=son[x][1];
}

bool exit(int pos,int x) {
    while(num[pos]!=x&&pos) {
        if(x<num[pos]) pos=lc;
        else pos=rc;
    }
    return pos;
}

int merge(int x,int y) {
    if((ll)x*y==0) return x^y;
    int pos=++tot;
    if(rnd[x]<rnd[y]) copy(pos,x),rc=merge(rc,y);
    else copy(pos,y),lc=merge(x,lc);
    return ud(pos),pos;
}

ll split(int p,int x) {
    if(!p) return 0;
    int pos=++tot; copy(pos,p);
    ll rs;
    if(sum[lc]>=x) {
        rs=split(lc,x);
        lc=se(rs);
        rs=pr(fi(rs),pos);
    }
    else {
        rs=split(rc,x-sum[lc]-1);
        rc=fi(rs);
        rs=pr(pos,se(rs));
    }
    return ud(pos),rs;
}

int qrank(int pos,int x) {
    if(!pos) return 1;
    if(num[pos]<x) return sum[lc]+1+qrank(rc,x);
    return qrank(lc,x);
}

int qnum(int pos,int x) {
    if(x==sum[lc]+1) return num[pos];
    if(x<=sum[lc]) return qnum(lc,x);
    return qnum(rc,x-sum[lc]-1);
}

void add(int &pos,int x) {
    int k=qrank(pos,x);
    ll p=split(pos,k-1);
    num[++tot]=x;sum[tot]=1;rnd[tot]=rand();
    pos=merge(fi(p),tot);
    pos=merge(pos,se(p));
}

void del(int &pos,int x) {
    if(!exit(pos,x)) return;
    int k=qrank(pos,x);
    ll p=split(pos,k-1);
    ll q=split(se(p),1);
    pos=merge(fi(p),se(q));
}

int q1(int pos,int x) {
    if(!pos) return -INF;
    if(num[pos]<x) return max(q1(rc,x),num[pos]);
    else return q1(lc,x);
}

int q2(int pos,int x) {
    if(!pos) return INF;
    if(num[pos]>x) return min(q2(lc,x),num[pos]);
    else return q2(rc,x);
}

int main() {
    srand(1031);
    read(n);int op,x,p;
    for(int i=1;i<=n;++i) {
        read(p); read(op); read(x);
        root[i]=root[p];
        switch(op) {
            case 1:add(root[i],x);break;
            case 2:del(root[i],x); break;
            case 3:printf("%d\n",qrank(root[i],x));break;
            case 4:printf("%d\n",qnum(root[i],x));break;
            case 5:printf("%d\n",q1(root[i],x));break;
            case 6:printf("%d\n",q2(root[i],x));break;
        }
        if(!root[i]) root[i]=merge(son[root[i]][0],son[root[i]][1]);
    }
    return 0;
}
/*
10
0 1 9
1 1 3
1 1 10
2 4 2
3 3 9
3 1 2
6 4 1
6 2 9
8 6 3
4 5 8
*/ 

 

树上莫队(uoj58糖果公园)

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
const int maxn=1e5+7,maxs=21;
int n,m,q,size;
ll nowans,v[maxn],w[maxn],c[maxn],ans[maxn];

char cc;ll ff;
template<typename T>void read(T& aa) {
	aa=0;ff=1; cc=getchar();
	while(cc!='-'&&(cc<'0'||cc>'9')) cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int dep[maxn],fa[maxn][maxs],dfn[maxn],dfsclock;
int bel[maxn],totb;
int zz[maxn],top;

struct Ask{
	int x,y,id,t;
	Ask(){}
	Ask(int x,int y,int id,int t):x(x),y(y),id(id),t(t){}
}ask[maxn];
int totask;
bool cmp(const Ask& a,const Ask& b) {
	return bel[a.x]==bel[b.x]? 
			(bel[a.y]==bel[b.y]? a.t<b.t:bel[a.y]<bel[b.y])
		:	bel[a.x]<bel[b.x];
}

struct Mdf{
	int pos,now,last;
	Mdf(){}
	Mdf(int pos,int now,int last):pos(pos),now(now),last(last){	}
}mdf[maxn];
int totmdf;

int fir[maxn],nxt[2*maxn],to[2*maxn],e=0;
void add(int x,int y) {
	to[++e]=y;nxt[e]=fir[x];fir[x]=e;
	to[++e]=x;nxt[e]=fir[y];fir[y]=e;
}

void together(int l,int &r) {
	++totb;  while(r>l) bel[zz[r--]]=totb;
}

void dfs(int pos,int d) {
	int bot=top,y,z;
	dfn[pos]=++dfsclock;
	dep[pos]=d;
	for(int i=1;i<20;++i) fa[pos][i]=fa[fa[pos][i-1]][i-1];
	for(y=fir[pos];y;y=nxt[y]) {
		if((z=to[y])==fa[pos][0]) continue;
		fa[z][0]=pos; dfs(z,d+1);
		if(top-bot>=size) together(bot,top);
	}
	zz[++top]=pos;
}

int get_lca(int x,int y) {
	if(dep[x]!=dep[y]) {
		if(dep[x]<dep[y]) swap(x,y);
		for(int i=19;(~i)&&dep[x]>dep[y];--i) 
			if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	}
	if(x==y) return x;
	for(int i=19;~i;--i) if(fa[x][i]!=fa[y][i]) {
		x=fa[x][i];y=fa[y][i];
	}
	return fa[x][0];
}

bool vis[maxn];int sum[maxn];
void chge(int pos,int p) {
	if(p) {
		sum[c[pos]]++;
		nowans+=(ll)w[sum[c[pos]]]*v[c[pos]];
	}
	else {
		nowans-=(ll)w[sum[c[pos]]]*v[c[pos]];
		sum[c[pos]]--;
	}
}

void fz(int x,int y) {
	int lca=get_lca(x,y);
	for(int i=x;i!=lca;i=fa[i][0]) chge(i,vis[i]^=1);
	for(int i=y;i!=lca;i=fa[i][0]) chge(i,vis[i]^=1);
}

int main() {
	read(n); read(m); read(q);
	size=0.5*pow(n,2.0/3);
	int op,x,y;
	for(int i=1;i<=m;++i) read(v[i]);
	for(int i=1;i<=n;++i) read(w[i]);
	for(int i=1;i<n;++i) {
		read(x); read(y);
		add(x,y);
	}
	dfs(1,1);
	while(top) bel[zz[top--]]=totb;
	for(int i=1;i<=n;++i) read(c[i]);
	for(int i=1;i<=q;++i) {
		read(op); read(x); read(y);
		if(!op) {mdf[++totmdf]=Mdf(x,y,c[x]);c[x]=y;}
		else {
			if(dfn[x]>dfn[y]) swap(x,y);
			++totask;
			ask[totask]=Ask(x,y,totask,totmdf);
		}
	}
	for(int i=totmdf;i;--i) c[mdf[i].pos]=mdf[i].last;
	
	sort(ask+1,ask+totask+1,cmp); 
	int nowt=0,pos,lastx=ask[1].x,lasty=ask[1].y;
	fz(lastx,lasty);
	for(int i=1;i<=totask;++i) {
		x=ask[i].x; y=ask[i].y;
		while(nowt>ask[i].t) {
			pos=mdf[nowt].pos;
			if(vis[pos]) chge(pos,0);
			c[pos]=mdf[nowt--].last;
			if(vis[pos]) chge(pos,1);
		}
		while(nowt<ask[i].t) {
			pos=mdf[++nowt].pos;
			if(vis[pos]) chge(pos,0);
			c[pos]=mdf[nowt].now;
			if(vis[pos]) chge(pos,1);
		}
		fz(lastx,x); fz(lasty,y);
		pos=get_lca(x,y);
		chge(pos,vis[pos]^=1);
		ans[ask[i].id]=nowans;
		chge(pos,vis[pos]^=1);
		lastx=x; lasty=y;
	}
	for(int i=1;i<=totask;++i) 	printf("%lld\n",ans[i]);
	return 0;
}
/*
4 3 5
1 9 2
7 6 5 1
2 3
3 1
3 4
1 2 3 2
1 1 2
1 4 2
0 2 1
1 1 2
1 4 2
*/

 

一些Achen几万年前就A过的模板题:

 

bzoj1208宠物收养所

treap/splay随便维护一下即可

注意bz上交不能有srand,为此debug2h

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+7;
const ll mod=1000000;
int n,root,now=-1,tot; ll ans,INF=1e14;
int son[maxn][2],rnd[maxn];ll num[maxn];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

void rotate(int &pos,int p) {
	int s=son[pos][p];
	son[pos][p]=son[s][!p];
	son[s][!p]=pos;
	pos=s;
}

void add(int &pos,ll x) {
	if(!pos) {
		pos=++tot;
		num[pos]=x;
		rnd[pos]=rand()*rand();
		return;
	}
	int p=x>num[pos];
	add(son[pos][p],x);
	if(rnd[son[pos][p]]<rnd[pos]) rotate(pos,p);
}

void del(int &pos,ll x) {
	if(!pos) return;
	if(num[pos]==x) {
		if((ll)lc*rc==0) {
			pos=lc^rc;
			return;
		}
		int p=rnd[rc]<rnd[lc];
		rotate(pos,p);
		del(son[pos][!p],x);
	}
	else del(son[pos][x>num[pos]],x);
}

ll q1(int pos,ll x) {
	if(!pos) return -INF;
	if(num[pos]<=x) return max(q1(rc,x),num[pos]);
	return q1(lc,x);
}

ll q2(int pos,ll x) {
	if(!pos) return INF;
	if(num[pos]>=x) return min(q2(lc,x),num[pos]);
	return q2(rc,x);
}

int main() {
	read(n); ll x,y,p1,p2;
	For(i,1,n) {
		read(x); read(y);
		if(x==now||(!root)) {
			add(root,y); now=x;
		}
		else {
			p1=q1(root,y); p2=q2(root,y);
			if(y-p1<=p2-y) del(root,p1),ans+=(y-p1);
			else del(root,p2),ans+=(p2-y);
			ans%=mod;
		}
	}
	printf("%lld",ans);
	return 0;
}

bzoj2733永无乡

splay+启发式合并

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+7;
int n,m;
char s[maxn];
int num[maxn],sum[maxn],fa[maxn],son[maxn][2];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int f[maxn];
int find(int x){return x==f[x]? x:f[x]=find(f[x]);}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p= son[x][1]==pos;
	fa[son[pos][!p]]=x; son[x][p]=son[pos][!p];
	fa[x]=pos; son[pos][!p]=x;
	fa[pos]=y; son[y][son[y][1]==x]=pos;
	sum[pos]=sum[x];
	sum[x]=sum[son[x][0]]+sum[son[x][1]]+1;
} 

void add(int &pos,int x,int f) {
	if(!pos) {
		pos=x;
		fa[pos]=f;
		return;
	}
	sum[pos]+=sum[x];
	add(son[pos][num[x]>num[pos]],x,pos);
}

void splay(int pos) {
	int x,y;
	for(;fa[pos];rotate(pos)) {
		y=fa[x=fa[pos]];
		if(y) (son[x][1]==pos ^ son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void insert(int root,int pos) {
	if(!pos) return;
	int lson=lc,rson=rc;
	fa[lc]=fa[rc]=0; lc=rc=0;
	sum[pos]=1; add(root,pos,0);
	insert(root,lson); insert(root,rson);
}

void lik(int x,int y) {
	if(find(x)==find(y)) return;
	f[find(x)]=find(y);
	splay(x);splay(y);
	if(sum[x]<sum[y]) swap(x,y);
	insert(x,y);
}

int qnum(int pos,int x) {
	if(x==sum[lc]+1) return pos;
	if(x<=sum[lc]) return qnum(lc,x);
	return qnum(rc,x-sum[lc]-1);
}

int q(int pos,int x) {
	splay(pos);
	if(x>sum[pos]) return -1;
	return qnum(pos,x);
}

int main() {
	read(n); read(m);
	int x,y;
	For(i,1,n) f[i]=i,read(num[i]=i),sum[i]=1;
	For(i,1,m) {
		read(x); read(y);
		lik(x,y);
	}
	read(m);
	For(i,1,m) {
		scanf("%s",s+1); read(x); read(y);
		if(s[1]=='B') lik(x,y);
		else printf("%d\n",q(x,y));
	}
	return 0;
}
/*
5 1           
4 3 2 5 1
1 2
7
Q 3 2
Q 2 1
B 2 3
B 1 5
Q 2 1
Q 2 4
Q 2 3
*/

bzoj2843极地旅行社

LCT

一开始1e5+7打+的时候按成了Backspace,结果打成1e7然后MLE了

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+7;
int n,m,f[maxn]; char s[maxn];
int fa[maxn],sum[maxn],son[maxn][2],num[maxn],laz[maxn];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int find(int x){return x==f[x]? x:f[x]=find(f[x]);}

bool isroot(int pos){return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}
void ud(int pos) {sum[pos]=sum[lc]+sum[rc]+num[pos];}
void pd(int pos) {
	if(!laz[pos]) return;
	swap(lc,rc); laz[pos]=0;
	if(lc) laz[lc]^=1;
	if(rc) laz[rc]^=1;
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[x][p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

int zz[maxn];
void splay(int pos) {
	int x,y,t=0;
	for(x=pos;!isroot(x);x=fa[x]) zz[++t]=x; zz[++t]=x;
	while(t) pd(zz[t--]);
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos ^ son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos); rc=t; ud(pos);
	}
}

void newroot(int pos) {
	access(pos);
	splay(pos);
	laz[pos]^=1;
}

void lik(int x,int y) {
	if(find(x)==find(y)) {
		printf("no\n"); return;
	}
	printf("yes\n");
	newroot(x); fa[x]=y; f[find(x)]=find(y);
}

void chge(int pos,int y) {
	splay(pos); 
	num[pos]=y; 
	ud(pos);
}

int q(int x,int y) {
	if(find(x)!=find(y)) return -1;
	newroot(x); access(y); splay(y);
	return sum[y];
}

int main() {
	read(n); int x,y,ans;
	For(i,1,n) read(num[i]),f[i]=i;
	read(m); 
	For(i,1,m) {
		scanf("%s",s+1); read(x); read(y);
		if(s[1]=='b') lik(x,y);
		else if(s[1]=='p') chge(x,y);
		else {
			ans=q(x,y);
			if(ans==-1) printf("impossible\n");
			else printf("%d\n",ans);
		}
	}
	return 0;
}
/*
5
4 2 4 5 6
10
excursion 1 1
excursion 1 2
bridge 1 2
excursion 1 2
bridge 3 4
bridge 3 5
excursion 4 5
bridge 1 3
excursion 2 4
excursion 2 5
*/

bzoj2002 弹飞绵羊

LCT:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=2e5+7;
int n,m,sum[maxn],fa[maxn],son[maxn][2];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

void ud(int pos) {sum[pos]=sum[lc]+sum[rc]+1;}
bool isroot(int pos){return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

void splay(int pos) {
	int x,y;
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos);
		rc=t;
		ud(pos);
	}
}

void cut(int x) {
	access(x); splay(x);
	fa[son[x][0]]=0; son[x][0]=0;
	ud(x);
}

int main() {
	read(n); int op,x,y;
	For(i,1,n+1) sum[i]=1;
	For(i,1,n) {
		read(x);
		fa[i]=min(i+x,n+1);
	}
	read(m);
	For(i,1,m) {
		read(op); read(x); x++;
		if(op==1) {
			access(x); 
			splay(x);
			printf("%d\n",sum[son[x][0]]);
		}
		else {
			read(y);
			y=min(y+x,n+1);
			cut(x); 
			fa[x]=y;
		}
	}
	return 0;
}

Splay维护括号序列:

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=4e5+7;
int n,m,root,num[maxn],sum[maxn],fa[maxn],son[maxn][2];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int fir[maxn],nxt[maxn],to[maxn],e=0;
void add(int x,int y) {
	to[++e]=y;nxt[e]=fir[x];fir[x]=e;
}

int dfn[maxn],end[maxn],dfn_clock;
void dfs(int pos) {
	num[dfn[pos]=++dfn_clock]=1;
	for(int y=fir[pos];y;y=nxt[y]) dfs(to[y]);
	num[end[pos]=++dfn_clock]=-1;
}

void ud(int pos) {sum[pos]=num[pos]+sum[lc]+sum[rc];}

void bld(int &pos,int l,int r,int f) {
	if(l>r) {pos=0;return;}
	pos=(l+r)>>1; fa[pos]=f;
	bld(lc,l,pos-1,pos);
	bld(rc,pos+1,r,pos);
	ud(pos);
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(y) son[y][son[y][1]==x]=pos; else root=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

void splay(int pos) {
	int x,y;
	for(;fa[pos];rotate(pos)) {
		y=fa[x=fa[pos]];
		if(y) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

int split(int pos) {
	splay(pos); int rs;
	fa[rs=rc]=0; rc=0;
	ud(pos); 
	return rs;
}

void merge(int x,int y) {
	splay(x); splay(y);
	while(son[x][1]) x=son[x][1];
	splay(x); fa[y]=x; son[x][1]=y;
	ud(x);
}

int q(int pos) {
	splay(pos);
	pos=lc;
	while(rc) pos=rc;
	return pos;
}

int main() {
	read(n); int op,x,y,l,r;
	For(i,1,n) {
		read(x);
		x=min(i+x,n+1);
		add(x,i);
	}
	dfs(n+1);
	bld(root,1,dfn_clock,0);
	
	son[2*n+3][0]=root; fa[root]=2*n+3;
	ud(2*n+3); root=2*n+3;
	son[2*n+4][1]=root; fa[root]=2*n+4;
	ud(2*n+4); root=2*n+4;
	
	read(m);
	For(i,1,m) {
		read(op); read(x); ++x;
		if(op==1) {
			splay(x=dfn[x]);
			printf("%d\n",sum[son[x][0]]);
		}
		else {
			read(y); y=min(y+x,n+1);
			l=q(dfn[x]);
			split(l);
			r=split(end[x]);
			merge(l,r);
			split(dfn[y]);
			merge(dfn[y],end[x]);
			merge(end[x],end[y]);
		}
	}
	return 0;
}

bzoj3786 星系探索

splay维护括号序列(下面是一份还T着的代码)

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=2e5+7;
int n,m,root;
char s[13];
ll w[maxn],d[maxn],tot[maxn],sum[maxn],num[maxn],laz[maxn];
int son[maxn][2],fa[maxn];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int fir[maxn],nxt[maxn],to[maxn],e=0;
void add(int x,int y) {
	to[++e]=y;nxt[e]=fir[x];fir[x]=e;
}

int dfn[maxn],end[maxn],dfn_clock;
void dfs(int pos) {
	num[dfn[pos]=++dfn_clock]=w[pos]; d[dfn_clock]=1;
	for(int y=fir[pos];y;y=nxt[y]) dfs(to[y]);
	num[end[pos]=++dfn_clock]=-w[pos]; d[dfn_clock]=-1;
}

void ud(int pos) {sum[pos]=num[pos]+sum[lc]+sum[rc];tot[pos]=d[pos]+tot[lc]+tot[rc];}
void add_tag(int pos,ll x) {
	sum[pos]+=tot[pos]*x;
	num[pos]+=d[pos]*x;//
	laz[pos]+=x;
}
void pd(int pos) {
	if(!laz[pos]) return;
	if(lc) add_tag(lc,laz[pos]);
	if(rc) add_tag(rc,laz[pos]);
	laz[pos]=0;
}

void bld(int &pos,int l,int r,int f) {
	if(l>r) {pos=0;return;}
	pos=(l+r)>>1; fa[pos]=f;
	bld(lc,l,pos-1,pos);
	bld(rc,pos+1,r,pos);
	ud(pos);
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(y) son[y][son[y][1]==x]=pos; else root=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

int zz[maxn];
void splay(int pos) {
	int t=0,x,y;
	for(int p=pos;p;p=fa[p]) zz[++t]=p;
	while(t) pd(zz[t--]);
	for(;fa[pos];rotate(pos)) {
		y=fa[x=fa[pos]];
		if(y) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

int split(int pos) {
	splay(pos); int rs;
	fa[rs=rc]=0; rc=0;
	ud(pos); return rs;
}

void merge(int x,int y) {
	splay(x); splay(y);//
	while(son[x][1]) pd(x),x=son[x][1];
	splay(x); fa[y]=x; son[x][1]=y;
	ud(x);
}

int q(int pos) {
	splay(pos); pos=lc;
	while(rc) pd(pos),pos=rc;
	return pos;
}

int main() {
	read(n); ll x,y,l,r;
	For(i,2,n) {
		read(x);
		add(x,i);
	}
	For(i,1,n) read(w[i]);
	dfs(1);
	bld(root,1,dfn_clock,0);
	
	son[2*n+3][0]=root; fa[root]=2*n+3;
	ud(2*n+3); root=2*n+3;
	son[2*n+4][1]=root; fa[root]=2*n+4;
	ud(2*n+4); root=2*n+4;
	
	read(m);
	For(i,1,m) {
		scanf("%s",s+1); read(x);
		if(s[1]=='Q') {
			splay(x=dfn[x]);
			printf("%lld\n",sum[son[x][0]]+num[x]);
		}
		else if(s[1]=='C') {
			read(y);
			l=q(dfn[x]);
			split(l);
			r=split(end[x]);
			split(dfn[y]);
			merge(dfn[y],dfn[x]);
			merge(dfn[x],end[y]);
		}
		else {
			read(y);
			l=q(dfn[x]);
			split(l);
			r=split(end[x]);
			splay(dfn[x]);
			add_tag(dfn[x],y);
			merge(l,dfn[x]);
			merge(dfn[x],r);
		}
	}
	return 0;
}
/*
3
1
1
4 5 7
5
Q 2
F 1 3
Q 2
C 2 3
Q 2
*/

bzoj4530 大融合

LCT维护子树信息,又把rotate的p和!p写错了,debug真烦

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=1e5+7;
int n,m,fa[maxn],sum[maxn],tot[maxn],son[maxn][2],laz[maxn];
char s[13];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

bool isroot(int pos) {return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}
void ud(int pos){sum[pos]=sum[lc]+sum[rc]+tot[pos]+1;}
void pd(int pos) {
	if(!laz[pos]) return;
	swap(lc,rc);
	laz[lc]^=1;
	laz[rc]^=1;
	laz[pos]=0;
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

int zz[maxn];
void splay(int pos) {
	int x,y,p,t=0;
	for(p=pos;!isroot(p);p=fa[p]) zz[++t]=p;
	zz[++t]=p;
	while(t) pd(zz[t--]);
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos);
		tot[pos]+=sum[rc]-sum[t];
		rc=t; ud(pos);
	}
}

void newroot(int pos) {
	access(pos);
	splay(pos);
	laz[pos]^=1;
}

void lk(int x,int y) {
	newroot(x); newroot(y);
	fa[y]=x; tot[x]+=sum[y];
	ud(x);
}

int main() {
	read(n); read(m); int x,y;
	For(i,1,n) sum[i]=1;
	For(i,1,m) {
		scanf("%s",s+1);
		read(x); read(y);
		if(s[1]=='A') lk(x,y);
		else {
			newroot(x); newroot(y);
			printf("%lld\n",(ll)(sum[y]-sum[x])*sum[x]);
		}
	}
	return 0;
}

UOJ207 共价大爷游长沙

LCT维护子树信息

num和sum没分清

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=3e5+7;
int n,m,son[maxn][2],fa[maxn],laz[maxn];
ll now,sum[maxn],num[maxn];

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

struct Node{
	int x,y;ll num;
	Node(){}
	Node(int x,int y,ll num):x(x),y(y),num(num){}
}node[maxn];
int tot;

bool isroot(int pos){return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}
void ud(int pos){sum[pos]=sum[lc]^sum[rc]^num[pos];}
void pd(int pos) {
	if(!laz[pos]) return;
	swap(lc,rc);
	if(lc) laz[lc]^=1;
	if(rc) laz[rc]^=1;
	laz[pos]=0;
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

int zz[maxn];
void splay(int pos) {
	int x,y,p,t=0;
	for(p=pos;!isroot(p);p=fa[p]) zz[++t]=p;
	zz[++t]=p;
	while(t) pd(zz[t--]);
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos);
		num[pos]^=(sum[rc]^sum[t]);//
		rc=t;
		ud(pos);
	}
}

void newroot(int pos) {
	access(pos);
	splay(pos);
	laz[pos]^=1;
}

void lk(int x,int y) {
	newroot(x); newroot(y);//
	fa[y]=x; num[x]^=sum[y];
	ud(x);
}

void del(int x,int y) {
	newroot(x); access(y); splay(y);
	son[y][0]=fa[x]=0;
	ud(y);
}

int main() {
	srand(1031);
	read(n); int op,x,y; ll r;
	read(n); read(m);
	For(i,1,n-1) {
		read(x); read(y);
		lk(x,y);
	}
	For(i,1,m) {
		read(op);
		if(op==1) {
			read(x); read(y);
			del(x,y);
			read(x); read(y);
			lk(x,y);
		}
		else if(op==2) {
			read(x); read(y);
			node[++tot]=Node(x,y,r=rand()*rand());
			newroot(x); num[x]^=r; ud(x);
			newroot(y); num[y]^=r; ud(y);
			now^=r;
		}
		else if(op==3) {
			read(r);
			x=node[r].x; y=node[r].y;
			r=node[r].num;
			newroot(x); num[x]^=r; ud(x);
			newroot(y); num[y]^=r; ud(y);
			now^=r;
		}
		else {
			read(x); read(y);
			newroot(x); newroot(y);
			printf(sum[x]==now?"YES\n":"NO\n");
		}
	}
	return 0;
}
/*
0
5 7
1 2
1 3
2 4
1 5
2 1 5
1 1 5 2 5
4 2 5
2 1 4
4 2 5
3 1
4 2 4
*/

bzoj4573大森林

感觉是一道非常有意思的题

当时看到这道题也觉得这种题肯定是把所有树的共同点一起维护了,然后对于不同的树的查询修改些地方

对于每次换生长节点就建一个虚点(权为0)

然后生长就在上一个建的虚点上长实点(权为1)

先把虚点串起来,考虑离线后,按照区间从左到右做

一个虚点会在一段时间有用,就把他连到对于实点的下面

否则他的父亲就是上一个虚点。

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=3e5+7;
int n,m,ans[maxn];
int son[maxn][2],fa[maxn],sum[maxn],num[maxn];
int L[maxn],R[maxn],id[maxn],tid,tot,last;

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

struct Ask{
	int pos,t,x,y;
	Ask(){}
	Ask(int pos,int t,int x,int y):pos(pos),t(t),x(x),y(y){}
	bool operator < (const Ask& b) const{return pos==b.pos? t<b.t:pos<b.pos;}
}ask[2*maxn];
int totask,totq;

void ud(int pos){sum[pos]=sum[lc]+sum[rc]+num[pos];}
int newnode(int x) {num[++tot]=x; ud(tot); return tot;}
bool isroot(int pos) {return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

void splay(int pos) {
	for(int x,y;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

int access(int pos) {
	int t=0;
	for(;pos;pos=fa[t=pos]) {
		splay(pos);
		rc=t;
		ud(pos);
	}
	return t;
}

void lk(int x,int y) {
	splay(x);
	fa[x]=y;
}

void cut(int pos) {
	access(pos);
	splay(pos);
	fa[lc]=0; lc=0;
	ud(pos);
}

int Yth(int x,int y) {
	int rs=0,lca;
	access(x); splay(x); rs+=sum[x];
	lca=access(y); splay(y); rs+=sum[y];
	access(lca); splay(lca); rs-=2*sum[lca];
	return rs;
}

int main() {
	read(n); read(m);
	id[++tid]=newnode(1);L[1]=1;R[1]=n;
	newnode(0);lk(2,1);last=2;
	int op,l,r,x;
	For(i,1,m) {
		read(op);
		if(op==0) {
			read(l); read(r);
			id[++tid]=newnode(1);
			L[tot]=l; R[tot]=r;
			lk(tot,last);
		}
		else if(op==1) {
			read(l); read(r); read(x);
			l=max(l,L[id[x]]); r=min(r,R[id[x]]);
			if(l>r) continue;
			lk(newnode(0),last);
			ask[++totask]=Ask(l,i-m,tot,id[x]);
			ask[++totask]=Ask(r+1,i-m,tot,last);
			last=tot;
		}
		else {
			read(x); read(l); read(r);
			ask[++totask]=Ask(x,++totq,id[l],id[r]);
		}
	}
	sort(ask+1,ask+totask+1);
	For(i,1,totask) {
		if(ask[i].t<=0) {
			cut(ask[i].x);
			lk(ask[i].x,ask[i].y);
		}
		else ans[ask[i].t]=Yth(ask[i].x,ask[i].y);
	}
	For(i,1,totq) printf("%d\n",ans[i]);
	return 0;
}
/*
5 5
0 1 5
1 2 4 2
0 1 4
2 1 1 3
2 2 1 3
*/

bzoj3669魔法森林

没有删除操作的动态MST问题

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=2e5+7,INF=0x3f3f3f3f;
int n,m,f[maxn],ans=INF;
int num[maxn],maxnum[maxn],son[maxn][2],fa[maxn],laz[maxn],tot;

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

int find(int x) {return x==f[x]? x:f[x]=find(f[x]);}

struct Node{
	int x,y,a,b;
	bool operator < (const Node& r) const{return a==r.a? b<r.b:a<r.a;}
}node[maxn];

void ud(int pos){
	maxnum[pos]=num[maxnum[lc]]>num[maxnum[rc]]? maxnum[lc]:maxnum[rc];
	if(num[pos]>num[maxnum[pos]]) maxnum[pos]=pos;
}
bool isroot(int pos) {return son[fa[pos]][0]!=pos&&son[fa[pos]][1]!=pos;}
int newnode(int r) {num[++tot]=r; ud(tot); return tot;}
void pd(int pos) {
	if(!laz[pos]) return ;
	swap(lc,rc);
	if(lc) laz[lc]^=1;
	if(rc) laz[rc]^=1;
	laz[pos]=0;
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(!isroot(x)) son[y][son[y][1]==x]=pos; fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

int zz[maxn];
void splay(int pos) {
	int x,y,t=0,p;
	for(p=pos;!isroot(p);p=fa[p]) zz[++t]=p; zz[++t]=p;
	while(t) pd(zz[t--]);
	for(;!isroot(pos);rotate(pos)) {
		y=fa[x=fa[pos]];
		if(!isroot(x)) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

void access(int pos) {
	for(int t=0;pos;pos=fa[t=pos]) {
		splay(pos);
		rc=t;
		ud(pos);
	}
}

void newroot(int pos) {
	access(pos);
	splay(pos);
	laz[pos]^=1;
}

void lk(int x,int y,int r) {
	int p=newnode(r);
	newroot(x); newroot(y);
	fa[x]=fa[y]=p;
}

void cut(int pos) {
	newroot(pos);
	splay(pos);
	fa[lc]=fa[rc]=0;
}

int Yth(int x,int y) {
	if(find(x)!=find(y)) return -1;
	newroot(x); 
	access(y); splay(y);
	return maxnum[y];
}

int main() {
	read(n); read(m);
	For(i,1,n) newnode(0),f[i]=i;
	For(i,1,m) {
		read(node[i].x);
		read(node[i].y);
		read(node[i].a);
		read(node[i].b);
	}
	sort(node+1,node+m+1);
	int x,y,pos;
	For(i,1,m) {
		x=node[i].x;
		y=node[i].y;
		if(find(x)!=find(y)) {
			f[find(x)]=find(y);
			lk(x,y,node[i].b);
		}
		else {
			pos=Yth(x,y);
			if(num[pos]>node[i].b) {
				cut(pos);
				lk(x,y,node[i].b);
			}
		}
		pos=Yth(1,n);
		if(~pos) ans=min(ans,num[pos]+node[i].a);
	}
	if(ans!=INF) printf("%d",ans);
	else printf("-1");
	return 0;
}
/*
4 5
1 2 19 1
2 3 8 12
2 4 12 15
1 3 17 8
3 4 1 17
*/

bzoj3595方伯伯的OJ

只要不卡常,什么都好,map+set+splay

//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<map>
#include<set>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define lc son[pos][0]
#define rc son[pos][1]
const int maxn=6e5+7;
const ll Bs=1e7;
int n,m,tot,root;
int L[maxn],R[maxn],sum[maxn],num[maxn],fa[maxn],son[maxn][2];
map<int,int> G,H;

struct Node{
	int l,r,id;
	Node(int l,int r,int id):l(l),r(r),id(id){}
	bool operator < (const Node& b) const{return l==b.l? r<b.r : l<b.l;}
};

set<Node> P;
set<Node>::iterator it;

char cc; ll ff;
template<typename T>void read(T& aa) {
	aa=0;cc=getchar();ff=1;
	while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
	if(cc=='-') ff=-1,cc=getchar();
	while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
	aa*=ff;
}

ll pr(ll x,ll y) {return x*Bs+y;}
ll fi(ll x){return x/Bs;}
ll se(ll x){return x%Bs;}

void ud(int pos) {num[pos]=R[pos]-L[pos]+1; sum[pos]=num[pos]+sum[lc]+sum[rc];}
int newnode(int l,int r,int f) {
	fa[++tot]=f;
	L[tot]=l; R[tot]=r;
	P.insert(Node(l,r,tot));
	ud(tot); 
	if(f) ud(f);
	return tot;
}
void del(int pos) {
	Node o=Node(L[pos],R[pos],pos);
	it=P.lower_bound(o);
	P.erase(it);
	lc=rc=fa[pos]=0;
	num[pos]=sum[pos]=0;
}
int get_id(int x) {
	Node o=Node(x,x+1,0);
	it=P.lower_bound(o);
	if(it->l>x) --it;
	return it->id;
}

void rotate(int pos) {
	int x,y,p; y=fa[x=fa[pos]]; p=son[x][1]==pos;
	if(y) son[y][son[y][1]==x]=pos; 
	if(x==root) root=pos;
	fa[pos]=y;
	son[x][p]=son[pos][!p]; fa[son[pos][!p]]=x;
	son[pos][!p]=x; fa[x]=pos;
	ud(x); ud(pos);
}

void splay(int pos) {
	for(int x,y;fa[pos];rotate(pos)) {
		y=fa[x=fa[pos]];
		if(y) (son[x][1]==pos^son[y][1]==x)? rotate(pos):rotate(x);
	}
}

int merge(int x,int y) {
	if((ll)x*y==0) return x+y;
	splay(x); splay(y);
	while(son[x][1]) x=son[x][1];
	splay(x);
	son[x][1]=y; fa[y]=x;
	ud(y); ud(x);
	return x;
}

ll split(int pos,int x) {
	splay(pos);
	int ld=lc,rd=rc,l,r;
	del(pos);
	l=newnode(L[pos],x,0);
	r=newnode(x+1,R[pos],0);
	son[l][0]=ld; fa[ld]=l;
	son[r][1]=rd; fa[rd]=r;
	ud(l); ud(r);
	root=merge(l,r);
	return pr(l,r);
}

void cut(int pos) {
	splay(pos);
	fa[lc]=fa[rc]=0;
	root=merge(lc,rc);
	lc=rc=0;
}

int sep(int pos,int x) {
	if(L[pos]==x&&R[pos]==x) return pos;
	if(L[pos]!=x) pos=se(split(pos,x-1));
	if(R[pos]!=x) pos=fi(split(pos,x));
	return pos;
}

int qrank(int x) {
	int pos=get_id(x);
	splay(pos);
	return sum[lc]+(x-L[pos]+1);
}

int qnum(int pos,int &x) {
	if(x>sum[lc]&&x<=sum[lc]+num[pos]) return x-=sum[lc],pos;
	if(x<=sum[lc]) return qnum(lc,x);
	x-=sum[lc]+num[pos]; return qnum(rc,x);
}

int main() {
	read(n); read(m);
	newnode(1,n,0); root=1;
	son[1][0]=newnode(1,0,1);
	son[1][1]=newnode(n+1,n,1);
	int op,x,y,pos,p,last=0;
	For(i,1,m) {
		read(op); read(x); x-=last;
		if(op==1) {
			read(y); y-=last;
			G[y]=G[x]? G[x]:x; G[x]=0;
			H[G[y]]=y;
			last=qrank(G[y]);
		}
		else {
			if(op!=4&&G[x]) x=G[x];
			if(op==2) {
				pos=get_id(x);
				last=qrank(x);
				pos=sep(pos,x);
				cut(pos); cut(2);
				root=merge(pos,root);
				root=merge(2,pos);
			}
			else if(op==3) {
				pos=get_id(x);
				last=qrank(x);
				pos=sep(pos,x);
				cut(pos); cut(3);
				root=merge(root,pos);
				root=merge(root,3);
			}
			else {
				pos=qnum(root,x);
				x=x+L[pos]-1;
				last=H[x]? H[x]:x;
			}
		}
		printf("%d\n",last);
	}
	return 0;
}
/*
10 10
1 2 11
3 13
2 5
3 7
2 8
2 10
2 11
3 14
2 18
4 9
*/

  

posted @ 2017-12-25 19:46  shixinyi  阅读(314)  评论(0编辑  收藏  举报