POJ 3468

题意：

树状数组区间更新区间求和

总结树状数组：

1、单点更新  区间求和（单点求值）

最简单最基础的操作  数组中存的东西就是要求和的数字

2、区间更新  单点求值

维护查分序列  例如要求a[i]  数组中存的东西应该是d[i] （d[i] = a[i] - a[i-1]）

那么 a[i] = d[1] + d[2] + d[3] + … + d[i]

如果要区间[l,r]每个元素都+v  那么只需让d[l] + v  再让d[r+1] - v

3、区间更新  区间求和

求 a[1] + a[2] + a[3] + … + a[x]  沿用d[i]思想

sum = sigma(a[i]) = sigma( sigma(d[i]) ) = sigma( d[i]*(x-i+1) ) = sigma( d[i]*(x+1) ) - sigma( d[i]*i )

因此需要维护d（即d[i]）和di（即d[i]*i）两个树状数组

求值时  ans(l,r) = sum(d,r)*(r+1) - sum(d,l-1)*l - ( sum(di,r) - sum(di,l-1) )

更新[l,r]时分开维护d树状数组和di树状数组

add(d,l,v); add(d,r+1,-v);  因为 d[l] = a[l](新的) - a[l-1] = a[l](原来的) + v - a[l-1] = d[l](原来的) + v 同理r+1处 - v

add(di,l,v*l); add(di,r+1,-v*(r+1)); 在d[i]的变化基础上乘i  一样是l处加  r+1处减

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
#define M 100010

__int64 a[M],d[M],di[M];
int n,m;

int lowbit(int x)
{
    return x&(-x);
}

void add(__int64 a[],int x,__int64 v)
{
    while(x<=n)
    {
        a[x]+=v;
        x+=lowbit(x);
    }
}

__int64 sum(__int64 a[],int x)
{
    __int64 res=0;
    while(x)
    {
        res+=a[x];
        x-=lowbit(x);
    }
    return res;
}

int main()
{
    int f,t;
    char op[5];
    __int64 val;
    scanf("%d%d",&n,&m);
    for(f=1;f<=n;f++)
    {
        scanf("%I64d",&a[f]);
        a[f]+=a[f-1];
    }
    while(m--)
    {
        scanf("%s",op);
        if(op[0]=='Q')
        {
            scanf("%d%d",&f,&t);
            printf("%I64d\n",a[t]-a[f-1]+sum(d,t)*(t+1)-sum(d,f-1)*f-sum(di,t)+sum(di,f-1));
        }
        else if(op[0]=='C')
        {
            scanf("%d%d%I64d",&f,&t,&val);
            add(d,f,val);
            add(d,t+1,-val);
            add(di,f,val*f);
            add(di,t+1,val*(-t-1));
        }
    }
	return 0;
}