树状数组

int lowbit(int t)
{
	return t&(-t);
}
void add(int x,int y)
{
	for(int i=x;i<=n;i+=lowbit(i))
		tree[i]+=y;
}
int getsum(int x)
{
	int ans=0;
	for(int i=x;i>0;i-=lowbit(i))
		ans+=tree[i];
	return ans;
}

 这篇笔记 会详细的讲解，使得队员们对树状数组彻底入门  而不是懵懵懂懂。

以上先给出 最常见的，三个函数。(单点更新，区间查询)
*************************************************分割线
树状数组 重点是在树状的数组
大家都知道二叉树吧
叶子结点代表A数组A[1]~A[8]

…
现在变形一下

现在定义每一列的顶端结点C[]数组
如下图

C[i]代表 子树的叶子结点的权值之和// 这里以求和举例
如图可以知道
C[1]=A[1];
C[2]=A[1]+A[2];
C[3]=A[3];
C[4]=A[1]+A[2]+A[3]+A[4];
C[5]=A[5];
C[6]=A[5]+A[6];
C[7]=A[7];
C[8]=A[1]+A[2]+A[3]+A[4]+A[5]+A[6]+A[7]+A[8];
下面观察如下图

将C[]数组的结点序号转化为二进制
1=(001) C[1]=A[1];
2=(010) C[2]=A[1]+A[2];
3=(011) C[3]=A[3];
4=(100) C[4]=A[1]+A[2]+A[3]+A[4];
5=(101) C[5]=A[5];
6=(110) C[6]=A[5]+A[6];
7=(111) C[7]=A[7];
8=(1000) C[8]=A[1]+A[2]+A[3]+A[4]+A[5]+A[6]+A[7]+A[8];
对照式子可以发现 C[i]=A[i-2^k+1]+A[i-2k+2]+…A[i]; （k为i的二进制中从最低位到高位连续零的长度）例如i=8时，k=3;
可以自行带入验证;
现在引入lowbit(x)
lowbit(x) 其实就是取出x的最低位1 换言之 lowbit(x)=2^k k的含义与上面相同 理解一下
下面说代码
int lowbit(int t)
{
return t&(-t);
}
//-t 代表t的负数 计算机中负数使用对应的正数的补码来表示
//例如 :

// t=6（0110） 此时 k=1
//-t=-6=(1001+1)=(1010)
// t&(-t)=(0010)=2=2^1
C[i]=A[i-2^k+1]+A[i-2^k+2]+......A[i];
C[i]=A[i-lowbit(i)+1]+A[i-lowbit(i)+2]+......A[i];

*************************************************分割线
区间查询
ok 下面利用C[i]数组，求A数组中前i项的和
举个例子 i=7;
sum[7]=A[1]+A[2]+A[3]+A[4]+A[5]+A[6]+A[7] ; 前i项和
C[4]=A[1]+A[2]+A[3]+A[4]; C[6]=A[5]+A[6]; C[7]=A[7];
可以推出: sum[7]=C[4]+C[6]+C[7];
序号写为二进制: sum[(111)]=C[(100)]+C[(110)]+C[(111)];

再举个例子 i=5
sum[7]=A[1]+A[2]+A[3]+A[4]+A[5] ; 前i项和
C[4]=A[1]+A[2]+A[3]+A[4]; C[5]=A[5];
可以推出: sum[5]=C[4]+C[5];
序号写为二进制: sum[(101)]=C[(100)]+C[(101)];

细细观察二进制 树状数组追其根本就是二进制的应用
结合代码

int getsum(int x)
{
	int ans=0;
	for(int i=x;i>0;i-=lowbit(i))
		ans+=C[i];
	return ans；
}
对于i=7 进行演示 
                                  7(111)          ans+=C[7]
lowbit(7)=001  7-lowbit(7)=6(110)    ans+=C[6]
lowbit(6)=010  6-lowbit(6)=4(100)    ans+=C[4]
lowbit(4)=100  4-lowbit(4)=0(000)    
对于i=5 进行演示 
                                  5(101)           ans+=C[5]
lowbit(5)=001  5-lowbit(5)=4(100)    ans+=C[4]
lowbit(4)=100  4-lowbit(4)=0(000)   

*************************************************分割线
单点更新

