替罪羊树（重量平衡树）入门
扯
学校清明竟然给放两天假期，心血来潮突然想去学习平衡树。
可是我太弱了学不会有旋转操作的treap和splay，这可怎么办啊qaq。
诶？我以前好像看过一种叫做替罪羊树的平衡树可以不用旋转操作，那还是学这个吧……
替罪羊树
用处
替罪羊树是一种平衡树，支持插入，删除，查找第k小元素，查找元素的排名等操作
什么数据结构优雅？暴力即是优雅！
替罪羊树就是一种暴力平衡树，旋转？不存在的！
替罪羊树保持平衡的方法就是暴力重构，即当树不平衡时拍扁重新建树，那么如何才能知道一棵树是否平衡呢？
在替罪羊树中选用了一种叫做平衡因子的东西，听起来很高端，其实就是一个0.5~1之间的任意浮点数，保持平衡的方法是这样的：
如果一棵树的左子树/右子树的存在的节点数量大于这棵树的存在的节点数量*旋转因子，那么就要重构这棵树
为什么我特意标出了是存在的节点数呢？是因为替罪羊树的删除不是真正的删除，而是惰性删除。
所以我们就可以写出代表替罪羊树的每个节点的结构体
1 const double alpha = 0.75; //旋转因子 2 struct Node { 3 Node* ch[2]; //左右子节点 4 int key,siz,cover; //key是值，siz是以该节点为根的树的存在的节点数，cover是所有节点数量 5 bool exist; //exist标志该节点是否被删除 6 void pushup() { //更新函数 7 this->siz=ch[0]->siz+ch[1]->siz+(int)exist; 8 this->cover=ch[0]->cover+ch[1]->cover+1; 9 } 10 int isbad() { //判断是否要重构 11 return (ch[0]->cover>this->cover*alpha+5)||(ch[1]->cover>this->cover*alpha+5); 12 } 13 };
（这种东西还是很好理解的吧，几乎所有平衡树都要写这种东西）
内部操作
内部操作指实现操作的函数需要调用的子函数，这些函数不会在主程序中调用
总的来说就只有四种内部操作：新节点，插入，删除，重构
新节点有什么可讲的？不就是一个new Node就行的事吗？
不不不，动态分配内存的速度可以说是非常慢了，所以我们要手写内存池。
1 Node mempol[maxn]; //内存池 2 Node *tail; //tail为指向内存池元素的指针 3 Node *bc[maxn]; //内存回收池（栈） 4 int bc_top; //内存回收池（栈）顶指针
分配机制很简单，如果*bc为空，那么从mempol里获得一个节点的内存，当删除节点时，把节点存入*bc，需要时再从*bc中取出即可。
另外还有两个不可缺少的东西，根节点和null节点，为Node*类型。
（null是为了代替NULL的，看到错误就像砸电脑，所以就自己模拟一个空节点）
有了这些东西，新节点操作也就不难写出：
1 Node* newnode(int key) { //返回一个新节点 2 Node* p=bc_top?bc[--bc_top]:tail++; //分配内存 3 p->ch[0]=p->ch[1]=null; 4 p->cover=p->siz=p->exist=1; 5 p->key=key; 6 return p; 7 }
接下来就是插入操作了，插入操作和二叉搜索树是一样的，如果当前节点小于等于要插入的节点权值，那么递归操作左子树，反之递归操作右子树
同时维护siz和cover
1 Node** insert(Node*& p,int val) { //返回指向距离根节点最近的一棵不平衡的子树的指针 2 if(p==null) { 3 p=newnode(val); 4 return &null; 5 } else { 6 p->siz++,p->cover++; //维护节点数 7 Node** res=insert(p->ch[val>=p->key],val); 8 if(p->isbad()) res=&p; 9 return res; 10 } 11 }
这里插入也是有返回值的，返回的是指向距离根节点最近的一棵不平衡的子树的指针，因为我们定义的树是指针实现的，所以这个地方要用到指向指针的指针，即双重指针。那么返回这么个拗口的东西又有什么用呢？前面有说过重构平衡树的方法，找到距离根节点最近的不平衡的子树就是为了能够一遍将树重构成平衡树，而不用多遍。
删除比较简单，因为不用像其他树一样真正的删除然后旋转，替罪羊树的删除操作只用将exist置为false然后更新siz值即可
还有一点就是这里的删除操作并不是删除值为k的节点，而是删除第k小的节点。
求第k小的节点要用到siz，根据二叉搜索树的性质，我们可以知道一个节点的左子树的值都要小于这个节点，所以当k小于左子树的节点时递归操作左子树，反之递归操作右子树，同时减去左子树的节点数即可
1 void erase(Node*& p,int k) { 2 p->siz--; //维护siz 3 int offset=p->ch[0]->siz+p->exist; //计算左子树的存在的节点总数 4 if(p->exist&&k==offset) { //判断当前节点权值是否第k小 5 p->exist=false; //删除节点 6 } else { 7 if(k<=offset) erase(p->ch[0],k); //如果k小于等于offset，递归操作左子树 8 else erase(p->ch[1],k-offset); //反之递归操作右子树，不要忘记将k减去offset 9 } 10 }
插入和删除都说完了，接下来是有些难度的重构操作。
重构需要用到三个函数，第一个是将树转化成有序序列的函数，第二个是建树的函数，而第三个就是将前两个函数连接在一起的函数
转化成序列的函数是怎么做到有序的呢？别忘了替罪羊树也是一颗二叉搜索树，所以我们可以跑一边中序遍历，把结果存入vector中，序列就一定有序
前面说过的删除标记exist在这里起到了作用，把树拍扁的同时判断节点是否存在，如果不存在直接扔到内存回收池中就好
所以说替罪羊树的删除操作只在重构操作中执行，erase函数只是惰性删除而已
1 void travel(Node* p,vector<Node*>& x) { //将一棵树转化成序列，保存在vector中 2 if(p==null) return; //如果是空树则退出 3 travel(p->ch[0],x); //递归操作左子树 4 if(p->exist) x.push_back(p); //如果该节点存在则放入序列中 5 else bc[bc_top++]=p; //回收内存，将不用的节点扔到内存回收池（栈）中 6 travel(p->ch[1],x); //递归操作右子树 7 }
建树函数就不用说啦，这个是基本操作吧，取序列的中间值作为树的根，然后递归操作左子树和右子树即可，显然能保证平衡
1 Node* divide(vector<Node*>& x,int l,int r) { //返回建好的树 2 if(l>=r) return null; //序列为空不用建树 3 int mid=(l+r)>>1; 4 Node* p=x[mid]; //mid保证平衡 5 p->ch[0]=divide(x,l,mid); //递归操作 6 p->ch[1]=divide(x,mid+1,r); //递归操作 7 p->pushup(); //维护节点信息 8 return p; 9 }
接下来用一个函数把前两个函数连接，就是一个完整的重构函数了
1 void rebuild(Node*& p) { 2 static vector<Node*> v; 3 v.clear(); 4 travel(p,v); //拍扁 5 p=divide(v,0,v.size()); //建树 6 }
至此所有内部函数已经全部实现，替罪羊树的全部操作函数可以由以上函数完成。
操作函数
哈哈，想不到吧，接下来才是实现替罪羊树操作的真正函数
有了上面的内部函数，操作函数的实现就显得非常简单
初始化
首先是初始化函数，其实没有什么可以初始化的，只是定义了一下空节点和根节点而已
把初始化函数放到构造函数中可以自动调用
1 void init() { 2 tail=mempol; //tail指向内存池的第一个元素 3 null=tail++; //为null指针分配内存 4 null->ch[0]=null->ch[1]=null; //null的两个儿子也是null 5 null->cover=null->siz=null->key=0; //null的所有标记都是0 6 root=null; //初始化根节点 7 bc_top=0; //清空栈 8 } 9 10 STree() { 11 init(); 12 }
插入
插入函数的实现只有两行，调用内部插入函数，然后判断是否需要重建，如果需要调用重建函数即可
1 void insert(int val) { 2 Node** res=insert(root,val); 3 if(*res!=null) rebuild(*res); 4 }
求值为val的节点的名次
和内部删除的算法差不多，但是因为只是求rank而不是删除，所以没必要用递归函数，用循环实现即可
1 int rank(int val) { 2 Node* now=root; 3 int ans=1; 4 while(now!=null) { 5 if(now->key>=val) now=now->ch[0]; 6 else { 7 ans+=now->ch[0]->siz+now->exist; 8 now=now->ch[1]; 9 } 10 } 11 return ans; 12 }
求第k小元素
由判断值变成判断名次，根据siz判断即可，思想和求rank差不多，用循环可以解决
1 int kth(int val) { 2 Node* now=root; 3 while(now!=null) { 4 if(now->ch[0]->siz+1==val&&now->exist) return now->key; 5 else if(now->ch[0]->siz>=val) now=now->ch[0]; 6 else val-=now->ch[0]->siz+now->exist,now=now->ch[1]; 7 } 8 }
删除
删除操作分两种，删除第k小和删除值为val的元素，但不是只删除就行的。
为了不让无用的节点过多，我们可以在有用节点与全部节点的比值小于平衡因子的的情况时重构整棵树，以提高效率。
1 void erase(int k) { //删除值为k的元素 2 erase(root,rank(k)); 3 if(root->siz<root->cover*alpha) rebuild(root); 4 } 5 6 void erase_kth(int k) { //删除第k小 7 erase(root,k); 8 if(root->siz<root->cover*alpha) rebuild(root); 9 }
替罪羊树的全部操作到此已经结束，其他操作就可以根据题目自己yy了。
模板
板子题：洛谷 P3369 【模板】普通平衡树（Treap/SBT)
其实把上面所有代码连成一串就是板子了qaq
代码：
1 #include <cstdio> 2 #include <vector> 3 using std::vector; 4 5 namespace Scapegoat_Tree { 6 const int maxn = 100000 + 10; 7 const double alpha = 0.75; //旋转因子 8 struct Node { 9 Node* ch[2]; //左右子节点 10 int key,siz,cover; //key是值，siz是以该节点为根的树的存在的节点数，cover是所有节点数量 11 bool exist; //exist标志该节点是否被删除 12 void pushup() { //更新函数 13 this->siz=ch[0]->siz+ch[1]->siz+(int)exist; 14 this->cover=ch[0]->cover+ch[1]->cover+1; 15 } 16 int isbad() { //判断是否要重构 17 return (ch[0]->cover>this->cover*alpha+5)||(ch[1]->cover>this->cover*alpha+5); 18 } 19 }; 20 struct STree { 21 protected: 22 Node mempol[maxn]; //内存池 23 Node *tail,*null,*root; //tail为指向内存池元素的指针 24 Node *bc[maxn]; //内存回收池（栈） 25 int bc_top; //内存回收池（栈）顶指针 26 27 Node* newnode(int key) { 28 Node* p=bc_top?bc[--bc_top]:tail++; 29 p->ch[0]=p->ch[1]=null; 30 p->cover=p->siz=p->exist=1; 31 p->key=key; 32 return p; 33 } 34 35 void travel(Node* p,vector<Node*>& x) { //将一棵树转化成序列，保存在vector中 36 if(p==null) return; //如果是空树则退出 37 travel(p->ch[0],x); //递归操作左子树 38 if(p->exist) x.push_back(p); //如果该节点存在则放入序列中 39 else bc[bc_top++]=p; //回收内存，将不用的节点扔到内存回收池（栈）中 40 travel(p->ch[1],x); //递归操作右子树 41 } 42 43 Node* divide(vector<Node*>& x,int l,int r) { //返回建好的树 44 if(l>=r) return null; //序列为空不用建树 45 int mid=(l+r)>>1; 46 Node* p=x[mid]; //mid保证平衡 47 p->ch[0]=divide(x,l,mid); //递归操作 48 p->ch[1]=divide(x,mid+1,r); //递归操作 49 p->pushup(); //维护节点信息 50 return p; 51 } 52 53 void rebuild(Node*& p) { 54 static vector<Node*> v; 55 v.clear(); 56 travel(p,v); //拍扁 57 p=divide(v,0,v.size()); //建树 58 } 59 60 Node** insert(Node*& p,int val) { //返回指向距离根节点最近的一棵不平衡的子树的指针 61 if(p==null) { 62 p=newnode(val); 63 return &null; 64 } else { 65 p->siz++,p->cover++; //维护节点数 66 Node** res=insert(p->ch[val>=p->key],val); 67 if(p->isbad()) res=&p; 68 return res; 69 } 70 } 71 72 void erase(Node*& p,int k) { 73 p->siz--; //维护siz 74 int offset=p->ch[0]->siz+p->exist; //计算左子树的存在的节点总数 75 if(p->exist&&k==offset) { //判断当前节点权值是否第k小 76 p->exist=false; //删除节点 77 } else { 78 if(k<=offset) erase(p->ch[0],k); //如果k小于等于offset，递归操作左子树 79 else erase(p->ch[1],k-offset); //反之递归操作右子树 80 } 81 } 82 83 void iod(Node* p) { 84 if(p!=null) { 85 iod(p->ch[0]); 86 printf("%d ",p->key); 87 iod(p->ch[1]); 88 } 89 } 90 public: 91 void init() { 92 tail=mempol; //tail指向内存池的第一个元素 93 null=tail++; //为null指针分配内存 94 null->ch[0]=null->ch[1]=null; //null的两个儿子也是null 95 null->cover=null->siz=null->key=0; //null的所有标记都是0 96 root=null; //初始化根节点 97 bc_top=0; //清空栈 98 } 99 100 STree() { 101 init(); 102 } 103 104 void insert(int val) { 105 Node** res=insert(root,val); 106 if(*res!=null) rebuild(*res); 107 } 108 109 int rank(int val) { 110 Node* now=root; 111 int ans=1; 112 while(now!=null) { 113 if(now->key>=val) now=now->ch[0]; 114 else { 115 ans+=now->ch[0]->siz+now->exist; 116 now=now->ch[1]; 117 } 118 } 119 return ans; 120 } 121 122 int kth(int val) { 123 Node* now=root; 124 while(now!=null) { 125 if(now->ch[0]->siz+1==val&&now->exist) return now->key; 126 else if(now->ch[0]->siz>=val) now=now->ch[0]; 127 else val-=now->ch[0]->siz+now->exist,now=now->ch[1]; 128 } 129 } 130 131 void erase(int k) { //删除值为k的元素 132 erase(root,rank(k)); 133 if(root->siz<root->cover*alpha) rebuild(root); 134 } 135 136 void erase_kth(int k) { //删除第k小 137 erase(root,k); 138 if(root->siz<root->cover*alpha) rebuild(root); 139 } 140 141 void iod() { //调试用的中序遍历 142 Node* p=root; 143 iod(p); 144 } 145 }; 146 } 147 148 using namespace Scapegoat_Tree; 149 STree st; 150 int main(int argc,char** argv) { 151 int n,opt,que; 152 scanf("%d",&n); 153 while(n--) { 154 scanf("%d%d",&opt,&que); 155 if(opt==1) st.insert(que); 156 if(opt==2) st.erase(que); 157 if(opt==3) printf("%d\n",st.rank(que)); 158 if(opt==4) printf("%d\n",st.kth(que)); 159 if(opt==5) printf("%d\n",st.kth(st.rank(que)-1)); 160 if(opt==6) printf("%d\n",st.kth(st.rank(que+1))); 161 if(opt==7) st.iod(); 162 } 163 return 0; 164 }
--完--
