树状数组（区间修改，区间查询模板题）

区间修改，区间查询
这是一道模板题。

给定数列 a[1],a[2],…,a[n]，你需要依次进行 q 个操作，操作有两类：

1 l r x：给定 l,r,x，对于所有 i∈[l,r]，将 a[i] 加上 x（换言之，将 a[l],a[l+1],…,a[r] 分别加上 x）；
2 l r：给定 l,r，求 ∑ri=la[i] 的值（换言之，求 a[l]+a[l+1]+⋯+a[r] 的值）。

Input

第一行包含 2 个正整数 n,q，表示数列长度和询问个数。保证 1≤n,q≤1e6。
第二行 n 个整数 a[1],a[2],…,a[n]，表示初始数列。保证 |a[i]|≤1e6。
接下来 q 行，每行一个操作，为以下两种之一：

1 l r x：对于所有 i∈[l,r]，将 a[i] 加上 x；
2 l r：输出 ${\sum }_{i=l}^{r}a[i]$的值。
保证 1≤l≤r≤n, |x|≤1e6。

Output
对于每个 2 l r 操作，输出一行，每行有一个整数，表示所求的结果。

Example

样例输入
5 10
2 6 6 1 1
2 1 4
1 2 5 10
2 1 3
2 2 3
1 2 2 8
1 2 3 7
1 4 4 10
2 1 2
1 4 5 6
2 3 4
样例输出
15
34
32
33
50

Hint

对于所有数据，1≤n,q≤1e6, |a[i]|≤1e6, 1≤l≤r≤n, |x|≤1e6。

代码

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=1e6+10;
ll b[N],c[N];
int q[N];
int n,m;
int lowbit(int x)
{
	return x&(-x);
}
void add(int u,int x)
{
	for(int i=u;i<=n;i+=lowbit(i))
	{
		b[i]+=x;
		c[i]+=(ll)u*x;
	}
}
ll query(int x)
{
	ll res=0;
	for(int i=x;i;i-=lowbit(i))
		res+=(x+1)*b[i]-c[i];
	return res;
 } 
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&q[i]);
		add(i,q[i]-q[i-1]);
	}
	while(m--)
	{
		int f,l,r,x;
		scanf("%d",&f);
		if(f==1)
		{
			scanf("%d%d%d",&l,&r,&x);
			add(l,x);
			add(r+1,-x);
		}
		else
		{
			scanf("%d%d",&l,&r);
			printf("%lld\n",query(r)-query(l-1));
		}
	}
	return 0;
 } 

推导证明有时间补上。