bzoj 1176(cdq分治)

1176: [Balkan2007]Mokia

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 1838  Solved: 810

[Submit][Status][Discuss]

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
分治。

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<iostream>

using namespace std;

int s,w,cnt,wen;

struct ss

 {

  int pan,x,y,z,mul,zhi;

 }q[210000];

int zhi[2100000];

int ans[200000];

inline int read()

{
char y; int x=0,f=1; y=getchar();
while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;

}

bool cmp(ss xg,ss yg)

 {

  if (xg.z!=yg.z)return xg.z<yg.z;else

  if (xg.x!=yg.x)return xg.x<yg.x;else 

   return xg.y<yg.y;

 }

void change(int now,int sug)

 {

  while (now<=w)

  {

  zhi[now]+=sug; now+=now&-now;
 }

 }

int query(int now)

 {

  int sug=0;

  while (now>0)

  {

  sug+=zhi[now]; now-=now&-now;
 }
return sug;

 }

bool cog(ss xg,ss yg)

 {

  if (xg.x!=yg.x) return xg.x<yg.x;else

  return xg.y<yg.y;

 }

void cdq(int l,int r)

 {

   if (l==r)return;

   int mid=(l+r)/2;

   cdq(l,mid); cdq(mid+1,r);

   sort(q+l,q+mid+1,cog); sort(q+mid+1,q+r+1,cog);

   int tail=l;

   for (int i=mid+1;i<=r;++i)

    if (q[i].pan!=0)
{
while (tail<mid+1 && q[i].x>=q[tail].x) 
{
change(q[tail].y,q[tail].zhi);
++tail;
 }

      int opp=query(q[i].y);

      ans[q[i].pan]+=opp*q[i].mul;
}

   for (int i=l;i<=tail-1;++i)

    {

    change(q[i].y,-q[i].zhi);
}

 }

int main()

{
s=read(); w=read(); cnt=0; wen=0; int pn,ogg=0;

    while (scanf("%d",&pn) && pn!=3)

     {

       ++ogg;

       if (pn==1)

        {

          ++cnt; 

          q[cnt].x=read(); q[cnt].y=read(); q[cnt].zhi=read(); q[cnt].z=ogg;
}else
{
++wen; int xg=read(); int yg=read(); int xt=read(); int yt=read();
if (xg!=1 && yg!=1) ++cnt,q[cnt].pan=wen,q[cnt].x=xg-1,q[cnt].y=yg-1,q[cnt].mul=1,q[cnt].z=ogg,ans[wen]+=s*(xg-1)*(yg-1);
++cnt; q[cnt].pan=wen; q[cnt].x=xt; q[cnt].y=yt; q[cnt].mul=1; ans[wen]+=q[cnt].x*q[cnt].y*s; q[cnt].z=ogg;
   if (xg!=1) ++cnt,q[cnt].pan=wen,q[cnt].x=xg-1,q[cnt].y=yt,q[cnt].mul=-1,ans[wen]-=q[cnt].x*q[cnt].y*s,q[cnt].z=ogg;
   if (yg!=1) ++cnt,q[cnt].pan=wen,q[cnt].x=xt,q[cnt].y=yg-1,q[cnt].mul=-1,ans[wen]-=q[cnt].x*q[cnt].y*s,q[cnt].z=ogg;
}
}
sort(q+1,q+cnt+1,cmp);
cdq(1,cnt);
for (int i=1;i<=wen;++i)
{
printf("%d\n",ans[i]);
}

}