[bzoj1176][cdq分治]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裸题
一个操作拆成四个随便做.
修改放在询问前面记得

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
struct ask{int x,y,gg;LL opt;}q[710000],tt[710000];int ln;
LL s[2100005],S;
int n;
int lowbit(int x){return x&-x;}
void change(int x,LL c){while(x<=n)s[x]+=c,x+=lowbit(x);}
LL findsum(int x){LL ret=0;while(x>=1)ret+=s[x],x-=lowbit(x);return ret;}
LL ans[11000],tlen;
//bit
void cdq(int l,int r)
{
    if(l==r)return ;
    int mid=(l+r)/2;cdq(l,mid);cdq(mid+1,r);
    int u=l,v=mid+1,ln=0;
    while(u<=mid && v<=r)
    {
        if(q[u].x<=q[v].x)
        {
            if(!q[u].gg)change(q[u].y,q[u].opt);
            tt[++ln]=q[u++];
        }
        else
        {   
            if(q[v].gg)ans[q[v].gg]+=q[v].opt*findsum(q[v].y);
            tt[++ln]=q[v++];
        }
    }
    while(u<=mid)
    {
        if(!q[u].gg)change(q[u].y,q[u].opt);
        tt[++ln]=q[u++];
    }
    while(v<=r)
    {
        if(q[v].gg)ans[q[v].gg]+=q[v].opt*findsum(q[v].y);
        tt[++ln]=q[v++];
    }
    for(int i=l;i<=mid;i++)if(!q[i].gg)change(q[i].y,-q[i].opt);
    for(int i=l;i<=r;i++)q[i]=tt[i-l+1];
}
int main()
{
//  freopen("10.in","r",stdin);
//  freopen("a.out","w",stdout);
    scanf("%d%d",&S,&n);
    int cnt;
    while(scanf("%d",&cnt)!=EOF)
    {
        if(cnt==3)break;
        int x1,y1,x2,y2;scanf("%d%d%d",&x1,&y1,&x2);
        if(cnt==2)
        {
            scanf("%d",&y2);tlen++;
            q[++ln].x=x2;q[ln].y=y2;q[ln].gg=tlen;q[ln].opt=1;
            if(x1-1>=1)q[++ln].x=x1-1,q[ln].y=y2,q[ln].gg=tlen,q[ln].opt=-1;
            if(y1-1>=1)q[++ln].x=x2,q[ln].y=y1-1,q[ln].gg=tlen,q[ln].opt=-1;
            if(x1-1>=1 && y1-1>=1)q[++ln].x=x1-1,q[ln].y=y1-1,q[ln].gg=tlen,q[ln].opt=1;
            ans[tlen]=(LL)(x2-x1+1)*(y2-y1+1)*S;
        }
        else q[++ln].x=x1,q[ln].y=y1,q[ln].opt=(LL)x2;
    }
    cdq(1,ln);
    for(int i=1;i<=tlen;i++)printf("%lld\n",ans[i]);
    return 0;
}

[bzoj1176][cdq分治]Mokia

猜你喜欢