Fast Matrix Operations (UVA - 11992,二维线段树 + 简单应用)

2024-03-30 14:08

本文主要是介绍Fast Matrix Operations (UVA - 11992,二维线段树 + 简单应用),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一.题目链接:

UVA-11992

二.题目大意:

有一个 r × c 大小的矩阵,最多有 10^{6} 个元素,r ≤ 20.

矩阵中元素的初始值均为 0.

现有三种操作

① 1 x1 y1 x2 y2 v :将 (x1, y1) 到 (x2, y2) 之间的元素都加上值 v.

② 2 x1 y1 x2 y2 v:将 (x1, y1) 到 (x2, y2) 之间的元素都置为值 v.

③ 3 x1 y1 x2 y2  :查询 (x1, y1) 到 (x2, y2) 之间的元素的 和、最大值、最小值,并将其输出.

操作总数 m ≤ 10^{6}.

三.分析:

这三种操作当然是线段树的基本操作啦.

由于 r ≤ 20,所以建造 r 个线段树即可.

把 左右区间、最小值、最大值、和、懒惰标记、set 值封装在结构体了就可以了.

这里注意一点:

在函数 down 中 set 和 add 的顺序

假设操作时先进行了 set 后进行了 add,那么这里肯定是要先 set 再 add 的.

相反,由于在 set 时,已经将 tree[row][k].add = 0,所以此时也是可以先进行 set的.

四.代码实现:

#include <set>
#include <map>
#include <ctime>
#include <queue>
#include <cmath>
#include <stack>
#include <vector>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define eps 1e-4
#define PI acos(-1.0)
#define ll long long int
using namespace std;const int M = (int)1e6 + 10;
const int inf = 0x3f3f3f3f;
struct node1
{int l, r;int sum, add;int Max, Min;int Set;
}tree[21][M * 3];struct node2
{int sum;int Min;int Max;
}tnode;void build(int k, int l, int r, int row)
{tree[row][k].l = l;tree[row][k].r = r;tree[row][k].sum = 0;tree[row][k].add = 0;tree[row][k].Max = 0;tree[row][k].Min = 0;tree[row][k].Set = -1;if(l == r)return;int m = (l + r) / 2;build(k * 2, l, m, row);build(k * 2 + 1, m + 1, r, row);
}void down(int row, int k)
{int lk = 2 * k;int rk = 2 * k + 1;if(~tree[row][k].Set){tree[row][lk].Set = tree[row][rk].Set = tree[row][k].Set;tree[row][lk].add = tree[row][rk].add = 0;tree[row][k].Set = -1;tree[row][lk].sum = (tree[row][lk].r - tree[row][lk].l + 1) * tree[row][lk].Set;tree[row][rk].sum = (tree[row][rk].r - tree[row][rk].l + 1) * tree[row][rk].Set;tree[row][lk].Max = tree[row][lk].Min = tree[row][lk].Set;tree[row][rk].Max = tree[row][rk].Min = tree[row][rk].Set;}if(tree[row][k].add > 0){tree[row][lk].add += tree[row][k].add;tree[row][rk].add += tree[row][k].add;tree[row][lk].sum += (tree[row][lk].r - tree[row][lk].l + 1) * tree[row][k].add;tree[row][rk].sum += (tree[row][rk].r - tree[row][rk].l + 1) * tree[row][k].add;tree[row][lk].Min += tree[row][k].add;tree[row][lk].Max += tree[row][k].add;tree[row][rk].Min += tree[row][k].add;tree[row][rk].Max += tree[row][k].add;tree[row][k].add = 0;}
}void add(int k, int y1, int y2, int v, int row)
{if(tree[row][k].l >= y1 && tree[row][k].r <= y2){tree[row][k].sum += (tree[row][k].r - tree[row][k].l + 1) * v;tree[row][k].add += v;tree[row][k].Min += v;tree[row][k].Max += v;return;}if(tree[row][k].add || ~tree[row][k].Set)down(row, k);int m  = (tree[row][k].l + tree[row][k].r) / 2;if(y1 <= m)add(k * 2, y1, y2, v, row);if(m < y2)add(k * 2 + 1, y1, y2, v, row);tree[row][k].sum = tree[row][k * 2].sum + tree[row][k * 2 + 1].sum;tree[row][k].Max = max(tree[row][k * 2].Max, tree[row][k * 2 + 1].Max);tree[row][k].Min = min(tree[row][k * 2].Min, tree[row][k * 2 + 1].Min);
}void Set(int k, int y1, int y2, int v, int row)
{if(tree[row][k].l >= y1 && tree[row][k].r <= y2){tree[row][k].sum = (tree[row][k].r - tree[row][k].l + 1) * v;tree[row][k].add = 0;tree[row][k].Set = v;tree[row][k].Min = v;tree[row][k].Max = v;return;}if(tree[row][k].add || ~tree[row][k].Set)down(row, k);int m  = (tree[row][k].l + tree[row][k].r) / 2;if(y1 <= m)Set(k * 2, y1, y2, v, row);if(m < y2)Set(k * 2 + 1, y1, y2, v, row);tree[row][k].sum = tree[row][k * 2].sum + tree[row][k * 2 + 1].sum;tree[row][k].Max = max(tree[row][k * 2].Max, tree[row][k * 2 + 1].Max);tree[row][k].Min = min(tree[row][k * 2].Min, tree[row][k * 2 + 1].Min);
}void query(int k, int y1, int y2, int row)
{if(tree[row][k].l >= y1 && tree[row][k].r <= y2){tnode.sum += tree[row][k].sum;tnode.Max = max(tnode.Max, tree[row][k].Max);tnode.Min = min(tnode.Min, tree[row][k].Min);return;}if(tree[row][k].add || ~tree[row][k].Set)down(row, k);int m = (tree[row][k].l + tree[row][k].r) / 2;if(y1 <= m)query(k * 2, y1, y2, row);if(m < y2)query(k * 2 + 1, y1, y2, row);
}int main()
{int r, c, m;while(~scanf("%d %d %d", &r, &c, &m)){for(int i = 1; i <= r; ++i)build(1, 1, c, i);while(m--){int t, x1, y1, x2, y2;scanf("%d %d %d %d %d", &t, &x1, &y1, &x2, &y2);if(t == 1){int v;scanf("%d", &v);for(int i = x1; i <= x2; ++i)add(1, y1, y2, v, i);}else if(t == 2){int v;scanf("%d", &v);for(int i = x1; i <= x2; ++i)Set(1, y1, y2, v, i);}else if(t == 3){int sum = 0;int Max = -inf;int Min = inf;for(int i = x1; i <= x2; ++i){tnode.sum = 0;tnode.Max = -inf;tnode.Min = inf;query(1, y1, y2, i);sum += tnode.sum;Max = max(Max, tnode.Max);Min = min(Min, tnode.Min);}printf("%d %d %d\n", sum, Min, Max);}}}return 0;
}

 

