可持久化treap

平衡树是一个有用的结构,但是常用的splay等由于均摊复杂度和旋转操作,无法方便的可持久化。

这时候跳出一个treap,他可以可持久化!

treap是一种笛卡尔树,有两个值:位置与权。这样构建出来的treap满足两个性质:1、左儿子位置小于当前节点小于右儿子,2、儿子的权值小于父亲。容易猜想,若权值随机,则树高是O(logn)。

使用类似二叉堆的旋转操作可以很方便的维护treap,但是还有一种更好的方法。

fhqtreap:平衡树中有两个操作非常关键:merge和split。其他大多数操作可以通过这两种操作实现。

(merge表示把两颗treap左右相接成一个新的treap,split表示把treap以某个元素为界分成两段)

不妨先考虑merge(a,b)。若val[a]>val[b],就递归合并rs[a]和b,否则合并a和ls[b]。

没了。。


如果需要可持久化,对每个操作建立新的节点,没了。。


#include
#define N 40005
#define M 4000005
#define ll long long
using namespace std;
int R,C,n,x,y,Rt;
int tg[M],ls[M],rs[M],size[M],val[M],cnt;ll sum[M],Ans;
vectorp[N];
int nd(int k=0)
{
	++cnt;
	tg[cnt]=tg[k];
	ls[cnt]=ls[k];
	rs[cnt]=rs[k];
	size[cnt]=size[k];
	val[cnt]=val[k];
	sum[cnt]=sum[k];
	return cnt;
}
void mdy(int k,int x)
{
	tg[k]+=x;
	val[k]+=x;
	sum[k]+=(ll)x*size[k]*(size[k]+1)/2;
}
void dn(int k)
{
	if (ls[k])mdy(ls[k],tg[k]);
	if (rs[k])mdy(rs[k],tg[k]);
	tg[k]=0;
}
void up(int k)
{
	size[k]=size[ls[k]]+size[rs[k]]+1;
	sum[k]=sum[ls[k]]+sum[rs[k]]+(ll)val[k]*size[k]*(size[k]+1)/2
	-(ll)val[k]*size[ls[k]]*(size[ls[k]]+1)/2
	-(ll)val[k]*size[rs[k]]*(size[rs[k]]+1)/2;
}
void build(int &k,int l,int r)
{
	if (l>r) {k=0;return;}
	k=nd();size[k]=r-l+1;
	int mid=l+r>>1;
	build(ls[k],l,mid-1);
	build(rs[k],mid+1,r);
}
void merge(int &k,int a,int b)
{
	if (!a){k=b;return;}
	if (!b){k=a;return;}
	dn(a);dn(b);
	if (val[a]


你可能感兴趣的:(算法)