cogs577 蝗灾

题意：W*W(W <= 500000)的矩阵，初始全是0，N(N <= 200000)个询问，单点加一个值，区间查询和。

分析：

一看这题，woc这不是二维线段树吗，然而线段树是O(W^2)的空间复杂度，开不下，怎么办呢？

于是一个叫陈丹琦的人发明了一种做法，叫cdq分治。

这个分治大概可以解决一类带修改的查询问题，是一种离线分治做法。

定义solve(l, r)：对于所有k∈(l, r)，若第k个操作是查询，则计算编号为l~k-1的修改操作对k的影响。

然后分治就行了。

计算过程大概就是按x排个序，然后每一列作为一个元素，按x方向建一个一维树状数组，乱搞就行。

第一道cdq分治...

#include <cstdio>
#include <algorithm>
using namespace std;
 
const int N = 200005;
int w,q,c[500005];
struct nd {int op,x1,y1,x2,y2,z,id,ans;}a[N],b[N];
bool cmp(const nd &a, const nd &b) {return a.x1<b.x1 || (a.x1==b.x1&&a.op<b.op);}
 
int lw(int x) {return x & -x;}
void add(int x, int y) {while(x <= w) c[x] += y, x += lw(x);}
int qry(int x) {
	int r = 0;
	while(x) r += c[x], x -= lw(x);
	return r;
}
 
void gzz(int l, int r) {
	if(l == r) return;
	int m = (l+r) >> 1, cnt = 0;
	gzz(l,m),gzz(m+1,r);
	for(int i = l; i <= m; i++) if(a[i].op == 1) b[cnt++] = a[i];
	for(int i = m+1; i <= r; i++) if(a[i].op == 2) {
		b[cnt++] = a[i];
		b[cnt++] = a[i];
		b[cnt-2].x1--, b[cnt-1].x1=a[i].x2,
		b[cnt-1].op = 3;
	}
	sort(b, b+cnt, cmp);
	for(int i = 0; i < cnt; i++)
	if(b[i].op == 1) add(b[i].y1, b[i].z);
	else if(b[i].op == 2) a[b[i].id].ans -= qry(b[i].y2)-qry(b[i].y1-1);
	else a[b[i].id].ans += qry(b[i].y2)-qry(b[i].y1-1);
	for(int i = 0; i < cnt; i++)
	if(b[i].op == 1) add(b[i].y1, -b[i].z);
}
 
int main() {
	freopen("locust.in","r",stdin);
	freopen("locust.out","w",stdout);
	scanf("%d%d",&w,&q);
	for(int i = 1; i <= q; i++) {
		scanf("%d%d%d",&a[i].op,&a[i].x1,&a[i].y1);
		if(a[i].op == 1) scanf("%d",&a[i].z); else scanf("%d%d",&a[i].x2,&a[i].y2);
		a[i].id = i;
	}
	gzz(1, q);
	for(int i = 1; i <= q; i++) if(a[i].op == 2) printf("%d\n", a[i].ans);
	return 0;
}