【bzoj1176】[Balkan2007]Mokia

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分治。
对于询问，可以拆分成对四个以原点为左下角矩形求值。
对于操作和询问按照x坐标从小到大排序。
按照操作先后分l - mid - r区间
计算操作次序在l-mid对mid+r的贡献
然后递归处理[l,mid],[mid+1,r]

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
int s,w,n;
int ans[10005];
int t[2000005];
struct node{int x,y,a,pos,id;bool flag;}q[200005],tmp[200005];
bool cmp(node a,node b)
{
    if (a.x!=b.x) return a.x<b.x;
    if (a.y!=b.y) return a.y<b.y;
    return a.flag<b.flag;
}
inline void Add()
{
    int X1=read(),Y1=read(),X2=read(),Y2=read();
    int pos=++ans[0];
    q[++n].pos=pos;q[n].id=n;q[n].x=X1-1;q[n].y=Y1-1;q[n].a=1;q[n].flag=1;
    q[++n].pos=pos;q[n].id=n;q[n].x=X2;q[n].y=Y2;q[n].a=1;q[n].flag=1;
    q[++n].pos=pos;q[n].id=n;q[n].x=X1-1;q[n].y=Y2;q[n].a=-1;q[n].flag=1;
    q[++n].pos=pos;q[n].id=n;q[n].x=X2;q[n].y=Y1-1;q[n].a=-1;q[n].flag=1;
}
void add(int x,int y)
{
    for (;x<=w;x+=x&-x)t[x]+=y;
}
int query(int x)
{
    int ans=0;
    for (;x;x-=x&-x)ans+=t[x];
    return ans;
}
void cdq(int l,int r)
{
    if (l==r) return;
    int mid=(l+r)>>1,l1=l,l2=mid+1;
    for (int i=l;i<=r;i++)
    {
        if (q[i].id<=mid&&!q[i].flag) add(q[i].y,q[i].a);
        if (q[i].id>mid&&q[i].flag) ans[q[i].pos]+=query(q[i].y)*q[i].a;
    }
    for (int i=l;i<=r;i++)
        if (q[i].id<=mid&&!q[i].flag) add(q[i].y,-q[i].a);
    for (int i=l;i<=r;i++)
        if (q[i].id<=mid) tmp[l1++]=q[i];else tmp[l2++]=q[i];
    for (int i=l;i<=r;i++) q[i]=tmp[i];
    cdq(l,mid);cdq(mid+1,r);
}
int main()
{
    s=read();w=read();
    while (1)
    {
        int opt=read();
        if (opt==1){q[++n].x=read();q[n].y=read();q[n].a=read();q[n].id=n;}
        else if (opt==2) Add();
        else break;
    }
    sort(q+1,q+n+1,cmp);
    cdq(1,n);
    for (int i=1;i<=ans[0];i++) printf("%d\n",ans[i]);
    return 0;
}