[BZOJ3600]没有人的算术

咕了这么久终于把这题补回来了...

关键在于怎么维护“数”的大小，直接用平衡树即可，每个平衡树节点维护的是一个实数区间$[l,r]$，那么$mid$就是这个节点所代表的“数”的大小，也即权值

如果我们使用的平衡树在插入后需要调整树的形态，那么被调整的部分的权值都需要被重新计算，如果平衡树深度太大，那么精度会爆，所以要用重量平衡树（修改影响到的子树大小和平衡树深度都不会太大）

这里用替罪羊树，原生支持重新计算权值，过几天补个旋转treap咕咕咕

btw...替罪羊树是不用记父亲的...是我智障了

#include<stdio.h>
typedef double du;
const du al=.75;
int max(int a,int b){return a>b?a:b;}
du rk[500010];
struct num{
	int l,r;
	num(int a=0,int b=0){l=a;r=b;}
}v[500010];
bool operator>(num a,num b){return rk[a.l]==rk[b.l]?rk[a.r]>rk[b.r]:rk[a.l]>rk[b.l];}
bool operator==(num a,num b){return rk[a.l]==rk[b.l]&&rk[a.r]==rk[b.r];}
int l[500010],r[500010],siz[500010],M,sc,rt,fs;
du sl,sr;
int insert(int&x,du lv,du rv,num d){
	du mv=(lv+rv)*.5;
	int t;
	if(x==0){
		x=++M;
		siz[x]=1;
		v[x]=d;
		rk[x]=mv;
		return x;
	}
	if(v[x]==d)return x;
	t=v[x]>d?insert(l[x],lv,mv,d):insert(r[x],mv,rv,d);
	siz[x]=siz[l[x]]+siz[r[x]]+1;
	if(sc&&(l[x]==sc||r[x]==sc))fs=x;
	if(siz[x]*al<max(siz[l[x]],siz[r[x]])){
		sc=x;
		sl=lv;
		sr=rv;
	}
	return t;
}
int tmp[500010],C;
void dfs(int x){
	if(l[x])dfs(l[x]);
	tmp[++C]=x;
	if(r[x])dfs(r[x]);
}
int build(int L,int R,du lv,du rv){
	int mid=(L+R)>>1,x=tmp[mid];
	du mv=(lv+rv)*.5;
	rk[x]=mv;
	l[x]=L<mid?build(L,mid-1,lv,mv):0;
	r[x]=mid<R?build(mid+1,R,mv,rv):0;
	siz[x]=siz[l[x]]+siz[r[x]]+1;
	return x;
}
int insert(num d){
	sc=fs=0;
	int p=insert(rt,0,1,d);
	if(sc){
		C=0;
		dfs(sc);
		(sc!=rt?(l[fs]==sc?l:r)[fs]:rt)=build(1,C,sl,sr);
	}
	return p;
}
int mx[400010],pos[100010];
void modify(int p,int l,int r,int x){
	if(l==r){
		mx[x]=l;
		return;
	}
	int mid=(l+r)>>1;
	if(p<=mid)
		modify(p,l,mid,x<<1);
	else
		modify(p,mid+1,r,x<<1|1);
	l=mx[x<<1];
	r=mx[x<<1|1];
	mx[x]=rk[pos[l]]>=rk[pos[r]]?l:r;
}
int query(int L,int R,int l,int r,int x){
	if(L<=l&&r<=R)return mx[x];
	int mid=(l+r)>>1,res=0,t;
	if(L<=mid){
		t=query(L,R,l,mid,x<<1);
		if(rk[pos[t]]>rk[pos[res]])res=t;
	}
	if(mid<R){
		t=query(L,R,mid+1,r,x<<1|1);
		if(rk[pos[t]]>rk[pos[res]])res=t;
	}
	return res;
}
int main(){
	int n,m,i,l,r,k;
	char s[5];
	scanf("%d%d",&n,&m);
	rk[0]=-1;
	insert(num());
	for(i=1;i<=n;i++)pos[i]=1;
	for(i=1;i<=n;i++)modify(i,1,n,1);
	while(m--){
		scanf("%s%d%d",s,&l,&r);
		if(s[0]=='C'){
			scanf("%d",&k);
			pos[k]=insert(num(pos[l],pos[r]));
			modify(k,1,n,1);
		}else
			printf("%d\n",query(l,r,1,n,1));
	}
}

转载于:https://www.cnblogs.com/jefflyy/p/9286077.html