详细解释见小结。http://blog.csdn.net/zmx354/article/details/31740985

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <stack>
#include <map> #pragma comment(linker, "/STACK:1024000000");
#define EPS (1e-8)
#define LL long long
#define ULL unsigned long long
#define _LL __int64
#define _INF 0x3f3f3f3f
#define Mod 9999991
#define lowbit(x) (x&(-x)) using namespace std; const int N = 2049; int a[N][N];
int ax[N][N];
int ay[N][N];
int axy[N][N]; void Updata(int a[N][N],int x,int y,int n,int m,int data)
{
for(; x <= n; x += lowbit(x))
{
for(int t = y; t <= m; t += lowbit(t))
{
a[x][t] += data;
}
}
} int Query(int a[N][N],int x,int y)
{
int sum = 0; while(x > 0)
{
for(int t = y; t > 0; t -= lowbit(t))
{
sum += a[x][t];
}
x -= lowbit(x);
} return sum;
} int sum(int x,int y)
{
return Query(a,x,y)*(x+1)*(y+1) - (y+1)*Query(ax,x,y) - (x+1)*Query(ay,x,y) + Query(axy,x,y);
} void Add(int x1,int y1,int x2,int y2,int n,int m,int data)
{
Updata(a,x1,y1,n,m,data);
Updata(a,x2+1,y1,n,m,-data);
Updata(a,x1,y2+1,n,m,-data);
Updata(a,x2+1,y2+1,n,m,data); Updata(ax,x1,y1,n,m,x1*data);
Updata(ax,x2+1,y1,n,m,-(x2+1)*data);
Updata(ax,x1,y2+1,n,m,-x1*data);
Updata(ax,x2+1,y2+1,n,m,(x2+1)*data); Updata(ay,x1,y1,n,m,y1*data);
Updata(ay,x2+1,y1,n,m,-y1*data);
Updata(ay,x1,y2+1,n,m,-(y2+1)*data);
Updata(ay,x2+1,y2+1,n,m,(y2+1)*data); Updata(axy,x1,y1,n,m,x1*y1*data);
Updata(axy,x2+1,y1,n,m,-(x2+1)*y1*data);
Updata(axy,x1,y2+1,n,m,-x1*(y2+1)*data);
Updata(axy,x2+1,y2+1,n,m,(x2+1)*(y2+1)*data);
} int main()
{
int n,m; char order[10]; int x1,x2,y1,y2,data; memset(a,0,sizeof(a));
memset(ax,0,sizeof(ax));
memset(ay,0,sizeof(ay));
memset(axy,0,sizeof(axy)); while(scanf("%s",order) != EOF)
{
if(order[0] == 'L')
{
scanf("%d %d %d %d %d",&x1,&y1,&x2,&y2,&data); Add(x1,y1,x2,y2,n,m,data);
}
else if(order[0] == 'k')
{
scanf("%d %d %d %d",&x1,&y1,&x2,&y2); printf("%d\n",sum(x2,y2) - sum(x1-1,y2) - sum(x2,y1-1) + sum(x1-1,y1-1));
}
else
{
scanf("%d %d",&n,&m);
}
} return 0;
}

