树状数组功能及实现介绍

树状数组的概念:
   树状数组(Fenwick Tree,又称二叉索引树)是一个查询和修改复杂度都为log(n)的数据结构,它是利用二进制的一些特点来实现。它的功能有局限性,主要是用来动态查询连续和(或者是前缀和)的问题。它利用O(n)的附加空间复杂度,将线性的数列结构转化成树状结构从而进行跨越扫描,高效完成查询连续和。

看图来推导一下他们的对应关系:

  C2=A1+A2   C4=C2+A3+A4=A1+A2+A3+A4   C6=A5+A6.......

  这里就不一一列出所有对应关系了,读者可以非常明白的看到是如何对应的。

经典例题与分析:
  那么这里给出一例:给出一个n个元素的数组A1,A2,A3...An,要求对他们进行两种操作: 1、Add操作:Add(x,d)让Ax加上d。2、Query操作:Query(L,R)计算Al+A(l+1)+...+Ar。

  这是一个很典型的动态查询区间和问题,如果按朴素前缀和递推的算法,那么每次修改一个参数后都要重新计算后面所有的前缀和数组S的值,修改m次值,每次递推重新计算前缀和数组S,时间复杂度为T(n)=O(m*n)(m为执行指令Add的次数),那么当m、n均大于10000时将无法承受时间复杂度,那么就要靠树状数组了!

 在学习树状数组前,要先介绍一个新函数lowbit(x),他的作用是将x转换为二进制之后从右往左第一个1的位置。这个函数很重要,是建树的索引,那么如何实现这样一个函数呢?给出公式:lowbit ( x ) = x & ( -x )

为什么x&(-x)可以实现这个功能?在计算机中负号具体运算规则为:按位取反(包括首位表示符号的一位),末尾加1之后的结果,还是给读者举上一例:8转化为二进制:0 1000  取反:1 0111  加一:1 1000 转回十进制:(符号1为负数)-8。这是‘-’取反的意义,那么‘&‘意义是将两个数字二进制进行比较,相同位数字都为则为1,否则为0,那么最后他们取&,算出的不就是从右往左1最先出现的位置了么。这个函数介绍完了,就可进入正题了,请读者观察上图联系刚刚介绍的lowbit函数,有什么发现呢?无论是在叶子上,还是在中间的节点上,对于左子节点,它的父亲节点编号是i+lowbit(i),对于右子节点,它的父亲节点编号是i-lowbit(i),而且还有一个很重要的发现,那就是C数组的覆盖范围:C[x]的区间是[x-lowbit(x)+1,x],读者可以自行验算。

  有了以上理论支持,问题便不难分析了,我要计算A1+A2+A3+...+A11,也就是A11的前缀和,如何计算?从11开始向左边走,用C对应的数组不重复也不漏地走到A1,所加起来的和便是所要的答案,具体可以见上图中红色标记,A11的前缀和=A11+C10+C8,明白了么?

  有了计算的方法,那么修改一个叶子上的值会变化哪些值呢?只需同刚才使用的方法,将上面的根节点都加上d即可完成修改。也就是说修改的值就很少了,时间复杂度T(n)=O(logn)。

  本题还有一个问题,非常简单,就是要求计算的是闭合区间[L,R]的和,而树状数组算出的却是前缀和,那如何才能转换过来呢?数学中学过的集合,用A表示Ar的前缀和,用B表示Al的前缀和,则B⊆A&&A!=B,那么[L,R]=A为全集的B的补集,用pre_sum(x)表示前缀和,则闭合区间[L,R]中元素的和为pre_sum[R]-pre_sum[L-1]。

  接下来上代码:

#include <bits/stdc++.h>
using namespace std;
const int MAXN=500005;
int C[MAXN],N,M,num;
int lowbit(int v){return v&(-v);}
void Add(int x,int value){
    while(x<=N){
        C[x]+=value;
        x+=lowbit(x);
    }
    return;
}
int pre_Query(int k){
    int sum=0;
    while(k>0){
        sum+=C[k];k-=lowbit(k);
    }
    return sum;
}
int pre_sum(int L,int R){
    return pre_Query(R)-pre_Query(L-1);
}
int main(){
    cin>>N>>M;
    for(int i=1;i<=N;i++){
        cin>>num;
        Add(i,num);
    }
    for(int i=1;i<=M;i++){
        int code,par1,par2;
        scanf("%d",&code);
        if(code==1){scanf("%d%d",&par1,&par2);Add(par1,par2);}
        if(code==2){
            scanf("%d%d",&par1,&par2);
            printf("%d
",pre_sum(par1,par2));
        }
    }
    return 0;
}

  

原文地址:https://www.cnblogs.com/xingchenOI/p/10424315.html