树状数组区间修改加区间查询

其实之前在K大数查询中就已经用到了,只是一直没有说明
所以今天就来补个欠账。
感觉单点修改、区间查询和区间修改、单点查询没什么必要讲,这里就只讲区间修改、区间查询(其实也不难)。
设原数组第\(i\)位的值为\(a_i\)\(d_i=a_i-a_{i-1}\),则有(这里认为\(a_0=0\)):

\[a_x=\sum_{i=1}^x d_i \]

所以有:

\[\sum_{i=1}^x a_i= \sum_{i=1}^x \sum_{j=1}^i d_j =\sum_{i=1}^x(x-i+1)d_i \]

于是我们得到了:

\[\sum_{i=1}^x a_i=(x+1)\sum_{i=1}^x d_i-\sum_{i=1}^x d_i \times i \]

于是我们把原数组差分后维护两个树状数组,一个维护\(d_i\),一个维护\(d_i \times i\)
这样区间求和时可以在两个树状数组中查询得到前缀和,区间修改时就是差分数组的修改,每次修改两个点即可。
具体代码如下:

void add(int x,int y){for(int i=x;i<=n;i+=i&(-i)) c1[i]+=y,c2[i]+=(long long)x*y;}//给差分数组中的位置x加上y
long long sum(int x){//查询前x项的和
    long long ans(0);
    for(int i=x;i;i-=i&(-i)) ans+=(x+1)*c1[i]-c2[i];
    return ans;
}

其中\(c{1_i}\)维护的是\(d_i\)\(c{2_i}\)维护的是\(d_i\times i\)
比线段树好写多了是不?

posted @ 2016-09-12 20:36  lcf2000  阅读(8980)  评论(8编辑  收藏  举报