树状数组和线段树都可以理解为前缀和,但是其的优点在于能够改变其中的某一个值在进行求前缀和。
树状数组
注意!:写树状数组的时候集合从1开始
w[],tr[]: tr[1]=w[1];
tr[2]=tr[1]+w[2]=w[1]+w[2];
tr[3]=w[3];
tr[4]=w[4]+tr[3]+tr[2]=w[1]+···+w[4];······
函数/方法
1、add(a,b):用来给集合的某一个数加上一个值;
当对于a位置上的数加上一个b,对tr[i]有影响的只有i+=lowbit上的数。
例如只改变a[5]的话,那么受改变影响的只有tr[5]、tr[6]、tr[8],他们就是从5这个位置开始,到6=5+1(5是奇数,lowbit返回的是1); 8=6+2(6的二进制数后面有一个0,返回2^1=2)。
void add(int a,int b){
for(int i=a;i<N;i+=lowbit(i)) tr[i]+=b;
}
2、lowbit(k):用来得到这个数的二进制数后面存在多少个0;
推导根据:https://blog.csdn.net/oyoung_2012/article/details/79932394
奇数返回1,偶数返回2^k,k表示二进制表示的数中后面0的个数;
int lowbit(int x){
return x&-x;
}
3、query(a):用来得到从0到a的前缀和;
例如求0~7的前缀和就是res=tr[7]+tr[6]+tr[4]+tr[0];
int query(int a){
int res=0;
for(int i=a;i;i-=lowbit(i)) res+=tr[i];
return res;
}
例子
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=1e5+10;
int n,m;
int w[N],tr[N];
int lowbit(int x){
return x&-x;
}
void add(int a,int b){
for(int i=a;i<N;i+=lowbit(i)) tr[i]+=b;
}
int query(int a){
int res=0;
for(int i=a;i;i-=lowbit(i)) res+=tr[i];
return res;
}
int main(){
scanf("%d %d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",&w[i]);
for(int i=1;i<=n;i++) add(i,w[i]);
while(m--){
int k,a,b;
scanf("%d %d %d",&k,&a,&b);
if(k==1) add(a,b);
else{
printf("%d",query(b)-query(a-1));
}
}
return 0;
}
线段树
注意!:线段树的时候集合从1开始
函数/方法
1、push_up(u):用来将两个子节点的和合并到父节点。
void push_up(int u){
tr[u].sum=tr[u*2].sum+tr[u*2+1].sum;
}
2、build(u,l,r):用来将w[]集合的数合并到tr[]中。
void build(int u,int l,int r){
if(l==r) tr[u]={l,r,w[r]};
else{
tr[u]={l,r};
int mid=l+r>>1;
build(u*2,l,mid);
build(u*2+1,mid+1,r);
push_up(u);
}
}
3、modify(u,a, b):从子节点递归修改。
void modify(int u,int a,int b){
if(tr[u].l==tr[u].r) tr[u].sum+=b;
else{
int mid=tr[u].l+tr[u].r>>1;
if(a<=mid) modify(2*u,a,b);
else modify(2*u+1,a,b);
push_up(u);
}
}
4、query(u,a,b):从在区间[a,b]的各个最大子节点递归求和。
int query(int u,int a,int b){
if(tr[u].l>=a&&tr[u].r<=b) return tr[u].sum;
int mid=tr[u].l+tr[u].r>>1;
int sum=0;
if(mid>=a) sum+=query(2*u,a,b);
if(mid+1<=b) sum+=query(2*u+1,a,b);
return sum;
}
例子
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=1e5+10;
int n,m;
int w[N];
struct Node{
int l,r;
int sum;
}tr[4*N];
void push_up(int u){
tr[u].sum=tr[u*2].sum+tr[u*2+1].sum;
}
void build(int u,int l,int r){
if(l==r) tr[u]={l,r,w[r]}; //这里记得给tr结构体的各个参数赋值
else{
tr[u]={l,r}; //这里同上,这里是经常出现debug的地方
int mid=l+r>>1;
build(u*2,l,mid);
build(u*2+1,mid+1,r);
push_up(u);
}
}
void modify(int u,int a,int b){
if(tr[u].l==tr[u].r) tr[u].sum+=b;
else{
int mid=tr[u].l+tr[u].r>>1;
if(a<=mid) modify(2*u,a,b);
else modify(2*u+1,a,b);
push_up(u);
}
}
int query(int u,int a,int b){
if(tr[u].l>=a&&tr[u].r<=b) return tr[u].sum;
int mid=tr[u].l+tr[u].r>>1;
int sum=0;
if(mid>=a) sum+=query(2*u,a,b);
if(mid+1<=b) sum+=query(2*u+1,a,b);
return sum;
}
int main(){
cin>>n>>m;
for(int i=1;i<=n;i++) scanf("%d",&w[i]);
build(1,1,n);
int k,a,b;
while(m--){
scanf("%d%d%d",&k,&a,&b);
if(k==1) modify(1,a,b);
else printf("%d\n",query(1,a,b));
}
return 0;
}