码迷,mamicode.com
首页 > 其他好文 > 详细

树状数组

时间:2014-06-16 07:24:03      阅读:190      评论:0      收藏:0      [点我收藏+]

标签:style   class   blog   code   java   http   

一、树状数组

树状数组是一个优美小巧的数据结构,在很多时候可以代替线段树。一句话概括就是,凡是树状数组可以解决的问题,线段树都可以解决,反过来线段树可以解决的问题,树状数组不一定能解决。

树状数组英文名称为Binary Index Tree,直译过来就是二进制索引树,我觉得二进制索引树更能说明其本质。树状数组的本质就是一种通过二进制位来维护一个序列前i和的数据结构。

bubuko.com,布布扣

 

1、更新元素

比如修改了a[3],则:

P1 = x = 3, L1 = 0;

P2 = p1+ 2(L1) = 3 + 1 = 4, L2 = 2;

P3 = p2 + 2(L2) = 4 + 4 = 8, L3 = 3;

....

2、子序列求和

比如求sum(6)

p1 = x = 6, L1 = 1; sum += c[6];

p2 = p1 - 2(L1) = 6 - 2 = 4; L2 = 2; sum += c[4];

p3 = p2 - 2(L2) = 4 - 4 = 0, return;

bubuko.com,布布扣

二、树状数组的维护

对于维护的序列A,定义C[i]=A[j+1]+...+A[i],其中ji的二进制表示中把最右边的1换成0的值。j的值可以通过lowbit求出,即i-lowbit(i)

lowbit(a)2^(a的二进制表示末尾0的个数)。可以用下面式子求出:

lowbit(a)=a&(a ^ (a-1))
//
lowbit(a)=a&(~a+1)
//或根据补码性质
lowbit(a)=a&(-a)

 

 修改方式如下

bubuko.com,布布扣
    void modify(int p,int delta)
    {
        while (p<=N)
        {
            C[p]+=delta;
            p+=lowbit(p);
        }
    }
bubuko.com,布布扣

 

求前缀和如下

bubuko.com,布布扣
    int sum(int p)
    {
        int rs=0;
        while (p)
        {
            rs+=C[p];
            p-=lowbit(p);
        }
        return rs;
    }
bubuko.com,布布扣

 

三、示例代码

bubuko.com,布布扣
package ads;

/**
 * 
 * @author loull
 * 下标index从1开始数
 */
public class BinaryIndexTree {

    int[] c;
    int n;
    public BinaryIndexTree(int n) {
        this.n = n;
        c = new int[n+1];
    }
    
    int lowbit(int x) {
        return x & (x ^ (x-1));
    }
    /**
     * 更新一个元素,初始数组c都是0,所以可以用这个方法初始化c,这时候增加与更新是等价的
     * @param p    更新第p个元素,索引从1开始
     * @param d    增加的值,不是更新后的值
     */
    void update(int p, int d) {
        while (p <= n) {
            c[p] += d;
            p += lowbit(p);
        }
    }
    
    /**
     * 计算从第一个元素到第p个元素的和
     * @param p
     * @return
     */
    int sum(int p) {
        int ret = 0;
        while (p > 0) {
            ret += c[p];
            p -= lowbit(p);
        }
        return ret;
    }
    
    public static void main(String[] args) {
        int[] numbers = {1,2,3,4,5,6,7,8,9};
        BinaryIndexTree bit = new BinaryIndexTree(numbers.length);
        for (int i=0; i<numbers.length; i++) {
            bit.update(i+1, numbers[i]);
        }
        
        System.out.println(bit.sum(9));
    }
}
bubuko.com,布布扣

 

 

四、C++模板

bubuko.com,布布扣
/***************************************
*                                      *
*    Template of Binary Indexed Tree   *
*                                      *
*           @author: Zyzsdy            *
*                                      *
****************************************
*   Thanks: shu_mj                     *
*           ACM Template of Tokyo U    *
***************************************/
template <class Type>
struct binaryIndexTree {
    Type *data;
    int size_n;
    binaryIndexTree (int _size) : size_n (_size + 1) {
        data = new Type[size_n]();
    }
    ~binaryIndexTree() {
        delete [] data;
    }
    void add (int _lower, Type _delta) {
        int pos = _lower + 1;
        while (pos < size_n) {
            data[pos] += _delta;
            pos += pos & -pos;
        }
    }
    Type sum (int _lower, int _upper) {
        if (_lower > 0) {
            return sum (0, _upper) - sum (0, _lower);
        } else {
            Type result = 0;
            int pos = _upper;
            while (pos > 0) {
                result += data[pos];
                pos -= pos & -pos;
            }
            return result;
        }
    }
};
typedef binaryIndexTree<int> BIT;
bubuko.com,布布扣

这份树状数组模板支持自动根据数据类型和数据量申请内存,单点更新和区间查询。就算不用模板,这份代码也很容易记忆。。

下面是HDU 1166,典型的树状数组题目,直接套上面的模板就可以解了。。因此不再放出完整代码,只贴上主函数,复制这两份代码再添上头文件即可运行。

bubuko.com,布布扣
int main(int agrc, char *agrv[]) {
    int T;
    scanf("%d",&T);for(int i=1;i<=T;i++){
        printf("Case %d:n",i);
        int n;
        scanf("%d",&n);
        BIT *ps=new BIT(n);
        for(int j=0;j<n;j++){
            int p;
            scanf("%d",&p);
            ps->updateDelta(j,p);
        }
        char buff[80];
        while(true){
            scanf("%s",buff);
            if(buff[0]==E) break;
            int posi,shuj;
            switch(buff[0]){
            case Q:
                int l,r;
                scanf("%d%d",&l,&r);
                printf("%dn",ps->sum(l-1,r));
                break;
            case A:
                scanf("%d%d",&posi,&shuj);
                ps->updateDelta(posi-1,shuj);
                break;
            case S:
                scanf("%d%d",&posi,&shuj);
                ps->updateDelta(posi-1,-shuj);
                break;
            }
        }
        delete ps;
    }
}
bubuko.com,布布扣

 

 

 

树状数组,布布扣,bubuko.com

树状数组

标签:style   class   blog   code   java   http   

原文地址:http://www.cnblogs.com/549294286/p/3783961.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!