bzoj1176 [Balkan2007]Mokia

http://www.elijahqi.win/2018/01/14/bzoj1176-balkan2007mokia/ ‎

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范围

因为cdq分治我现在还只会用来代替数据结构就暂时先划分到数据结构里了

三维偏序的题 常见可以用cdq来搞 时间+x+y 其中我选择y坐标用树状数组维护

树状数组求前缀和方便所以我容斥一下将整个矩形区间划分成4个区间加加减减即可 关于cdq我当时的学习笔记都写在这里了

http://www.elijahqi.win/2017/12/01/%e5%9e%83%e5%9c%be%e6%98%be%e7%a4%ba%e5%b1%8f/

重新复习cdq的时候发现需要注意几点 1、我永远只考虑左边对右边的影响 即左边右边内部在归并的时候不需要考虑对答案的贡献2、 发现如果都做完了 在进行归并排序的时候假如因为右边做完而停止的那么我接下来继续接过去的时候就不需要再把我的答案加入树状数组了 因为反正后面也会清空 并且对我后面不可能产生影响了 删去 减小常数

#include<cstdio>
#include<algorithm>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int x,y,v,id,f,op;
}qr[220000],tmp[220000];
int s[2200000],ss,w,ans[11000];
inline void add(int x,int v){while(x<=w){s[x]+=v;x+=x&-x;}}
inline int query(int x){int tmp=0;while(x){tmp+=s[x];x-=x&-x;}return tmp;}
inline void clear(int x){
    for (;x<=w;x+=x&-x) if (!s[x]) return;else s[x]=0;
}
inline void cdq(int l,int r){
    if (l==r) return;int mid=l+r>>1;
    cdq(l,mid);cdq(mid+1,r);int h1=l,h2=mid+1,op=l;
    while(h1<=mid&&h2<=r){
        if(qr[h1].x<qr[h2].x||(qr[h1].x==qr[h2].x&&qr[h1].op<qr[h2].op)) {
            if (qr[h1].op==1) add(qr[h1].y,qr[h1].v);tmp[op++]=qr[h1++];
        }else{if (qr[h2].op==2) ans[qr[h2].id]+=qr[h2].f*query(qr[h2].y);tmp[op++]=qr[h2++];}
    }
    while(h1<=mid) tmp[op++]=qr[h1++];
    while(h2<=r) {if (qr[h2].op==2) ans[qr[h2].id]+=qr[h2].f*query(qr[h2].y);tmp[op++]=qr[h2++];}
    for (int i=l;i<=r;++i) {
        qr[i]=tmp[i];if (qr[i].op==1) clear(qr[i].y);
    }
}
int main(){
    freopen("bzoj1176.in","r",stdin);
    ss=read();w=read();int nm=0,num=0;
    while(1){
        int op=read();if (op==3) break;
        if (op==1){
            qr[++num].x=read();qr[num].y=read();qr[num].op=1;qr[num].v=read();
        }else{
            int x1=read(),y1=read(),x2=read(),y2=read();
            ans[++nm]=(x2-x1)*(y2-y1)*ss;
            qr[++num].x=x2;qr[num].y=y2;qr[num].id=nm;qr[num].f=1;qr[num].op=2;
            qr[++num].x=x2;qr[num].y=y1-1;qr[num].id=nm;qr[num].f=-1;qr[num].op=2;
            qr[++num].x=x1-1;qr[num].y=y2;qr[num].id=nm;qr[num].f=-1;qr[num].op=2;
            qr[++num].x=x1-1;qr[num].y=y1-1;qr[num].id=nm;qr[num].f=1;qr[num].op=2;
        }
    }cdq(1,num);
    for (int i=1;i<=nm;++i) printf("%d\n",ans[i]);
    return 0;
}