poj 3468 树状数组 区间修改 区间求和

这是利用树状数组进行 区间求和区间修改，写起来比线段树轻松许多，具体的原理参考大佬博客：https://blog.csdn.net/fsahfgsadhsakndas/article/details/52650026

代码如下

/*
POJ 3468 树状数组  区间修改 和 区间查询
原理：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/52650026 
c1 和 c2 是两个差分数组，所以输入数据时  c1[i]=a[i]-a[i-1]，
 
*/
#include <stdio.h>
#include <iostream>
#include <string.h>
using namespace std;
#define MAX 100020

__int64 c1[MAX],c2[MAX];//两个差分树状数组，这题数据太大所以__int64 
__int64 n,q;
__int64 lowbit(int x){ return x&(-x);}

__int64 sum(__int64* r,__int64 x) //求和 
{
   __int64 ans = 0;
   for (__int64 i = x;i>0;i-=lowbit(i))
    ans += r[i];
   return ans;
}
void add(__int64 *r,__int64 x,__int64 val)//单点更新 
{
  for (__int64 i=x;i<=n;i+=lowbit(i))
     r[i]+=val;
}

int main ()
{ 
  __int64 a,b,c,s1,s2;
  memset(c1,0,sizeof(c1));
  memset(c2,0,sizeof(c2));
  b = 0;
  scanf ("%I64d%I64d",&n,&q);
  for (int i=1;i<=n;++i)
   {
        scanf ("%I64d",&a);
     add(c1,i,a-b);//差分数组c1[n]，c1[i]=a[i]-a[i-1]
     add(c2,i,(i-1)*(a-b));//C2[I] = C1[I]*(I-1);
     b = a;
   }
    char ch = 0;
   for (int i=0;i<q;++i)
    {
      cin >> ch;
      if (ch =='C')
       {
            scanf ("%I64d%I64d%I64d",&a,&b,&c) ;//区间更新 
            add(c1,a,c);add(c1,b+1,-c) ; //一正一负 
            add(c2,a,(a-1)*c);add(c2,b+1,(b)*(-c));
       }
      else
      {
          scanf ("%I64d%I64d",&a,&b) ;
          s1 = (a-1)*sum(c1,a-1) - sum(c2,a-1);//注意这里是a - 1,因为区间a --b 包含a  
          s2 =  (b)*sum(c1,b) - sum(c2,b);
          printf ("%I64d
",s2-s1);
      }
    }
 return 0;
}