POJ 2352（顺路讲解一下树状数组）

接触到的第一道树状数组的题，AC之后感觉对树状数组思想的理解明显清晰了很多，入门必备呀。

先来讲讲树状数组吧。

上个图

如果要求区间和，例如求a3->a8，我们先用数组c来记录下来各个区间的和，那么就可以直接c8-c2就可以得到a3->a8的区间和了。

c[n]数组的下标是用来记录数组a[1]->a[n]的区间和。

再想，我们如果要求s->t的和可以怎么求？

是不是可以用(1->t)减去(1->s-1)？

同理，树状数组也可以用这种方式求和。比如求a[2]的话就可以用c[2]-c[1]来得到a[2]。

关于c数组的记录方法，就不得不提到x&(-x)这个表达式。

x&(-x)是用来计算每次查询或者更新时候c数组的偏移量。

通常写成函数的形式。

int lowbit(int x){    return x & (-x);}

注意：计算机中负数是用的补码表示的，正数是用的原码表示，可以自己拿张草稿纸来实现一下x&(-x)的计算

关于树状数组的求区间和，如果求区间[1,8]，那么直接c[8]就好，但是如果要求区间[4,8]该怎么办？

那么可以用区间[1,8]的和减去区间[1,3]的和就得到了区间[4,8]的和。

那又如何求区间[1,3]呢？我们可以c[2]+c[3]就得到了区间[1,3]的和了。

下面是求区间[1,x]和的函数：

int sum(int x){    int s = 0;    while (x > 0)    {        s += c[x];        x -= lowbit(x);    }    return s;}

以求区间[1,3]为例，首先x等于3代入方程，3 > 0，则s = c[3]，然后执行x-=lowbit(x);

让我们进入lowbit函数，首先3的二进制原码为0000 0011，-3为3的补码，则-3的二进制码为1111 1101，进行&运算之后为1，所以x-=1，此时x为2。由于2>0，所以s此时等于c[3] + c[2]。执行x-=lowbit(x)之后x = 0。循环结束，此时s已经是区间[1,3]的和了。

我们再来看树状数组的更新函数：

int update(int x， int num){    while (x <= MAX)    {        c[x] += num;        x += lowbit(x);    }}

和线段树一样，树状数组也是需要对节点所影响到的所有节点进行更新，采取从根到顶的方式。也就是对每个影响到的节点都加上更改信息。

树状数组时间复杂度O(log n)

————————————————————————分割线——————————————————————————————

题意：有n个星星节点，存在星星节点左下角（包括正左和正下）的其他星星节点，则该星星节点比它左下角的星星节点大，level 0表示该星星节点没有比他还小的节点，level 1表示存在一个比该星星节点小的点。输出统计好的每个level等级存在多少星星节点。

解题思路：为什么要用树状数组呢，因为如果你用for循环统计的话，由于数据很大然后又有很多组数据需要统计，那么肯定是会超时的，所以此时需要一种高效的数据结构（感觉像一句废话TAT），树状数组类似于线段树，能够很高效的解决区间问题，将这道题提炼一下其实也就是一个统计区间和的问题，线段树写起来好麻烦的=。=于是乎用了树状数组。

直接上代码：

/*因为是按照y升序输入，所以后面输入对前面输入并无影响**前x与后x如果相同，那么肯定后x是包含前x的，因为是按照y升序**即前后x虽是在同一列，但后x肯定在前x上面，即包含前x**PS:树状数组下标从1开始*/#include <iostream>#include <cstdio>#include <cstring>using namespace std;int level[32001];int c[32001];int lowbit(int x){    return x & (-x);}int sum(int x){    int s = 0;    while (x > 0)    {        s += c[x];        x -= lowbit(x);    }    return s;}int update(int x){    while (x <= 32001)    {        c[x]++;        x += lowbit(x);    }}int main(){    int n;    int x, y;    while(~scanf("%d", &n))    {        int N = n;        memset(level, 0, sizeof(level));        memset(c, 0, sizeof(c));        while (n--)        {            scanf("%d %d", &x, &y);            level[sum(x+1)]++;            update(x+1);        }        for (int i = 0; i <= N - 1; i++)            printf("%d\n", level[i]);    }    return 0;}