二维树状数组【信息技术】.pptVIP

二维树状数组 胡伟栋 树状数组 问题： 有一个数列{an}，给出一种操作序列，每次可以改变数列中的一个元素的值 要求动态维护，使得任何时刻都能用较快速度求出数列的部分和 树状数组 令t[i]=a[x-lowbit(i)..i]的和 则S[i]=t[i]+S[lowbit(i)] S[i]表示a[i]的部分和 S[i]=sigma{t[x(i)k]} x(i)0=i x(i)1=x(i)0-lowbit(x(i)0) x(i)2=x(i)1-lowbit(x(i)1) …… 树状数组 当更新a[i]时，需要更新 t[y(i)0] y(i)0=i t[y(i)1] y(i)1=y(i)0+lowbit(y(i)0) t[y(i)2] y(i)2=y(i)1+lowbit(y(i)1) …… 树状数组 二维树状数组 问题 给出矩阵{aij} 可以随时改变矩阵某一个位置的值 要求能随时求出矩阵中某一个连续子矩阵的和 二维树状数组 此问题和原来的问题很像，只是由一维变成了二维 需要对原来的数据结构进行改进，以适应新的要求 二维树状数组 改进 用t[i,j]表示a[i-lowbit(i)+1..i, j-lowbit(j)+1..j]的和 S[i,j]=? 二维树状数组 计算S[15,15] 二维树状数组 S[i,j]= sigma{t[x(i)k, x(j)l]} 二维树状数组 更新a[i, j] 需要更新t[y(i)k, y(j)l] 二维树状数组 二维树状数组 程序实现 update(i, j, v) add(i, j, v-old[i, j]); add(i, j, v) j0?j while (imaxi) j?j0 while (jmaxj) t[i,j]+=v j?(j|(j-1))+1 i?(i|(i-1))+1 程序实现 sum(i, j) sum?0 j0?j while (i0) j?j0 while (j0) sum+=t[i,j] j=j-1 i=i-1 二维树状数组 复杂度 单个操作　O(log2max) 例 有一个线段盒子 A x y 表示插入一条[x,y]线段 D x y 表示删除一条[x,y]线段，如果当前盒子里没有，就什么也不做 S x y 表示要求出完全在[x,y]区间的线段数 命令数105 1=x,y=1000 * * 13 14 15 16 9 10 11 12 5 6 7 8 1 2 3 4 14 14 1 2 3 4 5 6 7 8 9 10 11 12 13 15 13 12 11 10 9 8 7 6 5 4 3 2 1 14 14 1 2 3 4 5 6 7 8 9 10 11 12 13 15 13 12 11 10 9 8 7 6 5 4 3 2 1 更新a[1,1]