2683: 简单题
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1779 Solved: 720
[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。
/* CDQ分治。 用CDQ分治做了线段树练习之后感觉大体脉络了解了一点,再做这个题,发现基本思路是一样的,就是会出现其他算法的插入,比如说加一点数据结构维护多维什么的。 这个题用了树状数组,注意树状数组的插入和询问是在哪个位置使用的 */ #include<iostream> #include<cstdio> #include<algorithm> #define maxn 800010 using namespace std; int n,tot,T,ans[maxn],tr[maxn]; struct node{ int t;//第几次操作 int x,y;//操作的位置 int z;//增加量 int belong;//属于哪次查询 int op;//操作类型 bool operator < (const node & q1)const{ if(x!=q1.x)return x<q1.x; if(y!=q1.y)return y<q1.y; return op<q1.op; } }q[maxn],tmp[maxn]; void modify(int x,int v){ while(x<=n){ tr[x]+=v; x+=x&(-x); } } int query(int x){ int sum=0; while(x){ sum+=tr[x]; x-=x&(-x); } return sum; } void work(int l,int r){ if(l==r)return; int mid=(l+r)>>1; for(int i=l;i<=r;i++){ if(q[i].t<=mid&&q[i].op==1)modify(q[i].y,q[i].z); if(q[i].t>mid&&q[i].op==2) ans[q[i].belong]+=query(q[i].y)*q[i].z; } for(int i=l;i<=r;i++) if(q[i].t<=mid&&q[i].op==1)modify(q[i].y,-q[i].z);//进行恢复操作 int l1=l,l2=mid+1; for(int i=l;i<=r;i++){ if(q[i].t<=mid)tmp[l1++]=q[i]; else tmp[l2++]=q[i]; } for(int i=l;i<=r;i++)q[i]=tmp[i]; work(l,mid);work(mid+1,r); } int main(){ freopen("Cola.txt","r",stdin); scanf("%d",&n); int opt,x1,y1,x2,y2,x,y,z; while(1){ scanf("%d",&opt); if(opt==3)break; if(opt==1){ scanf("%d%d%d",&x,&y,&z); q[++tot].x=x;q[tot].y=y;q[tot].t=tot;q[tot].op=1;q[tot].z=z; } if(opt==2){ scanf("%d%d%d%d",&x1,&y1,&x2,&y2); q[++tot].x=x1-1;q[tot].y=y1-1;q[tot].t=tot;q[tot].op=2;q[tot].z=1;q[tot].belong=++T; q[++tot].x=x1-1;q[tot].y=y2;q[tot].t=tot;q[tot].op=2;q[tot].z=-1;q[tot].belong=T; q[++tot].x=x2;q[tot].y=y1-1;q[tot].t=tot;q[tot].op=2;q[tot].z=-1;q[tot].belong=T; q[++tot].x=x2;q[tot].y=y2;q[tot].t=tot;q[tot].op=2;q[tot].z=1;q[tot].belong=T; } } sort(q+1,q+tot+1); work(1,tot); for(int i=1;i<=T;i++)printf("%d\n",ans[i]); return 0; }