BZOJ1176---[Balkan2007]Mokia （CDQ分治 + 树状数组）

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1176

CDQ第一题，warush了好久。。

CDQ分治推荐论文：

1 《从＜Cash＞谈一类分治算法的应用》 陈丹琦

2 《浅谈数据结构题的几个非经典解法》  许昊然

关于CDQ分治，两种要求：①操作不相互影响  ②可以离线处理

题目描述是有问题的，，初始时 全部为0，不是s

题意：二维平面内，两种操作，1 x y v ，位于（x,y）的值加上v.。。2 x1,y1,x2,y2，，（x1,y1） (x2,y2)分别矩形的左上角右下角，查询矩形内所有元素的和。

大致思路，在x这一维上进行分治，然后y这一维直接就可以用树状数组乱搞了。

  1 #include <cstdio>
  2 #include <cstdlib>
  3 #include <cstring>
  4 #include <algorithm>
  5 using namespace std;
  6 const int maxb = 2e6+10;
  7 struct Node
  8 {
  9     int x,y,delt;
 10     int flag,idx;
 11     Node(){}
 12     Node(int _x,int _y,int _delt,int _flag,int _idx):
 13         x(_x),y(_y),delt(_delt),flag(_flag),idx(_idx){};
 14     bool operator < (const Node &rhs)const
 15     {
 16         return x < rhs.x || (x == rhs.x && y < rhs.y);
 17     }
 18 }a[400000];
 19 struct BIT
 20 {
 21     int c[maxb],MAX;
 22     inline int lowbit(int x)
 23     {
 24         return x & -x;
 25     }
 26     void add(int x,int val)
 27     {
 28         while (x <= MAX)
 29         {
 30             c[x] += val;
 31             x += lowbit(x);
 32         }
 33     }
 34     int sum(int x)
 35     {
 36         int res = 0;
 37         while (x > 0)
 38         {
 39             res += c[x];
 40             x -= lowbit(x);
 41         }
 42         return res;
 43     }
 44 }arr;
 45
 46 //---------------------------------------------
 47 int ans[200000];
 48 void CDQ(int l,int r)
 49 {
 50     if (l  == r)
 51         return;
 52     int mid = (l + r) >> 1;
 53     CDQ(l,mid);
 54     CDQ(mid+1,r);
 55     int j = l;
 56     for (int i = mid + 1; i <= r; i++)
 57     {
 58         if (a[i].flag == 2)
 59         {
 60             for ( ; j <= mid && a[j].x <= a[i].x; j++)
 61             {
 62                 if (a[j].flag == 1)
 63                 {
 64                     arr.add(a[j].y,a[j].delt);
 65                 }
 66             }
 67             ans[a[i].idx] += arr.sum(a[i].y) * a[i].delt;
 68         }
 69     }
 70     for (int i = l; i < j; i++)
 71         if (a[i].flag == 1)
 72             arr.add (a[i].y,-a[i].delt);
 73     inplace_merge (a+l,a+mid+1,a+r+1);                   // 合并区间 [l,mid+1) [mid+1,r+1)
 74 }
 75 int main(void)
 76 {
 77     #ifndef ONLINE_JUDGE
 78         freopen("in.txt","r",stdin);
 79     #endif // ONLINE_JUDGE
 80     int s,w;
 81     while (~scanf ("%d%d",&s,&w))
 82     {
 83         int op;
 84         int tot = 1,totq = 1;
 85         arr.MAX = w;
 86         memset(ans,0,sizeof(ans));
 87         while (scanf ("%d",&op), op != 3)
 88         {
 89             if (op == 1)
 90             {
 91                 int x,y,delt;
 92                 scanf ("%d%d%d",&x,&y,&delt);
 93                 a[tot] = Node(x,y,delt,1,tot);
 94                 tot++;
 95             }
 96             if (op == 2)
 97             {
 98                 int x1,y1,x2,y2;
 99                 scanf ("%d%d%d%d",&x1,&y1,&x2,&y2);
100                 a[tot] = Node(x1-1, y1-1,  1, 2, totq); tot++;
101                 a[tot] = Node(x2,   y1-1, -1, 2, totq); tot++;
102                 a[tot] = Node(x1-1, y2,   -1, 2, totq); tot++;
103                 a[tot] = Node(x2,   y2,    1, 2, totq); tot++;
104                 totq++;
105             }
106         }
107         CDQ (1,tot-1);
108         for (int i = 1; i < totq; i++)
109             printf("%d\n",ans[i]);
110     }
111     return 0;
112 }