BZOJ 1176: [Balkan2007]Mokia

1176: [Balkan2007]Mokia

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 2012  Solved: 896
[Submit][Status][Discuss]

Description

维护一个W*W的矩阵，初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.

Input

第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小

接下来每行为一下三种输入之一(不包含引号):

"1 x y a"

"2 x1 y1 x2 y2"

"3"

输入1:你需要把(x,y)(第x行第y列)的格子权值增加a

输入2:你需要求出以左下角为(x1,y1),右上角为(x2,y2)的矩阵内所有格子的权值和,并输出

输入3:表示输入结束

Output

对于每个输入2,输出一行,即输入2的答案

Sample Input

0 4
1 2 3 3
2 1 1 3 3
1 2 2 2
2 2 2 3 4
3

Sample Output

3
5

HINT

保证答案不会超过int范围

Source

 

[Submit][Status][Discuss]

CDQ分治。

 #include <bits/stdc++.h>
 const int maxn = ;
 const int maxw = ;
 int tr[maxw], w;
 inline void add(int t, int k) {
     for (; t <= w; t += t&-t)
         tr[t] += k;
 }
 inline int ask(int t) {
     int ret = ;
     for (; t; t -= t&-t)
         ret += tr[t];
     return ret;
 }
 inline int ask(int a, int b) {
     return ask(b) - ask(a - );
 }
 struct data {
    int k, a, b, c, d, t, p;
 }s[maxn]; int tot = ;
 int ans[maxn], cnt = ;
 inline void add(void) {
     scanf("%d%d%d", &s[tot].a, &s[tot].b, &s[tot].c), s[tot].t = tot, s[tot++].k = ;
 }
 inline void qry(void) {
     int x1, x2, y1, y2; scanf("%d%d%d%d", &x1, &y1, &x2, &y2); --x1;
     s[tot].k = , s[tot].a = x1, s[tot].b = y1, s[tot].c = y2, s[tot].d = -, s[tot].t = tot, s[tot++].p = cnt;
     s[tot].k = , s[tot].a = x2, s[tot].b = y1, s[tot].c = y2, s[tot].d = +, s[tot].t = tot, s[tot++].p = cnt++;
 }
 inline bool cmp(const data &a, const data &b) {
     if (a.a != b.a)return a.a < b.a;
     else return a.k < b.k;
 }
 void solve(int l, int r) {
     if (l >= r)return;
     int mid = (l + r) >> ;
     solve(l, mid); solve(mid + , r);
     std::sort(s + l, s + r + , cmp);
     for (int i = l; i <= r; ++i)
         if (s[i].k ==  && s[i].t <= mid)
             add(s[i].b, s[i].c);
         else if (s[i].k ==  && s[i].t > mid)
             ans[s[i].p] += ask(s[i].b, s[i].c) * s[i].d;
     for (int i = l; i <= r; ++i)
         if (s[i].k ==  && s[i].t <= mid)
             add(s[i].b, -s[i].c);
 }
 signed main(void) {
     scanf("%*d%d", &w);
     int k; while (scanf("%d", &k), k != )
         (k == ) ? add() : qry();
     solve(, tot - );
     for (int i = ; i < cnt; ++i)
         printf("%d\n", ans[i]);
 }

@Author: YouSiki