bzoj 1176 [Balkan2007]Mokia

1176: [Balkan2007]Mokia

Time Limit: 30 Sec   Memory Limit: 162 MB
Submit: 2355   Solved: 1060
[ Submit][ Status][ Discuss]

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范围

Source

【分析】

Orz...双倍经验题

【代码】

//bzoj 1176
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=2000005;
int S,n,m,tot,T,num;
int c[mxn],ans[mxn<<1];
struct query
{
	int opt,x,y,A,id,be;
}q[640005],nq[640005];
inline bool comp(query v,query w)
{
	if(v.x==w.x && v.y==w.y) return v.opt<w.opt;
	if(v.x==w.x) return v.y<w.y;
	return v.x<w.x;
}
inline 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;
}
inline int lowbit(int x) {return x&-x;}
inline void add(int x,int v)
{
	while(x<=n) c[x]+=v,x+=lowbit(x);
}
inline int getsum(int x) 
{
	int sum=0; 
	while(x) sum+=c[x],x-=lowbit(x);
	return sum; 
}
inline void CDQ(int l,int r)  //x轴排序，y树状数组，CDQ时间线
{
	if(l==r) return;
	int i,j,mid=(l+r)>>1,l1=l,l2=mid+1;
	fo(i,l,r)
	{
		if(q[i].id<=mid && q[i].opt==1) add(q[i].y,q[i].A);
		if(q[i].id>mid && q[i].opt==2) ans[q[i].be]+=q[i].A*getsum(q[i].y);
	}
	fo(i,l,r) if(q[i].id<=mid && q[i].opt==1) add(q[i].y,-q[i].A);
	fo(i,l,r) if(q[i].id<=mid) nq[l1++]=q[i];else nq[l2++]=q[i];
	fo(i,l,r) q[i]=nq[i];
	CDQ(l,mid),CDQ(mid+1,r);
}
int main()
{
	int i,j,opt,x1,y1,x2,y2;
	S=read(),n=read();
	while(scanf("%d",&opt) && opt!=3)
	{
		if(opt==1) q[++tot].opt=opt,q[tot].x=read(),q[tot].y=read(),q[tot].A=read(),q[tot].id=tot;
		else
		{
			x1=read(),y1=read(),x2=read(),y2=read();
			q[++tot].x=x1-1,q[tot].y=y1-1,q[tot].opt=2,q[tot].A=1,q[tot].id=tot,q[tot].be=++num;
			q[++tot].x=x2,q[tot].y=y2,q[tot].opt=2,q[tot].A=1,q[tot].id=tot,q[tot].be=num;
			q[++tot].x=x1-1,q[tot].y=y2,q[tot].opt=2,q[tot].A=-1,q[tot].id=tot,q[tot].be=num;
			q[++tot].x=x2,q[tot].y=y1-1,q[tot].opt=2,q[tot].A=-1,q[tot].id=tot,q[tot].be=num;
			ans[num]+=S*(x2-x1+1)*(y2-y1+1);
		}
	}
	sort(q+1,q+tot+1,comp);
	CDQ(1,tot);
	fo(i,1,num) printf("%d\n",ans[i]);
	return 0;
}