树状数组区间更新区间查询

如果给区间[l,r]同时加上x的话，每个节点的值将会如何变化呢？如果令
s(i)=加上x之前的Σ_j=1ⁱa_j
s’(i)=加上x之后的Σ_j=1ⁱa_j
那么就有
is’(i)=s(i)
l<=i<=r->s’(i)=s(i)+x ✖ (i-l+1)=s(i)+x✖i-x✖(l-1)
rs’(i)=s(i)+x✖(r-l+1)
下面记sum(bit,i(为树状数组bit的前i项和。我们构建两个树状数组bit0和bit1，并且设
Σ_j=1ⁱ=sum(bit1,i)*i+sum(bit0,i);
那么在[l,r]区间同时加上x就可以看作是

• 在bit0的l位置上加上-x(l-1)
• 在bit1的l位置上加上x
• 在bit0的r+1的位置上加上rx
• 在bit1的r+1位置上加上-x
  这四个操作，因此，查询和更新操作都可以在O(logn)的时间内完成。
  更一般地，如果操作得到的结果可以用i的n次多项式表示，那么就可以使用n+1个树状数组来进行维护了。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int N,Q;
const int maxn=100005;
long long A[maxn];
long long bit0[maxn],bit1[maxn];
long long sum(long long * b,int i){
	long long res=0;
	while(i>0){
		res+=b[i];
		i-=i&(-i);
	}
	return res;
}
void add(long long * b,int i,int x){
	while(i<=N){
		b[i]+=x;
		i+=i&(-i);
	}
}
void solve(){
	for(int i=1;i<=N;i++){
		add(bit0,i,A[i]);
	}
	int l,r,x;
	char QQ[10];
	for(int i=0;i<Q;i++){
		scanf("%s%d%d",QQ,&l,&r);
		if(QQ[0]=='C'){
			scanf("%d",&x);
			add(bit0,l,-x*(l-1));
			add(bit1,l,x);
			add(bit0,r+1,x*r);
			add(bit1,r+1,-x);
		}
		else{
			long long res=0;
			res+=sum(bit0,r)+sum(bit1,r)*r;//这里是区间求和所以要减去之前的数
			res-=sum(bit0,l-1)+sum(bit1,l-1)*(l-1);
			cout << res << endl; 
		}
	}
}
int main(){
	scanf("%d%d",&N,&Q);
	for(int i=1;i<=N;i++){
		scanf("%d",&A[i]);
	}
	solve();
	return 0;
} 

例题：POJ 3468