DESCRIPTION
C国国土辽阔，地大物博......但是最近却在闹蝗灾.....
我们可以把C国国土当成一个W×W的矩阵，你会收到一些诸如（X，Y，Z）的信息，代表（X，Y）这个点增多了
Z只蝗虫，而由于C国政府机关比较臃肿，为了批复消灭蝗虫的请求需要询问一大堆的问题......每个询问形如
（X1,Y1,X2,Y2），询问在（X1,Y1,X2,Y2）范围内有多少蝗虫（请注意询问不会改变区域内的蝗虫数），
你作为一个C国的程序员，需要编一个程序快速的回答所有的询问。
NOTICE
C国一开始没有蝗虫。
INPUT
输入文件的第一行包括一个整数W，代表C国国土的大小。第二行有一个整数N，表示事件数。接下来有N行表示N个事件，以（1 X Y Z）的形式或（2,X1,Y1,X2,Y2）的形式给出，分别代表蝗虫的增加和询问。
OUTPUT
对于每个询问输出一个整数表示需要的结果。
SAMPLE INPUT
locust.in
5
8
2 4 1 4 2
1 3 1 8
1 4 4 4
2 1 3 4 4
1 1 5 1
1 4 4 5
2 2 2 5 4
2 3 2 4 4
SAMPLE OUTPUT
locust.out
0
4
9
9
数据范围：
10%的数据满足W<=100,N<=100;
30%的数据满足W<=2000,N<=5000;
50%的数据满足W<=100000,N<=50000;
100%的数据满足W<=500000,N<=200000,每次蝗虫增加数不超过1000;
时间限制：
2s
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
CDQ分治+树状数组~
x坐标排序，y坐标树状数组，时间CDQ分治。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;int n,m,k,x,y,x1,y1,x2,y2,val,c[500001],cnt,tot,ans[200001];struct node{bool flag;int x,y,val,pos,id;
}a[800001],b[800001];int read()
{int x=0,f=1;char ch=getchar();while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;
}bool operator < (node u,node v)
{return u.x==v.x ? u.pos<v.pos:u.x<v.x;
}void add(int u,int v)
{for(;u<=n;u+=u&(-u)) c[u]+=v;
}int cal(int u)
{int now=0;for(;u;u-=u&(-u)) now+=c[u];return now;
}void sol(int l,int r)
{if(l>=r) return;int mid=l+r>>1,l1=l,l2=mid+1;for(int i=l;i<=r;i++)if(!a[i].flag && a[i].pos<=mid) add(a[i].y,a[i].val);else if(a[i].flag && a[i].pos>mid) ans[a[i].id]+=cal(a[i].y)*a[i].val;for(int i=l;i<=r;i++) if(!a[i].flag && a[i].pos<=mid) add(a[i].y,-a[i].val);for(int i=l;i<=r;i++)if(a[i].pos<=mid) b[l1++]=a[i];else b[l2++]=a[i];for(int i=l;i<=r;i++) a[i]=b[i];sol(l,mid);sol(mid+1,r);
}int main()
{freopen("locust.in","r",stdin);freopen("locust.out","w",stdout);n=read();m=read();for(int i=1;i<=m;i++){k=read();if(k==1) x=read(),y=read(),val=read(),a[++tot]=(node){0,x,y,val,tot,0};else{x1=read()-1;y1=read()-1;x2=read();y2=read();cnt++;a[++tot]=(node){1,x1,y1,1,tot,cnt};a[++tot]=(node){1,x2,y1,-1,tot,cnt};a[++tot]=(node){1,x1,y2,-1,tot,cnt};a[++tot]=(node){1,x2,y2,1,tot,cnt};}}sort(a+1,a+tot+1);sol(1,tot);for(int i=1;i<=cnt;i++) printf("%d\n",ans[i]);return 0;
}