bzoj 2683 简单题

2683: 简单题

Time Limit: 50 Sec  Memory Limit: 128 MB
Submit: 1375  Solved: 554
[ Submit][ Status][ Discuss]

Description

你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作：

 

命令	参数限制	内容
1 x y A	1<=x,y<=N，A是正整数	将格子x,y里的数字加上A
2 x1 y1 x2 y2	1<=x1<= x2<=N 1<=y1<= y2<=N	输出x1 y1 x2 y2这个矩形内的数字和
3	无	终止程序

Input

输入文件第一行一个正整数N。

接下来每行一个操作。

 

Output

对于每个2操作，输出一个对应的答案。

 

Sample Input

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

1<=N<=500000,操作数不超过200000个，内存限制20M。

对于100%的数据，操作1中的A不超过2000。

【分析】

x排序，y树状数组，时间线CDQ分治

【代码】

//bzoj 2683 简单题 
#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=500005;
int n,m,tot,T,num;
int c[mxn],tim[mxn],ans[mxn<<1];
struct query
{
	int opt,x,y,A,id,be;
}q[mxn<<1],nq[mxn<<1];
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;
	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;
		}
	}
	sort(q+1,q+tot+1,comp);
	CDQ(1,tot);
	fo(i,1,num) printf("%d\n",ans[i]);
	return 0;
}