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=2000100;
inline int read(){
    int sum=0,flag=1;char ch=getchar();
    while(ch!='-'&&(ch>'9'||ch<'0'))ch=getchar();
    if(ch=='-'){flag=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){sum=sum*10+ch-'0';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=0,ict=0;
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=0;
    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)>>1;
    int l1=l,l2=mid+1;
    for(int i=l;i<=r;i++){
        if(a[i].flag==1 && a[i].ti<=mid)add(a[i].y,a[i].a);
        else if(a[i].flag==2 && 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==1 && 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+1,r);
    return;
}
int main(){
    int i,j,S;
    S=read();n=read();
    int op,X1,Y1,X2,Y2,w;
    while(1){
        op=read();
        if(op==3)break;
        if(op==1){
            X1=read();Y1=read();w=read();
            a[++cnt].x=X1;a[cnt].y=Y1;a[cnt].a=w;a[cnt].flag=1;a[cnt].ti=cnt;
        }
        else{
            X1=read()-1;Y1=read()-1;X2=read();Y2=read();
            a[++cnt].x=X1;a[cnt].y=Y1;a[cnt].a=1;a[cnt].flag=2;a[cnt].ti=cnt;a[cnt].id=++ict;
            a[++cnt].x=X1;a[cnt].y=Y2;a[cnt].a=-1;a[cnt].flag=2;a[cnt].ti=cnt;a[cnt].id=ict;
            a[++cnt].x=X2;a[cnt].y=Y1;a[cnt].a=-1;a[cnt].flag=2;a[cnt].ti=cnt;a[cnt].id=ict;
            a[++cnt].x=X2;a[cnt].y=Y2;a[cnt].a=1;a[cnt].flag=2;a[cnt].ti=cnt;a[cnt].id=ict;
            ans[ict]=(X2-X1)*(Y2-Y1)*S;
        }
    }
    sort(a+1,a+cnt+1,cmp);
    solve(1,cnt);
    for(i=1;i<=ict;i++){
        printf("%lld\n",ans[i]);
    }
    return 0;
}