详细解释见小结。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. MFC学习之EDIT控件初始化

    //四种方法为EDIT控件初始化     //调用系统API     HWND hEidt = ::GetDlgItem(m_hWnd,IDC_EDIT1);     ::SetWindowText( ...

  2. Appium+python自动化7-输入中文【转载】

    前言 在做app自动化过程中会踩很多坑,咱们都是用的中文的app,所以首先要解决中文输入的问题! 本篇通过屏蔽软键盘,绕过手机的软键盘方法,解决中文输入问题. 一.定位搜索 1.打开淘宝点搜索按钮,进 ...

  3. input 限制输入字段总结

    把输入框中 输入的字符串含有中文逗号 改成 英文逗号 举例: <input type="text" id="keywords" style="w ...

  4. PHP获取不带后缀的文件名方法

    $filename = "test.txt"; $houzhui = substr(strrchr($filename, '.'), 1); $result = basename( ...

  5. [BZOJ1295][SCOI2009]最长距离 最短路+枚举

    1295: [SCOI2009]最长距离 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1683  Solved: 912[Submit][Statu ...

  6. SpringMVC + Hibernate + MySQL 的简易网页搭建(Dao层 实现篇)

    首先在搭建一个网站后台前,需要明确出一个合理的网页搭建的结构框架,即从DB  ==>  Dao层 ==>Service层 ==>Control层 ==>View层(该层严格意义 ...

  7. HDU 4417 Super Mario (划分树)(二分)

    Super Mario Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total ...

  8. 四. Java继承和多态4. 多态和动态绑定

    在Java中,父类的变量可以引用父类的实例,也可以引用子类的实例. 请读者先看一段代码: public class Demo { public static void main(String[] ar ...

  9. iOS 代理 重定向消息 forwardInvocation

    今天简单研究一下iOS的重定向消息forwardInvocation: 首先看看Invocation类: @interface NSInvocation : NSObject { @private _ ...

  10. hash算法散列算法

    Hash,一般翻译做“散列”,也有直接音译为“哈希”的,就是把任意长度的输入(又叫做预映射, pre-image),通过散列算法,变换成固定长度的输出,该输出就是散列值.这种转换是一种压缩映射,也就是 ...