COGS 577 蝗灾

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;
}