这篇关于Fast Matrix Operations (UVA - 11992,二维线段树 + 简单应用)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/861690

相关文章

中文分词jieba库的使用与实景应用(一)

知识星球:https://articles.zsxq.com/id_fxvgc803qmr2.html 目录 一.定义: 精确模式(默认模式): 全模式: 搜索引擎模式: paddle 模式(基于深度学习的分词模式): 二 自定义词典 三.文本解析   调整词出现的频率 四. 关键词提取 A. 基于TF-IDF算法的关键词提取 B. 基于TextRank算法的关键词提取

水位雨量在线监测系统概述及应用介绍

在当今社会,随着科技的飞速发展,各种智能监测系统已成为保障公共安全、促进资源管理和环境保护的重要工具。其中,水位雨量在线监测系统作为自然灾害预警、水资源管理及水利工程运行的关键技术,其重要性不言而喻。 一、水位雨量在线监测系统的基本原理 水位雨量在线监测系统主要由数据采集单元、数据传输网络、数据处理中心及用户终端四大部分构成,形成了一个完整的闭环系统。 数据采集单元:这是系统的“眼睛”,

poj2576(二维背包)

题意:n个人分成两组,两组人数只差小于1 , 并且体重只差最小 对于人数要求恰好装满,对于体重要求尽量多,一开始没做出来,看了下解题,按照自己的感觉写,然后a了 状态转移方程:dp[i][j] = max(dp[i][j],dp[i-1][j-c[k]]+c[k]);其中i表示人数,j表示背包容量,k表示输入的体重的 代码如下: #include<iostream>#include<

hdu2159(二维背包)

这是我的第一道二维背包题,没想到自己一下子就A了,但是代码写的比较乱,下面的代码是我有重新修改的 状态转移:dp[i][j] = max(dp[i][j], dp[i-1][j-c[z]]+v[z]); 其中dp[i][j]表示,打了i个怪物,消耗j的耐力值,所得到的最大经验值 代码如下: #include<iostream>#include<algorithm>#include<

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

zoj3820(树的直径的应用)

题意:在一颗树上找两个点,使得所有点到选择与其更近的一个点的距离的最大值最小。 思路:如果是选择一个点的话,那么点就是直径的中点。现在考虑两个点的情况,先求树的直径,再把直径最中间的边去掉,再求剩下的两个子树中直径的中点。 代码如下: #include <stdio.h>#include <string.h>#include <algorithm>#include <map>#