UVa 11992 Fast Matrix Operations 线段树

2024-09-05 15:58

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

UVa 11992 Fast Matrix Operations

题目大意:有一个r行c列的全0矩阵,支持三种操作:

1 x1 y1 x2 y2 v 子矩阵(x1,y1,x2,y2)的所有元素增加v(v > 0)。

2 x1 y1 x2 y2 v 子矩阵(x1,y1,x2,y2)的所有元素设为v(v > 0)。

3 x1 y1 x2 y2    查询子矩阵(x1,y1,x2,y2)的元素和、最小值和最大值。

子矩阵(x1,y1,x2,y2)是指满足x1 <= x <= x2,y1 <= y <= y2,的所有元素(x,y)。输入保证任何时刻矩阵内的所有元素之和不超过10 ^ 9。

其中,矩阵不超过20行,元素总数不超过10 ^ 6,且操作数不超过20000。

题目分析:

本题即区间修改,区间查询问题。

本来是想写一棵树套树的。。不过发现实在是写不出来,不会写QAQ,然后转头扑倒一维线段树。

矩阵不超过20行,所以我们可以为每行建一棵线段树,则本题便可以化为一维线段树。

其他的就注意一下值的传递即可。


代码如下:


#include <cstdio>
#include <cstring>#define clear( A , X ) memset ( A , X , sizeof A )
#define lson l , m , o << 1
#define rson m + 1 , r , o << 1 | 1const int maxR = 22 ;
const int maxC = 200005 ;
const int oo = 0x3f3f3f3f ;int max[maxR][maxC] , min[maxR][maxC] , sum[maxR][maxC] ;
int add[maxR][maxC] , set[maxR][maxC] ;int Min ( const int X , const int Y ) {if ( X < Y ) return X ;return Y ;
}int Max ( const int X , const int Y ) {if ( X > Y ) return X ; return Y ;
}void PushUp ( int x , int o ) {sum[x][o] = sum[x][o << 1] + sum[x][o << 1 | 1] ;min[x][o] = Min ( min[x][o << 1] , min[x][o << 1 | 1] ) ;max[x][o] = Max ( max[x][o << 1] , max[x][o << 1 | 1] ) ;
}void PushDown ( int x , int l , int r , int o ) {int m = ( l + r ) >> 1 ;if ( set[x][o] ) {int tmp = set[x][o] ;set[x][o << 1] = tmp ;min[x][o << 1] = tmp ;max[x][o << 1] = tmp ;sum[x][o << 1] = tmp * ( m - l + 1 ) ;add[x][o << 1] =   0 ;set[x][o << 1 | 1] = tmp ;min[x][o << 1 | 1] = tmp ;max[x][o << 1 | 1] = tmp ;sum[x][o << 1 | 1] = tmp * ( r - m ) ;add[x][o << 1 | 1] =   0 ;set[x][o] = 0 ;}if ( add[x][o] ) {int tmp = add[x][o] ;add[x][o << 1] += tmp ;min[x][o << 1] += tmp ;max[x][o << 1] += tmp ;sum[x][o << 1] += tmp * ( m - l + 1 ) ;add[x][o << 1 | 1] += tmp ;min[x][o << 1 | 1] += tmp ;max[x][o << 1 | 1] += tmp ;sum[x][o << 1 | 1] += tmp * ( r - m ) ;add[x][o] = 0 ;}
}void Build ( int x , int l , int r , int o ) {set[x][o] = add[x][o] = 0 ;if ( l == r ) {min[x][o] = max[x][o] = sum[x][o] = 0 ;return ;}int m = ( l + r ) >> 1 ;Build ( x , lson ) ;Build ( x , rson ) ;
}void Set ( int L , int R , int v , int x , int l , int r , int o ) {if ( L <= l && r <= R ) {add[x][o] = 0 ;set[x][o] = v ;min[x][o] = v ;max[x][o] = v ;sum[x][o] = v * ( r - l + 1 ) ;return ;}PushDown ( x , l , r , o ) ;int m = ( l + r ) >> 1 ;if ( L <= m ) Set ( L , R , v , x , lson ) ;if ( m <  R ) Set ( L , R , v , x , rson ) ;PushUp ( x , o ) ;
}void Add ( int L , int R , int v , int x , int l , int r , int o ) {if ( L <= l && r <= R ) {add[x][o] += v ;min[x][o] += v ;max[x][o] += v ;sum[x][o] += v * ( r - l + 1 ) ;return ;}PushDown ( x , l , r , o ) ;int m = ( l + r ) >> 1 ;if ( L <= m ) Add ( L , R , v , x , lson ) ;if ( m <  R ) Add ( L , R , v , x , rson ) ;PushUp ( x , o ) ;
}int ansmin , ansmax , anssum ;void Query ( int L , int R , int x , int l , int r , int o ) {if ( L <= l && r <= R ) {ansmin = Min ( ansmin , min[x][o] ) ;ansmax = Max ( ansmax , max[x][o] ) ;anssum += sum[x][o] ;return ;}PushDown ( x , l , r , o ) ;int m = ( l + r ) >> 1 ;if ( L <= m ) Query ( L , R , x , lson ) ;if ( m <  R ) Query ( L , R , x , rson ) ;
}void work () {int n , m , q , ch , x1 , x2 , y1 , y2 , v ;while ( ~scanf ( "%d%d%d" , &n , &m , &q ) ) {for ( int i = 1 ; i <= n ; ++ i ) Build ( i , 1 , m , 1 ) ;while ( q -- ) {scanf ( "%d" , &ch ) ;if ( 1 == ch ) {scanf ( "%d%d%d%d%d" , &x1 , &y1 , &x2 , &y2 , &v ) ;for ( int i = x1 ; i <= x2 ; ++ i ) {Add ( y1 , y2 , v , i , 1 , m , 1 ) ;}}if ( 2 == ch ) {scanf ( "%d%d%d%d%d" , &x1 , &y1 , &x2 , &y2 , &v ) ;for ( int i = x1 ; i <= x2 ; ++ i ) {Set ( y1 , y2 , v , i , 1 , m , 1 ) ;}}if ( 3 == ch ) {scanf ( "%d%d%d%d" , &x1 , &y1 , &x2 , &y2 ) ;ansmin =  oo ;ansmax =   0 ;anssum =   0 ;for ( int i = x1 ; i <= x2 ; ++ i ) {Query ( y1 , y2 , i , 1 , m , 1 ) ;}printf ( "%d %d %d\n" , anssum , ansmin , ansmax ) ;}}}
}int main () {work () ;return 0 ;
}


