Bzoj1176 [Balkan2007]Mokia

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 2000  Solved: 890

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

CDQ分治

 #include<iostream>
 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 #include<algorithm>
 #define LL long long
 using namespace std;
 const int mxn=;
 inline int read(){
     int sum=,flag=;char ch=getchar();
     while(ch!='-'&&(ch>''||ch<''))ch=getchar();
     if(ch=='-'){flag=-;ch=getchar();}
     while(ch<=''&&ch>=''){sum=sum*+ch-'';ch=getchar();}
     return sum*flag;
 }
 int n;
 int t[mxn];
 struct opt{
     int flag,ti;int x,y,a;int id;
 }a[mxn];
 int cnt=,ict=;
 LL ans[mxn];
 int cmp(opt a,opt b){
     return (a.x<b.x || (a.x==b.x && a.ti<b.ti));
 }
 inline void add(int x,int v){
     while(x<=n){t[x]+=v;x+=x&-x;}
 }
 inline int ask(int x){
     int res=;
     while(x){res+=t[x];x-=x&-x;}
     return res;
 }
 opt p[mxn];
 void solve(int l,int r){
     if(l>=r)return;
     int mid=(l+r)>>;
     int l1=l,l2=mid+;
     for(int i=l;i<=r;i++){
         if(a[i].flag== && a[i].ti<=mid)add(a[i].y,a[i].a);
         else if(a[i].flag== && a[i].ti>mid) ans[a[i].id]+=ask(a[i].y)*a[i].a;
     }
     for(int i=l;i<=r;i++)if(a[i].flag== && a[i].ti<=mid)add(a[i].y,-a[i].a);
     for(int i=l;i<=r;i++){
         if(a[i].ti<=mid) p[l1++]=a[i];
         else p[l2++]=a[i];
     }
     for(int i=l;i<=r;i++)a[i]=p[i];
     solve(l,mid);solve(mid+,r);
     return;
 }
 int main(){
     int i,j,S;
     S=read();n=read();
     int op,X1,Y1,X2,Y2,w;
     while(){
         op=read();
         if(op==)break;
         if(op==){
             X1=read();Y1=read();w=read();
             a[++cnt].x=X1;a[cnt].y=Y1;a[cnt].a=w;a[cnt].flag=;a[cnt].ti=cnt;
         }
         else{
             X1=read()-;Y1=read()-;X2=read();Y2=read();
             a[++cnt].x=X1;a[cnt].y=Y1;a[cnt].a=;a[cnt].flag=;a[cnt].ti=cnt;a[cnt].id=++ict;
             a[++cnt].x=X1;a[cnt].y=Y2;a[cnt].a=-;a[cnt].flag=;a[cnt].ti=cnt;a[cnt].id=ict;
             a[++cnt].x=X2;a[cnt].y=Y1;a[cnt].a=-;a[cnt].flag=;a[cnt].ti=cnt;a[cnt].id=ict;
             a[++cnt].x=X2;a[cnt].y=Y2;a[cnt].a=;a[cnt].flag=;a[cnt].ti=cnt;a[cnt].id=ict;
             ans[ict]=(X2-X1)*(Y2-Y1)*S;
         }
     }
     sort(a+,a+cnt+,cmp);
     solve(,cnt);
     for(i=;i<=ict;i++){
         printf("%lld\n",ans[i]);
     }
     return ;
 }