[BZOJ1176][Balkan2007]Mokia

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
char buf[10000000], *ptr = buf - 1;
inline int readint() {
    int f = 1, n = 0;
    char ch = *++ptr;
    while (ch < ‘0‘ || ch > ‘9‘) {
        if (ch == ‘-‘) f = -1;
        ch = *++ptr;
    }
    while (ch <= ‘9‘ && ch >= ‘0‘) {
        n = (n << 1) + (n << 3) + ch - ‘0‘;
        ch = *++ptr;
    }
    return f * n;
}
int s, w;
int C[2000000 + 10], T[2000000 + 10] = { 0 }, flag = 0;
inline int Query(int pos) {
    int ret = 0;
    for (int i = pos; i; i -= i & -i)
        if (T[i] == flag) ret += C[i];
    return ret;
}
inline void Update(int pos, int val) {
    for(int i = pos; i <= w; i += i & -i)
        if (T[i] != flag) {
            T[i] = flag;
            C[i] = val;
        }
        else C[i] += val;
}
struct Que {
    int type, x, y, val;
    Que() {}
    Que(int _type, int _x, int _y, int _val) : type(_type), x(_x), y(_y), val(_val) {}
}q[200000 + 10], tp[200000 + 10];
int ans[10000 + 10] = {0};
void CDQ(int l, int r) {
    if (l == r) return;
    int mid = l + r >> 1, ll = l, rr = mid + 1, tcnt = 0;
    CDQ(l, mid); CDQ(mid + 1, r);
    flag++;
    while (ll <= mid && rr <= r) {
        if (q[ll].x <= q[rr].x) {
            if (q[ll].type == 1) Update(q[ll].y, q[ll].val);
            tp[++tcnt] = q[ll++];
        }
        else {
            if (q[rr].type == 2) ans[q[rr].val] += Query(q[rr].y);
            else if(q[rr].type == 3) ans[q[rr].val] -= Query(q[rr].y);
            tp[++tcnt] = q[rr++];
        }
    }
    while (ll <= mid) tp[++tcnt] = q[ll++];
    while (rr <= r) {
        if (q[rr].type == 2) ans[q[rr].val] += Query(q[rr].y);
        else if(q[rr].type == 3) ans[q[rr].val] -= Query(q[rr].y);
        tp[++tcnt] = q[rr++];
    }
    for (int i = 1; i <= tcnt; i++) q[l + i - 1] = tp[i];
}
int sz = 0;
int main() {
    fread(buf, sizeof(char), sizeof(buf), stdin);
    s = readint();
    w = readint();
    int opt, x1, y1, x2, y2, a, qtime = 0;
    while ((opt = readint()) != 3) {
        x1 = readint();
        y1 = readint();
        if (opt == 1) {
            a = readint();
            q[++sz] = Que(1, x1, y1, a);
        }
        else {
            qtime++;
            x2 = readint();
            y2 = readint();
            q[++sz] = Que(2, x2, y2, qtime);
            q[++sz] = Que(2, x1 - 1, y1 - 1, qtime);
            q[++sz] = Que(3, x1 - 1, y2, qtime);
            q[++sz] = Que(3, x2, y1 - 1, qtime);
        //    ans[qtime] = (x2 - x1 + 1) * (y2 - y1 + 1) * s;
        }
    }
    CDQ(1, sz);
    for (int i = 1; i <= qtime; i++)
        printf("%d\n", ans[i]);
    return 0;
}