二维树状数组的区间加减及查询 tyvj 1716 上帝造题的七分钟的更多相关文章

  1. POJ 2155 Matrix (二维线段树入门,成段更新,单点查询 / 二维树状数组,区间更新,单点查询)

    题意: 有一个n*n的矩阵,初始化全部为0.有2中操作: 1.给一个子矩阵,将这个子矩阵里面所有的0变成1,1变成0:2.询问某点的值 方法一:二维线段树 参考链接: http://blog.csdn ...

  2. POJ2155【二维树状数组,区间修改,点查询?】【又被输入输出坑】

    这题反反复复,到现在才过. 这道题就是树状数组的逆用,用于修改区间内容,查询点的值. 如果单纯就这个奇偶数来判的话,似乎这个思路比较好理解. 看了一下国家集训队论文(囧),<关于0与1在信息学奥 ...

  3. POJ2155 Matrix 【二维树状数组】+【段更新点查询】

    Matrix Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 17766   Accepted: 6674 Descripti ...

  4. 【bzoj5173】[Jsoi2014]矩形并 扫描线+二维树状数组区间修改区间查询

    题目描述 JYY有N个平面坐标系中的矩形.每一个矩形的底边都平行于X轴,侧边平行于Y轴.第i个矩形的左下角坐标为(Xi,Yi),底边长为Ai,侧边长为Bi.现在JYY打算从这N个矩形中,随机选出两个不 ...

  5. 牛客网 暑期ACM多校训练营(第二场)J.farm-STL(vector)+二维树状数组区间更新、单点查询 or 大暴力?

    开心.jpg J.farm 先解释一下题意,题意就是一个n*m的矩形区域,每个点代表一个植物,然后不同的植物对应不同的适合的肥料k,如果植物被撒上不适合的肥料就会死掉.然后题目将每个点适合的肥料种类( ...

  6. SCOI2014 bzoj3594 方伯伯的玉米田(二维树状数组+dp)

    3594: [Scoi2014]方伯伯的玉米田 Time Limit: 60 Sec  Memory Limit: 128 MBSubmit: 1971  Solved: 961[Submit][St ...

  7. 【poj2155】Matrix(二维树状数组区间更新+单点查询)

    Description Given an N*N matrix A, whose elements are either 0 or 1. A[i, j] means the number in the ...

  8. POJ 2155 Matrix【二维树状数组+YY(区间计数)】

    题目链接:http://poj.org/problem?id=2155 Matrix Time Limit: 3000MS   Memory Limit: 65536K Total Submissio ...

  9. POJ2155 Matrix(二维树状数组||区间修改单点查询)

    Given an N*N matrix A, whose elements are either 0 or 1. A[i, j] means the number in the i-th row an ...

随机推荐

  1. 使用bottle进行web开发(2):http request

    我们知道,http request有多个方法,比如get,post,delete,patch,put等.对用的,bottle都定义了相应的装饰器,目前定义了五个: get(),post(),put() ...

  2. hdu 5158(水题)

    Have meal Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total S ...

  3. Android 轻松实现后台搭建+APP版本更新

    http://blog.csdn.net/u012422829/article/details/46355515 (本文讲解了在Android中实现APP版本更新,文末附有源码.) 看完本文,您可以学 ...

  4. csu1812

    csu1812 题意 求三角形和矩形交的面积. 分析 半平面交.把三角形的三条边当作直线去切割矩形,最后求切割后的多边形面积即可. code #include<algorithm> #in ...

  5. 拓扑排序【p1137】 旅行计划

    Description 小明要去一个国家旅游.这个国家有\(N\)个城市,编号为\(1\)至\(N\),并且有\(M\)条道路连接着,小明准备从其中一个城市出发,并只往东走到城市\(i\)停止. 所以 ...

  6. 数学【P2524】 Uim的情人节礼物·其之弐 (康托展开)

    因为某人@ZAGER挖坑让我讲一下康托展开,所以发现了这个题,顺便说一下康托展开是个什么东西 题目概括 给定n与一个数列,要求求出给定数列在n的全排列中的排名(按照字典序从小到大排列) 康托展开 先放 ...

  7. Frequency

    题目描述 Snuke loves constructing integer sequences. There are N piles of stones, numbered 1 through N. ...

  8. 三. Java类与对象9. 源文件的声明规则

    当在一个源文件中定义多个类,并且还有import语句和package语句时,要特别注意这些规则: 一个源文件中只能有一个public类. 一个源文件可以有多个非public类. 源文件的名称应该和pu ...

  9. Linux下Shell的复制粘贴快捷键

    [Shift]+[Insert]:复制 [Ctrl]+[Insert]:粘贴

  10. 【SQL】用Sql Server自动生产html格式的数据字典

    原文:[SQL]用Sql Server自动生产html格式的数据字典 本文软件环境:Sql Server 2008. 1.打开sql server管理器,给选定的表添加描述信息,给指定的字段添加描述信 ...