COJS 1752. [BOI2007]摩基亚Mokia

1752. [BOI2007]摩基亚Mokia

★★★   输入文件：mokia.in   输出文件：mokia.out   简单对比
时间限制：5 s   内存限制：128 MB

【题目描述】

摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。

在定位系统中，世界被认为是一个W*W的正方形区域，由1*1的方格组成。每个方格都有一个坐标(x,y)，1<=x,y<=W。坐标的编号从1开始。对于一个4*4的正方形，就有1<=x<=4,1<=y<=4（如图）：

请帮助Mokia公司编写一个程序来计算在某个矩形区域内有多少名用户。

【输入格式】

有三种命令，意义如下：

命令	参数	意义
0	W	初始化一个全零矩阵。本命令仅开始时出现一次。
1	x y A	向方格(x,y)中添加A个用户。A是正整数。
2	x1 y1 x2 y2	查询X1<=x<=X2，Y1<=y<=Y2所规定的矩形中的用户数量
3	无参数	结束程序。本命令仅结束时出现一次。

【输出格式】

对所有命令2，输出一个一行整数，即当前询问矩形内的用户数量。

【输入样例】

0 4

1 2 3 3

2 1 1 3 3

1 2 2 2

2 2 2 3 4

3

【输出样例】

3

5

【数据规模】

1<=W<=2000000

1<=X1<=X2<=W

1<=Y1<=Y2<=W

1<=x,y<=W

0<A<=10000

命令1不超过160000个。

命令2不超过10000个。

【来源】

Balkan Olypiad in Informatics 2007,Mokia

解题：CDQ分治

 #include <bits/stdc++.h>
 using namespace std;
 const int maxn = ;
 struct QU {
     int x1,x2,y,f,id;
     bool operator<(const QU &t) const {
         return y == t.y?id < t.id:y < t.y;
     }
 } A[maxn<<],B[maxn<<],Q[maxn<<];
 int C[maxn<<],Li[maxn<<],ans[maxn<<],tot,cnt,ask;
 void add(int i,int val) {
     for(; i < maxn; i += i&(-i)) C[i] += val;
 }
 int sum(int i) {
     int ret = ;
     for(; i > ; i -= i&(-i)) ret += C[i];
     return ret;
 }
 void cdq(int L,int R) {
     if(R <= L) return;
     int mid = (L + R)>>;
     cdq(L,mid);
     cdq(mid + ,R);
     int a = ,b = ,j = ;
     for(int i = L; i <= mid; ++i)
         if(Q[i].id == -) A[a++] = Q[i];
     for(int i = mid+; i <= R; ++i)
         if(Q[i].id != -) B[b++] = Q[i];
     sort(A,A + a);
     sort(B,B + b);
     for(int i = ; i < b; ++i) {
         for(; j < a && A[j].y <= B[i].y; ++j) add(A[j].x1,A[j].f);
         ans[B[i].id] += B[i].f*sum(B[i].x2);
         ans[B[i].id] -= B[i].f*sum(B[i].x1);
     }
     for(int i = ; i < j; ++i)
         add(A[i].x1,-A[i].f);
 }
 int main() {
     int n,op;
     freopen ( "mokia.in", "r", stdin ) ;
     freopen ( "mokia.out", "w", stdout ) ;
     while(~scanf("%*d%d",&n)) {
         memset(C,,sizeof C);
         memset(ans,,sizeof ans);
         tot = cnt = ask = ;
         while(scanf("%d",&op),op < ) {
             if(op == ) {
                 scanf("%d%d%d",&Q[cnt].x1,&Q[cnt].y,&Q[cnt].f);
                 Li[tot++] = Q[cnt].x1;
                 Q[cnt++].id = -;
             } else if(op == ) {
                 scanf("%d%d%d%d",&Q[cnt].x1,&Q[cnt+].y,&Q[cnt+].x2,&Q[cnt].y);
                 Li[tot++] = Q[cnt].x1-;
                 Li[tot++] = Q[cnt+].x2;
                 --Q[cnt+].y;
                 Q[cnt].x2 = Q[cnt+].x2;
                 Q[cnt+].x1 = Q[cnt].x1;
                 Q[cnt].id = Q[cnt+].id = ask++;
                 Q[cnt++].f = ;
                 Q[cnt++].f = -;
             }
         }
         sort(Li,Li + tot);
         tot = unique(Li,Li + tot) - Li;
         for(int i = ; i < cnt; ++i) {
             if(Q[i].id == -) Q[i].x1 = lower_bound(Li,Li + tot,Q[i].x1) - Li + ;
             else {
                 Q[i].x1 = lower_bound(Li,Li + tot,Q[i].x1 - ) - Li + ;
                 Q[i].x2 = lower_bound(Li,Li + tot,Q[i].x2) - Li + ;
             }
         }
         cdq(,cnt-);
         for(int i = ; i < ask; ++i)
             printf("%d\n",ans[i]);
     }
     return ;
 }