线段树简介:
线段树是一种高效的数据结构,他的应用非常广泛,主要可以求区间和,求区间最小值,求区间最大值,求区间中某个数据出现的次数……为什么它效率这么高呢?因为它用到了二分的思想,将一个完整的线段的区间答案保存在根节点,左边是线段的左半部分,右边是线段的右半部分,以此类推。如下面就是一个线段树(求区间和):
线段树的性质:
由上图,我们可以发现以下性质:1.一个由N个数据构成的线段树,其结点的数量为2N-1(不信你画)。
2.一个由N个数据构成的线段树,其树的深度(以根节点深度为1算)为log2N(向上取整)+1(不信你画)。
3.一个节点所管理的区间为L~R,那么它左孩子管理范围是L~(L+R)/2,右孩子是(L+R)/2+1~R(人家就这样定义
的)
4.每个线段树的叶子节点都是他的原有的数据(自己看图)。
5.线段树是一个完全二叉树。
线段树的基本操作:
线段树的操作主要有三种:建树,查询,修改。其中修改又能细分为修改单个节点和{【修改整条线段】<-主要要讨论的内容}
线段树的建造:
首先我们都知道,树是递归定义的,那我们仔细观察观察上面那棵线段树,我们可以发现什么? 没错,对一棵线段树按层序遍历编号:对于一个父亲节点而言,它的两个孩子节点之和就是它自己的值。我们可以推广这个结论,对一棵关于x操作的线段树,它的非叶子节点i的值等于他的孩子节点i*2和i*2+1通过x操作所得到的值,举个例子:设一个x操作为求区间最小值的线段树,它的第i*2节点的值为25,它的第i*2+1节点的值为28,那么第i个节点的值就是min(25,28) = 25;又比如上面那棵求和树,它的第二个节点的值(105)=它的第2*2=4个节点的值(45)+它的第2*2+1个节点的值(60);
由性质5知,线段树是一个完全二叉树,那么我们就可以用数组来保存它(像保存堆那样):
1 struct node{ 2 int value; 3 }tree[10000 + 20];
这样,我们就不难写出建树操作了,下面以求和树为栗子:
1 void make_tree(int begin, int end, int pla){ //构造线段树的过程,pla为当前处理数组的哪个位置 2 if(begin == end){ //叶子节点 3 tree[pla].value = a[begin]; 4 return; 5 }else{ //递归建树 6 make_tree(begin, (begin + end) / 2, pla * 2); 7 make_tree((begin + end) / 2 + 1, end, pla * 2 + 1); 8 } 9 tree[pla].value = tree[pla * 2].value + tree[pla * 2 + 1].value; //求和 10 }
从上面这个代码我们可以看出,建树操作的复杂度是O(n)。
线段树的查找:
关于线段树的查找,各位应该看到树的样式就差不多明白了一些,比如上面那个求和树,查找1-7应该会吧,1-4应该也会吧,只要沿着树二分查找就行了。但是,看到这里各位肯定会产生到一丝疑惑,怎么查找不在这个线段树节点区间上的值,比如1-6,2-5,2-6?
我们不妨这样看:如下图,将1-6看成1-4和5-6,这样我们只需要进行两次查找,然后再进行一次那样的操作,就能得到答案了,比如上面那棵线段树,1-4为105,5-6为81,然后进行一次相加操作,得到答案186,不信你看图。
(图:查找线段1-6)
以此类推,2-5可以看成2,3-4,5,根据上面那个二叉树,2为9,3-4为60,5为11,相加80,自己验证看看对不对。
2-6就略了,你们可以自己算算。
那么我们如何具体实现这个操作呢,我们不妨用集合的思想来想想:
当我们从根节点向下查找时,我们可以设计一个比较,比较当前要查找的这个节点begin~end和查找范围kaishi~jieshu的关系,我们设A={x | x在begin~end所代表的元素之间},B={x | x在kaishi~jieshu所代表的元素之间}。如当begin=1,end=4时,A={36,9,31,29},当kaishi=2,jieshu=5时,B={9,31,29,11};
这样,我们就可以发现:当A∩B≠∅时,也就是说A这个范围内有我想要的东西,那我就进一步查找这个区间内的左右子树。
当A∩B=∅时,返回。
当A⊆B时,说明我们不需要进一步查找了,就可以存下来这个值,等到完成操作时再进行操作,当然我们也可以即时操作。
那么代码应该也不难写了,下面是查找操作的代码:
1 int sum = 0; 2 void lookup(int begin,int end,int kaishi,int jieshu,int pla){ 3 if(jieshu < begin || kaishi > end) //表示线段集合begin~end交线段集合kaishi~jieshu为空集合的情况 4 return; 5 else{ 6 if(begin >= kaishi && end <= jieshu){ //线段集合begin~end子集于线段集合kaishi~jieshu 7 sum += tree[pla].value; 8 }else{ //交集不为空 9 lookup(begin, (begin + end) / 2, kaishi, jieshu , pla * 2); 10 lookup((begin + end) / 2 + 1, end, kaishi, jieshu, pla * 2 + 1); 11 } 12 } 13 }
由上面的代码,我们不难发现,当存在着一棵2N-1个节点的二叉树,最坏情况为查找2~N-1(自己去证)。那么它的时间复杂度为O(log2N);(你可以试着用笔画一下,会发现经过节点左右对称,数量大概是2log2N,所以是O(log2N));
线段树的修改操作:
线段树的修改分为两部分,一个是单个节点的修改,一个是线段的修改,相信大家看完线段树的查找后一定有了灵感,单个节点的修改肯定是没问题了。但是,如果按照单个节点的修改方法修改线段,假设线段长度为K,那么它的时间复杂度就是O(Klog2N),这个效率是很差的,因为我们有方法可以使时间优化到O(log2N),这就要使用二叉树的最难的东西了——延迟标记;
那么我们要改一下代表线段树的数组的定义:
1 struct node{ 2 int value; 3 int biaoji; 4 node(){biaoji = 0;} 5 }tree[10000 + 20];
延迟标记,顾名思义,就是延迟某项东西的处理时间,换句话说就是将一些事情拖到后面去做,到底是将什么拖到后面去做呢?既然是修改节点,那么肯定就是把修改节点的操作拖到后面去做。看到这里有些人就会疑惑了,既然是要修改线段,那么为什么不立即修改了呢,不修改那还叫什么修改线段呢?其实,我的意思是说不修改线段上的每一个具体的点,而是将树上某一些代表某条线段的节点给修改,然后存到延迟标记中,查找的时候消除掉它,应该说是下移。
说到这里肯定有人没理解,因为我也是看了好多次才懂的,举个实例吧,就给上面那棵线段树中1-4这一段每个节点增加3吧。这时,我们就先二分查找到代表1-4的2号节点,这时按照修改单个节点的方法就是继续向下查找,修改1,2,3,4,然后递归修改它们的父节点,这样的话效率就是O(4log27),当要修改多个线段时这种时间复杂度是无法承受的,但我们引入延迟标记后,我们可以直接修改2号节点的值,并不向下继续修改1,2,3,4的值,然后更改2号节点的延迟标记增加3,意思是说1-4这个区间每个数我都欠了3给他们加。这时候肯定有人会问查找怎么办?我们仔细观察一下查找的步骤,都是由根往叶子查找的,而且查找没修改的点肯定会遇到被修改的那个点,比如说上面给1-4增加了3,那么我们要查找1-2的值时,必然会路过1-4,这时候延迟标记的作用又来了,我们将延迟标记向下移,将原先属于2号节点(1-4)的延迟标记送给它的两个孩子,即使它两个孩子的延迟标记增加3,然后给它的左右孩子的值增加(R-L+1)*3点。这样使得查找出来的结果正确,如果查找的范围不正好是节点所表示的区间的话,就用查找的方法二分。
代码如下:
1 void change(long long begin,long long end,long long kaishi,long long jieshu,long long pla,long long add) { 2 if(jieshu < begin || kaishi > end) 3 return; 4 else{ 5 if(begin >= kaishi && end <= jieshu){ 6 tree[pla].biaoji += add; 7 tree[pla].value += add * (end - begin + 1); 8 }else{ 9 if(tree[pla].biaoji){ 10 tree[pla * 2].biaoji += tree[pla].biaoji; 11 tree[pla * 2].value += ((begin + end) / 2 - begin + 1) * tree[pla].biaoji; 12 tree[pla * 2 + 1].biaoji += tree[pla].biaoji; 13 tree[pla * 2 + 1].value += (end - (begin+end) / 2 ) * tree[pla].biaoji; 14 tree[pla].biaoji=0; 15 } 16 change(begin, (begin + end) / 2, kaishi, jieshu , pla * 2,add); 17 change((begin + end) / 2 + 1, end, kaishi, jieshu, pla * 2 + 1,add); 18 if(begin!=end){ 19 tree[pla].value = tree[pla*2].value + tree[pla*2+1].value; 20 } 21 } 22 } 23 }
反复读这个代码你就会知道到底延迟标记有什么用了。最坏情况O(log2n),证明方法很简单,画张图就行了,画一张1~n的图,然后选择改变2~n-1,改变的是每层两个。一个二叉树有log2n层,所以复杂度为O(log2n);
上面我们说过,查找的方式就是将延迟标记下移,这是坠吼的,所以我们要对查找进行一点PY交易。
1 long long sum = 0; 2 void lookup(long long begin,long long end,long long kaishi,long long jieshu,long long pla){ 3 if(jieshu < begin || kaishi > end) 4 return; 5 else{ 6 if(begin >= kaishi && end <= jieshu){ 7 sum += tree[pla].value; 8 }else{ 9 if(tree[pla].biaoji){ 10 tree[pla * 2].biaoji += tree[pla].biaoji; 11 tree[pla * 2].value += ((begin + end) / 2 - begin + 1) * tree[pla].biaoji; 12 tree[pla * 2 + 1].biaoji += tree[pla].biaoji; 13 tree[pla * 2 + 1].value += (end - (begin+end) / 2 ) * tree[pla].biaoji; 14 tree[pla].biaoji=0; 15 } 16 lookup(begin, (begin + end) / 2, kaishi, jieshu , pla * 2); 17 lookup((begin + end) / 2 + 1, end, kaishi,jieshu, pla * 2 + 1); 18 } 19 } 20 }
进行完整棵肮脏的交易后,线段树就打完了,后面的两个延迟标记的我没测试,请看了的帮忙测试下,非常感谢!