详细解释见小结。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. Java微信公众平台开发_01_本地服务器映射外网

    做微信开发之前,我们需要先做一个内网穿透,让我们的工程可以在公网上被访问. 一.工具列表 内网穿透的相关工具有: (1)natapp 官网 :https://natapp.cn/ (2)花生壳 官网: ...

  2. Educational Codeforces Round 33 (Rated for Div. 2) C. Rumor【并查集+贪心/维护集合最小值】

    C. Rumor time limit per test 2 seconds memory limit per test 256 megabytes input standard input outp ...

  3. 学习sphinx

    在我们php开发中如何使用sphinx技术? 答:对应MySQL数据库中的字段,就是将字段中的中文信息拆分成多个词语,然后对这些词语建立索引.以后查询的时候,先去查询这些索引文件,然后返回这些满足条件 ...

  4. 以最简单的方式了解--Github

    大概是从寒假的时候开始正式的赚取github,从github上面学习一些开源的文档,我记得我注册github账号到现在已经9个月了,但只有最近的2个月才发现github这个新世界,写这篇文章是为了刚入 ...

  5. 网络爬虫框架Webmagic

    1 谈谈网络爬虫 1.1 什么是网络爬虫 在大数据时代,信息的采集是一项重要的工作,而互联网中的数据是海量的,如果单纯靠人力进行信息采集,不仅低效繁琐,搜集的成本也会提高.如何自动高效地获取互联网中我 ...

  6. 警告Conversion specifies type'int' but the argument has type'size_t'

    代码: #import<Foundation/Foundation.h> int main(int argc,const char * argv[]){ const char *words ...

  7. 来自Jakob Jenkov的Jackson教程

    Jakob Jenkov是Java界的牛人. 下面是收集的与Jackson有关的教程: JSON: http://tutorials.jenkov.com/java-json/index.html J ...

  8. source insight研究——快捷键篇

    转:http://blog.csdn.net/ison81/article/details/3510426 关于键盘和鼠标谁更快捷之争,是一个永远被程序员争论的话题.我想大多数人都不会极端的信奉一种操 ...

  9. 利用Yum彻底移除docker

    yum remove docker \ docker-client \ docker-client-latest \ docker-common \ docker-latest \ docker-la ...

  10. jquery 遍历取值问题

    方法一:将a定义在外边 var a = new Array; var i = 0 ; $("img").each(function(){ a[i] = $(this).attr(& ...