这篇关于UVa 11992 Fast Matrix Operations 线段树的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

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

uva 10055 uva 10071 uva 10300(水题两三道)

情歌两三首,水题两三道。 好久没敲代码了为暑假大作战热热身。 uva 10055 Hashmat the Brave Warrior 求俩数相减。 两个debug的地方,一个是longlong,一个是输入顺序。 代码: #include<stdio.h>int main(){long long a, b;//debugwhile(scanf("%lld%lld", &

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

uva 10061 How many zero's and how many digits ?(不同进制阶乘末尾几个0)+poj 1401

题意是求在base进制下的 n!的结果有几位数,末尾有几个0。 想起刚开始的时候做的一道10进制下的n阶乘末尾有几个零,以及之前有做过的一道n阶乘的位数。 当时都是在10进制下的。 10进制下的做法是: 1. n阶位数:直接 lg(n!)就是得数的位数。 2. n阶末尾0的个数:由于2 * 5 将会在得数中以0的形式存在,所以计算2或者计算5,由于因子中出现5必然出现2,所以直接一

uva 568 Just the Facts(n!打表递推)

题意是求n!的末尾第一个不为0的数字。 不用大数,特别的处理。 代码: #include <stdio.h>const int maxn = 10000 + 1;int f[maxn];int main(){#ifdef LOCALfreopen("in.txt", "r", stdin);#endif // LOCALf[0] = 1;for (int i = 1; i <=

uva 575 Skew Binary(位运算)

求第一个以(2^(k+1)-1)为进制的数。 数据不大,可以直接搞。 代码: #include <stdio.h>#include <string.h>const int maxn = 100 + 5;int main(){char num[maxn];while (scanf("%s", num) == 1){if (num[0] == '0')break;int len =