题解 P3374 【【模板】树状数组 1】
恩,这是AC的第一道树状数组呢。
本蒟蒻以前遇到RMQ问题一般都用线段树或ST表,可惜ST表不支持在线修改,而线段树代码量又太大。
如今终于找到了折中方案:树状数组!!!!
代码量小,还支持修改!
树状数组也就是二叉索引树,又被称为Fenwick树,然而我个人认为它不能被严谨地成为树,因为充其量只是借用的树形结构的思想,于实现上有着较大的区别。
树状数组虽然运用范围没有线段树那么广,但是它的效率要高很多,比如线段树是$nlogn$,但树状数组是$logn$。
还有一点需要注意的是:树状数组可以区间查询,但不能运用于任意区间查询。这一点在后面会提到。
那么这个树状数组的基本思路就是
用节点ci储存和,比如:
- c1=a1
- c2=a1+a2
- c3=a3
- c4=a1+a2+a3+a4
- c5=a5
- c6=a5+a6
- c7=a7
- c8=a1+...+a8
当然这样子可能不是很容易看出内在的联系,因此不妨将其转化为二进制来观察:
- c0001=a0001
- c0010=a0001+a0010
- c0011=a0011
- c0100=a0001+a0010+a0011+a0100
- c0101=a0101
- c0110=a0101+a0110
- c0111=a0111
- c1000=a0001+...+a1000
是不是发现了什么?
没有吗?好吧。
事实上这里的规律就是cn=a(n–2^k+1)+...+an,这里的k指的是n二进制末尾0的数量。
获取2^k的操作我们称之为lowbit,其实现如下:
int lowbit(int k){ return k&(-k); }
有了lowbit操作之后,求和就很好写了:
1 int query(int x){ 2 int ans=0; 3 while(x!=0){ 4 ans+=tree[x]; 5 x-=lowbit(x); 6 } 7 return ans; 8 }
要注意一点,这里求的ans是区间[1,x]的和,想要[y,x]的和只能$query(x)-query(y-1)$。
因此我们回到了之前那个问题:树状数组不能解决所有区间查询,它只能解决如上的有关联的区间查询。
emmmm.....还有update操作:
1 void update(int x,int k){ 2 while(x<=n){ 3 tree[x]+=k; 4 x+=lowbit(x); 5 } 6 }
这个在明白了树状数组的本质之后也很好理解,就不多做叙述了。
总的来说,树状数组挺好用的,值得一学。但切记,无论如何都必须掌握线段树,因为能用树状数组解决的都能用线段树,而反之不一定如此。
另附AC代码见下:
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 using namespace std; 5 6 const int maxn=500500; 7 8 int n,m; 9 int tree[maxn<<2]; 10 11 int lowbit(int k){ 12 return k&(-k); 13 } 14 15 void update(int x,int k){ 16 while(x<=n){ 17 tree[x]+=k; 18 x+=lowbit(x); 19 } 20 } 21 22 int query(int x){ 23 int ans=0; 24 while(x!=0){ 25 ans+=tree[x]; 26 x-=lowbit(x); 27 } 28 return ans; 29 } 30 31 int main(){ 32 scanf("%d%d",&n,&m); 33 for(int i=1;i<=n;i++){ 34 int a; 35 scanf("%d",&a); 36 update(i,a); 37 } 38 for(int i=1;i<=m;i++){ 39 int a,b,c; 40 scanf("%d%d%d",&a,&b,&c); 41 if(a==1)update(b,c); 42 else printf("%d\n",query(c)-query(b-1)); 43 } 44 }