当我们修改A[]数组中的某一个值时 应当如何更新C[]数组呢？
回想一下 区间查询的过程，再看一下上文中列出的图

结合代码分析
void add(int x,int y)
{
for(int i=x;i<=n;i+=lowbit(i))
tree[i]+=y;
}
//可以发现 更新过程是查询过程的逆过程
//由叶子结点向上更新C[]数组

如图：
当更新A[1]时 需要向上更新C[1] ,C[2],C[4],C[8]
C[1], C[2], C[4], C[8]
写为二进制 C[(001)],C[(010)],C[(100)],C[(1000)]
1(001) C[1]+=A[1]
lowbit(1)=001 1+lowbit(1)=2(010) C[2]+=A[1]
lowbit(2)=010 2+lowbit(2)=4(100) C[4]+=A[1]
lowbit(4)=100 4+lowbit(4)=8(1000) C[8]+=A[1]

附上一道二维树状数组的例题：
//情人节的电灯泡
题目描述：
链接：https://ac.nowcoder.com/acm/problem/15172
来源：牛客网

情人节到了，小芳和小明手牵手，打算过一个完美的情人节，但是小刚偏偏也来了，当了一个明晃晃的电灯泡，小明很尴尬，就和小刚说，我交给你个任务，你完成了我俩就带你玩，否则你就回家吧。小刚很有当单身狗的觉悟，他坚决不想让小明过好情人节，同为单身狗的你能帮帮他吗？现在有一个n×n（1 <= n <= 1000）的格子，每一个格子都有一个电灯泡，可能是亮的，也可能是灭的（1代表亮， 0代表灭），现在有两种操作，一种是给你一个坐标，对于那个坐标上的灯泡，如果他是亮的，那么熄灭他，反之如果他是灭的，那么打开它。第二种操作是给你两个坐标，第一个坐标代表一个子矩阵的左上角，另一个坐标是右下角，请你求出当前子矩阵中有多少个灯泡是亮着的。燥起来吧！！！单身狗们！！！！

输入描述：

第一行两个整数，n（1 <= n <= 1000）和m（1 <= m <= 100000），分别代表正方形格子的边长和询问次数。
接下来n行，每一行有n个bool形数字（0或1），代表灯泡的状态。
接下来m行，每一行第一个数字f（1或2）代表操作的类型，如果f是1，那么接下来输入一个坐标（x, y）(1 <= x, y <= n)，对于当前位置的灯泡状态进行改变，如果是2，那么接下来输入两个坐标（x1, y1）(1 <= x1, y1 <= n)， （x2, y2）(1 <= x2, y2 <= n)，确定子矩阵的位置，输出子矩阵中亮着的灯泡数量并换行。

输出描述：
对于每一个2操作，输出子矩阵中亮着的灯泡数量并换行。

代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include<bits/stdc++.h>
using namespace std;
const int N =1010,INF=0x7f7f7f7f;
int n , m;
int c[N][N], a[N][N];
int lowbit(int t)
{
	return t&(-t);
}
void add(int x,int y,int z)
{
	for(int i=x;i<=n;i+=lowbit(i))
	{
		for(int j = y ;j <= n ; j+=lowbit(j))
		{
			c[i][j] += z;
		}
	}
}
int getsum(int x , int y)
{
	int ans=0;
	for(int i=x;i>0;i-=lowbit(i))
	{
		for(int j = y ; j > 0 ; j-=lowbit(j))
		{
			ans = ans + c[i][j];
		}
	}
	return ans;
}
int main()
{
	cin >> n >> m;
	for(int i = 1 ; i <= n ; i++)
	{
		for(int j = 1 ; j <= n ; j++)
		{
			cin >> a[i][j];
			add(i,j,a[i][j]);
		}
	}
	while(m--)
	{
		int k , x1 , x2 , y1 , y2;
		cin >> k;
		if(k==1)
		{
			cin >> x1 >> y1;
			add(x1, y1, a[x1][y1] ? -1 : 1);
			a[x1][y1] ^= 1;
		}
		else
		{
			scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
			int ans = getsum(x2, y2) - getsum(x1-1, y2) - getsum(x2, y1-1) + getsum(x1-1, y1-1);
			printf("%d\n", ans);
		}
	}
	return 0;
}