hoj 1640 Mobile phones //poj 1195 Mobile phones 二维树状数组

/*

(x1,y2)   ____________    (x2,y2)

     　　 |                      |

     　　 |                      |

      　　|                      |

     　　 |____________|

(x1,y1)            　　　　　　(x2,y1)

 

在上面的矩形中，本题用树状数组来计算的，所以sum(x2,y2)计算的是它的左下方所有的

手机用户数，但是它多算了sum(x2,y1)以下的以及sum(x1,y2)以下的，若减掉即

ans = sum(x2,y2) - sum(x2,y1) -sum(x1,y2)的话，又减多了sum(x1,y1)部分，那只好加

上这一部分了。

 

另外，在update()，sum()时，无非算多了一维，很容易理解的，具体看代码，新学到了

      i -= lowbit(i)   ==      i = i&(i-1)

      i += lowbit(i)   ==      i = (i|(i-1))+1

输入的数据大时速度快很多，比如这题未改前我的代码在hoj上要1.39s，改完后要1.24s

 

*/

#include <iostream>

#include <cstdio>

using namespace std;

const int X = 1026;

int c[X][X],x,y,n,a,q,row,col;

/*int lowbit(int xx)

{

   return xx & -xx;

}*/

int sum(int i,int j)

{

   int temp,ans = 0;

   while(i>0)

   {

      temp = j;

      while(temp>0)

      {

         ans += c[i][temp];

         //temp -= lowbit(temp);

         temp = temp&(temp-1);

      }

      //i -= lowbit(i);

      i = i&(i-1);

   }

   return ans;

}

void update(int i,int j,int num)

{

   int temp;

   while(i<=row)

   {

      temp = j;

      while(temp<=col)

      {

         c[i][temp] += num;

         //temp += lowbit(temp);

         temp = (temp|(temp-1))+1;

      }

      //i += lowbit(i);

      i = (i|(i-1))+1;

   }

}

int main()

{

   freopen("sum.in","r",stdin);

   freopen("sum.out","w",stdout);

   int x1,x2,y1,y2;

   while(scanf("%d",&q)!=EOF)     //多case输入

   {

      if(!q)                //指令为0时，重新载入矩形大小，数组c[][]清零

      {

         scanf("%d",&n);

         row = col = n;

         for(int i=0;i<=n;i++)

            for(int j=0;j<=n;j++)

                c[i][j] = 0;

      }

      else if(q==1)      //指令为1时，在点(x,y)加上新的数，即update(x,y,num)

      {

         scanf("%d%d%d",&x,&y,&a);

         x++;        //给x,y增加一，以免为0时出现死循环

         y++;

         update(x,y,a);

      }

      else if(q==2)      //指令为2时，计算并输出矩形(x1,y1) (x2,y2)上所有的手机数量(大概意思)

      {

         scanf("%d%d%d%d",&x1,&y1,&x2,&y2);

         x2++;

         y2++;

         printf("%d\n",sum(x2,y2)+sum(x1,y1)-sum(x1,y2)-sum(x2,y1));

      }

   }

   return 0;

}