2683: 简单题
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 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
1 2 3 3
2 1 1 3 3
1 2 2 2
2 2 2 3 4
3
Sample Output
3
5
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;
}