題目：

維護一個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

Sample Output

題解：

又一道用cdq來解決三維偏序的題,三維分別是時間，x，y;

首先肯定要拆掉詢問····我們用sum（x,y）來表示包括該點的左上角的矩形值之和,然後針對一個詢問x1,y1,x2,y2,答案明顯是sum(x1-1,y1-1)-sum(x1-1,y2)-sum(x2,y1-1)+sum(x2,y2),因此可以將詢問拆成對應x,y的四個點,然後就可以先以默認的時間為第一順序,每次cdq分治時兩邊以x為第二順序,最後用樹狀數組維護y的第三順序,計算左邊部分的修改對右邊部分詢問的貢獻從而解決問題

代碼：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cctype>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int N=2e5+5;
const int M=2e6+5;
struct node
{
  int type,x,y,val,pos,id;
}query[N],temp[N];
inline 
 
int R()
{
  char c;int f=0,i=1;
  for(c=getchar();(c<‘0‘||c>‘9‘)&&c!=‘-‘;c=getchar());
  if(c==‘-‘)  i=-1,c=getchar();
  for(;c<=‘9‘&&c>=‘0‘;c=getchar())
    f=(f<<3)+(f<<1)+c-‘0‘;
  return f*i;
}
int w,s,tree[M],n,m,X1,Y1,X2,Y2,ans[10005],tim,tag[M];
bool cmp(node a,node b)
{
  if(b.x==a.x)  return a.id<b.id;
  else return a.x<b.x;
}
inline void insert(int u,int v)
{
  for(int i=u;i<=w;i+=(i&(-i)))
    if(tag[i]!=tim)  tag[i]=tim,tree[i]=v;
    else tree[i]+=v;
}
inline int ask(int u)
{
  int temp=0;
  for(int i=u;i;i-=(i&(-i)))
    if(tag[i]!=tim)  continue;
    else temp+=tree[i];
  return temp;
}
inline void solve(int l,int r)
{
  if(l==r)  return;
  int mid=(l+r)/2;
  solve(l,mid),solve(mid+1,r);
  int i=l,j=mid+1,k=l;tim++;
  while(i<=mid&&j<=r)
  {
    if(cmp(query[i],query[j]))
    {
      if(!query[i].pos)
        insert(query[i].y,query[i].val);
      temp[k++]=query[i++];
    }
    else
    {
      if(query[j].pos)
        ans[query[j].pos]+=ask(query[j].y)*query[j].val;
      temp[k++]=query[j++];
    }
  }
  while(i<=mid)  temp[k++]=query[i++];
  while(j<=r)  
  {
    if(query[j].pos)
      ans[query[j].pos]+=ask(query[j].y)*query[j].val;
    temp[k++]=query[j++];
  }
  for(i=l;i<=r;i++)  query[i]=temp[i];
}
int main()
{
  //freopen("a.in","r",stdin);
  s=R(),w=R();
  while(true)
  {
    query[++n].type=R();
    if(query[n].type==3)  break;
    if(query[n].type==1)
    {
      query[n].x=R(),query[n].y=R(),query[n].val=R(),query[n].id=n;
    }
    else
    {
      int X1=R()-1,Y1=R()-1,X2=R(),Y2=R();m++;
      query[n].x=X1,query[n].y=Y1,query[n].val=1,query[n].pos=m,query[n].id=n;
      query[++n].x=X1,query[n].y=Y2,query[n].val=-1,query[n].pos=m,query[n].id=n;
      query[++n].x=X2,query[n].y=Y1,query[n].val=-1,query[n].pos=m,query[n].id=n;
      query[++n].x=X2,query[n].y=Y2,query[n].val=1,query[n].pos=m,query[n].id=n;
      ans[m]=(X2-X1)*(Y2-Y1)*s;
    }
  }
  solve(1,n);
  for(int i=1;i<=m;i++)
    printf("%d\n",ans[i]);
  return 0;
}

刷題總結——mokia（bzoj1176）