【BZOJ1176】[Balkan2007]Mokia

Description

维护一个W*W的矩阵，初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.

Input

第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小
接下来每行为一下三种输入之一(不包含引号):
"1 x y a"
"2 x1 y1 x2 y2"
"3"
输入1:你需要把(x,y)(第x行第y列)的格子权值增加a

输入2:你需要求出以左下角为(x1,y1),右上角为(x2,y2)的矩阵内所有格子的权值和,并输出

输入3:表示输入结束

Output

对于每个输入2,输出一行,即输入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

保证答案不会超过int范围

【BZOJ2683】简单题

题意：和上面一样

题解：依旧是分治思路

我们将询问差分，拆成四个点分别求前缀和，然后将询问和修改一起按横坐标排序；

然后在跑cdq分治的时候，我们一边递归一边按照修改的时间进行归并排序。根据cdq分治的思想，我们在执行区间[l,r]时已经处理完[l,mid]和[mid+1,r]了，那么我们只需要处理[l,mid]中的修改对[mid+1,r]中的询问造成的影响。于是我们在跑归并排序的时候，将[l,mid]中的修改按照纵坐标插入树状数组，然后将[mid+1,r]中的询问按照纵坐标在树状数组上查找，再更新一下答案就好了

这代码我已经不怎么压行了，是不是可读性高了好多。。。

#include <cstdio>

#include <cstring>

#include <iostream>

#include <algorithm>

using namespace std;

int n,num,m;

const int maxn=200010;

struct node

{

	int x,y,org,v,k;

}s[maxn];

int p[maxn],q[maxn],tr[maxn*10],ans[maxn];

//p是归并排序的顺序，q是p的复制

bool cmp(node a,node b)

{

	if(a.x==b.x&&a.org==b.org)	return a.y<b.y;

	if(a.x==b.x)	return a.org<b.org;

	return a.x<b.x;

}

void updata(int a,int val)

{

	for(int i=a;i<=n;i+=i&-i)	tr[i]+=val;

}

int query(int a)

{

	int i,ret=0;

	for(i=a;i;i-=i&-i)	ret+=tr[i];

	return ret;

}

void dfs(int l,int r)

{

	if(l==r)

	{

		p[l]=l;

		return ;

	}

	int i,mid=l+r>>1,h1=l,h2=mid+1;

	dfs(l,mid),dfs(mid+1,r);

	for(i=l;i<=r;i++)

	{

		if(h2<=r&&(h1>mid||s[p[h1]].org>=s[p[h2]].org))

		{

			if(!s[p[h2]].k)	ans[s[p[h2]].org]+=s[p[h2]].v*query(s[p[h2]].y);

			q[i]=p[h2++];

		}

		else

		{

			if(s[p[h1]].k)	updata(s[p[h1]].y,s[p[h1]].v);

			q[i]=p[h1++];

		}

	}

	for(i=l;i<=mid;i++)	if(s[p[i]].k)	updata(s[p[i]].y,-s[p[i]].v);

	for(i=l;i<=r;i++)	p[i]=q[i];

}

int main()

{

	int i,a,b,c,d,e;

	scanf("%d%d",&a,&n);

	for(i=1;e!=3;i++)

	{

		scanf("%d",&e);

		if(e==1)

		{

			scanf("%d%d%d",&a,&b,&c);

			s[++num].x=a,s[num].y=b,s[num].org=m,s[num].k=1,s[num].v=c;

		}

		if(e==2)

		{

			scanf("%d%d%d%d",&a,&b,&c,&d),a--,b--;

			s[num+1].org=s[num+2].org=s[num+3].org=s[num+4].org=++m;

			s[num+1].k=s[num+2].k=s[num+3].k=s[num+4].k=0;

			s[++num].x=a,s[num].y=b,s[num].v=1;

			s[++num].x=a,s[num].y=d,s[num].v=-1;

			s[++num].x=c,s[num].y=b,s[num].v=-1;

			s[++num].x=c,s[num].y=d,s[num].v=1;

		}

	}

	sort(s+1,s+num+1,cmp);

	dfs(1,num);

	for(i=1;i<=m;i++)	printf("%d\n",ans[i]);

	return 0;

}