利用树状数组解决几类问题讲述.doc

利用树状数组解决几类问题 浙江省温岭中学 林希 树状数组作为一种实现简单、应用较广的高级数据结构，在OI界的位越来越重要，下面我来简单介绍一下树状数组和它的简单应用。 一、树状数组简介 树状数组(Binary Indexed Trees，简称BIT)可以方便地查询出一段区间中的数字之和。其查询和修改的时间复杂度均为O(N),并且是一个在线的数据结构，可以随时修改并查询。我接下来介绍树状数组的几个基本操作（s，t都是输入的数，s=t） 这道题目用线段树很容易可以解决，我就不多说了，那么如何用树状数组来解决呢？我们新增一个数组C[]，其中C[i]=a[i-2k+1]+……+a[i](k为i在二进制形式下末尾0的个数)。根据定义我们可以得出一下这张表格： i 二进制 K 1 (1)2 0 c[1]=a[1] 2 (10) 2 1 c[2]=a[1]+a[2]=c[1]+a[2] 3 (11) 2 0 c[3]=a[3] 4 (100) 2 2 c[4]=a[1]+a[2]+a[3]+a[4]=c[2]+c[3]+a[4] 5 (101) 2 0 c[5]=a[5] 6 (110) 2 1 c[6]=a[5]+a[6]=c[5]+a[6] …… 在这里，我们会发现k值的求解会有一些难度，这就引出了树状数组的第一个操作：低位技术Lowbit(k)。 对于Lowbit这里我提三种不同的写法1.Lowbit(k)=k and (k or (k-1)) 2.Lowbit(k)=k and not (k-1) 3.Lowbit(x)=k and –k 然后我来分析引例的树状数组解法，为了可以更好地理解这种方法，读者可以根据以下这幅图来加以理解。 【操作2】修改操作：可以从C[i]往根节点一路上溯，调整这条路上的所有C[]即可，这个操作的复杂度在最坏情况下就是树的高度即O(lb)。?……C[pm]， 则p1=k，而（li为pi在二进制中末尾０的个数）。 ?例如A[1]……A[8]中，a[3] 添加x； p1=k=3 p2=3+20=4 p3=4+22=8 p4=8+23=168 由此得出，c[3]、c[4]、c[8]亦应该添加x。 定理的证明如下： 【引 理】若A[k]所牵动的序列为C[p1]，C[p2] C[pm]， 且p1 p2…pm，则有l1 l2…lm（li为pi在二进制中末尾０的个数）。 证明：若存在某个i有li≥li+1，则，，即： (1) 而由Li=Li+1、Pi Pi+1可得 (2) (2)矛盾，可知l1 l2 ……lm 定理2：p1=k，而p i+1=pi+2li　 证明：p1 p2 …… pm　且C[p1]，C[p2] C[pm]中包含a[k]，因此p1=k。在p序列中，pi+1=pi +2li是pi后最小的一个满足li+1 li的数（若出现Pi+x比pi+1更小，则x2li ，与x在二进制中的位数小于li相矛盾）。Pi+1=pi+2li，li+1≥li+1。 由pi-2li+1≤K≤Pi可知，Pi+1-2li+1+1≤Pi+2li–2*2li+1=Pi– 2li+1≤K≤Pi≤P i+1 ，故Pi与pi+1之间的递推关系式为 P i+1=Pi+2li 【操作3】求数列的前n项和只需找到n以前的所有最大子树，把其根节点的C加起来即可。不难发现，这些子树的数目是n在二进制时1的个数，或者说是把n展开成2的幂方和时的项数,?因此，求和操作的复杂度也是O(l)。… +a[k] (l为k在二进制数中末尾０的个数)，我们从k1=k出发，按：ki+1=ki-2lki（lki为ki在二进制数中末尾0的个数） 递推k2，k3，…，km (km+1=0)。 由此得出 S=c[k1]+c[k2]+c[k3] + … + c[km] 相信看了这两个操作的理论部分，读者应该有了一定的理解，下面给出这两种操作和其他一些重要操作的Pascal代码。 【操作1】低位技术：Lowbit(k) Lowbit(k)：即2进制数中从最低位开始连续0的位数的关于2的幂(k:Longint) :Longint; Begin Lowbit:=k and –k; End; 【操作2】修改操作：Modify(k, D) Modify(k, D)：对数组 a[] 中的第个元素加上。为了维护 [] 数组，我就必须要把 [] 中所有管着 a[] 的 c[] 全部加上，这样才能随时以 O(l) 的复杂度进行 um(i) 的操作。而 owbit(k) 正是依次访问所有包含 a[] 的 c[] 的