摩基亚Mokia

P1948 - 【BOI2007】摩基亚Mokia

Description

摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。
在定位系统中，世界被认为是一个W * W的正方形区域，由1 * 1的方格组成。每个方格都有一个坐标(x,y)，1<=x,y<=W。坐标的编号从1开始。对于一个4 * 4的正方形，就有1<=x<=4,1<=y<=4（如图）：

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

Input

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

Output

对所有命令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

样例提示：
输入             输出         意义
0 4                             大小为4 * 4的全零正方形
1 2 3 3                         向(2,3)方格加入3名用户
2 1 1 3 3                     查询矩形1<=x<=3,1<=y<=3内的用户数量
                    3                查询结果
1 2 2 2                         向(2,2)方格加入2名用户
2 2 2 3 4                     查询矩形2<=x<=3,2<=y<=4内的用户数量
                    5                 查询结果
3                                 终止程序

数据规模：
1<=W<=2000000
1<=X1<=X2<=W
1<=Y1<=Y2<=W
1<=x,y<=W
0＜A＜=10000
命令1不超过160000个。
命令2不超过10000个。

Source

Balkan Olypiad in Informatics 2007,Mokia
陈丹琦分治， 分治， 树状数组

Accepted entrance

CDQ分治，把一个矩阵的询问拆成四个，运用容斥原理，给每个询问一个时间顺序。可以发现，修改只是单个格子修改，只有修改的这个各自的横坐标小于询问的横坐标，才会对这个询问产生贡献，刚好是三维的（x,y,time），于是time排序，xCDQ分治，y树状数组维护，还是一个样的，计算左区间[l,mid]对右区间[mid+1,r]的贡献，统计答案。

参考博客:http://wulala.logdown.com/posts/207262-boi-2007-mokia

#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define lowbit(p) p & (-p)
#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<queue>
#include<ctime>
#include<cmath>
#include<stack>
#include<map>
#include<set>
using namespace std;
const int N=160010*4+10000;
int tot,W,C[2000000+10],q;
struct Q{
    int x,y,k,q,v,od;//  k 种类  q Index_ans   v val or -+   id _ time
}que[N],np[N];
int Ans[N];
bool comp(const Q &a,const Q &b) {
    return a.x<b.x;
}
void add(int p,int x) {
    if(!p) return;
    while(p<=W) C[p]+=x,p+=lowbit(p);
}
int getsum(int p) {
    int sum=0;
    while(p) sum+=C[p],p-=lowbit(p);
    return sum;
}
int gi() {
    int res=0,f=1;
    char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if(ch=='-') ch=getchar(),f=-1;
    while(ch>='0'&&ch<='9') res=res*10+ch-'0',ch=getchar();
    return res*f;
}
void CDQ(int l,int r)  {
    if(l==r) return;
    int mid=(l+r)>>1;
    int l1=l,l2=mid+1;
    rep(i,l,r) if(que[i].od<=mid) np[l1++]=que[i];else np[l2++]=que[i];
    rep(i,l,r) que[i]=np[i];
    l1=l;
    rep(i,mid+1,r) if(que[i].k==2){
    for(;l1<=mid&&que[l1].x<=que[i].x;l1++) {
        if(que[l1].k==1)
        add(que[l1].y,que[l1].v);
    }
    Ans[que[i].q]+=getsum(que[i].y)*que[i].v;    
    }
    rep(i,l,l1-1) if(que[i].k==1) add(que[i].y,-que[i].v);
    CDQ(l,mid) , CDQ(mid+1,r);
}
void work() {
    while(1) {
    int Type=gi(),x,y;
    if(Type==3) break;
    if(Type==1) {
        x=gi(),y=gi();int c=gi();
        que[++tot]=(Q){x,y,1,0,c,tot};
    }
    else {
        ++q;
        x=gi(),y=gi();int u2=gi(),v2=gi();
        que[++tot]=(Q){x-1,y-1,2,q,1,tot};
        que[++tot]=(Q){x-1,v2,2,q,-1,tot};
        que[++tot]=(Q){u2,y-1,2,q,-1,tot};// y-1 v2-1
        que[++tot]=(Q){u2,v2,2,q,1,tot};
    }
    }
  
    sort(que+1,que+tot+1,comp);
    CDQ(1,tot);
    return;
}
int main() {
        freopen("mokia.in","r",stdin);
        freopen("mokia.out","w",stdout);
    
    int Type=gi();
    W=gi();
    work();
    rep(i,1,q) printf("%d
",Ans[i]);
    return 0;
}