LA-11992-Fast Matrix Operations 快速矩阵操作（线段树成段更新）

最新推荐文章于 2019-07-09 10:33:44 发布
最新推荐文章于 2019-07-09 10:33:44 发布
阅读量854 收藏
点赞数
分类专栏： 数据结构 树 结构 数据结构 线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/viphong/article/details/48247503
版权
这篇博客介绍了一个使用线段树处理矩阵操作的问题，矩阵大小限制为r*c，r≤20, r*c≤1e6，初始元素全为0。提供了三种操作：1) 子矩阵元素加val；2) 子矩阵元素设为val；3) 查询子矩阵元素的sum、max和min。解决方案是为每一行建立一个线段树，维护sum、max和min信息，并确保在add操作前处理set标记。博主提到代码在VJ上运行时间为600ms，表现良好。" 102646127,8538151,Ubuntu 18.04 安装及配置 Oracle VirtualBox 和 phpVirtualBox,"['Ubuntu虚拟化安装', 'Oracle VirtualBox', 'phpVirtualBox']
摘要由CSDN通过智能技术生成

题意：

给定一个r*c(r<=20,r*c<=1e6)的矩阵，其元素都是0，现在对其子矩阵进行操作。

1 x1 y1 x2 y2 val 表示将(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素add上val；

2 x1 y1 x2 y2 val 表示将(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素set为val；

3 x1 y1 x2 y2 val 表示输出(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素的sum，最大最小值max,min

线段树基本功能啦。。开r个线段树..（最多20）.每次对r个线段树操作就好了

这里有个问题就是，set和add的处理，  要知道set和add不可能同时存在于一个节点中

所以，每次add操作之前都要把当前节点的set标记 推送到 他的儿子中

维护好 sum ，max，min三个信息就可以了

对于查询，写了三个函数，

注意  query_max中 非法区间要返回0，而query_min中非法区间要返回inf 

一开始memset SET数组为-1 表示当前节点无set标记


全程无槽点..1A了...   跑了600ms、留在了VJ第一页、算不错啦。。。上一个线段树代码跑了全VJ倒数第二大哭


#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
#include <queue>
#include <map> 
#include <vector>
using namespace std;
#define inf 2147483647
const int N = 1000050;
int  sum[21][2*N], add[21][2*N], SET[21][2*N];
int  maxx[21][2*N], minn[21][2*N];

int max(int a,int b)
{return a<b?b:a;}
int min(int a,int b)
{return a>b?b:a;}	
 
void pushDown(int i, int l, int r,int num,int op)		//把i节点的延迟标记传递到左右儿子节点
{
	if (op==2)
	{
		if (SET[num][i]!=-1)
		{
			int mid = (l + r) >> 1;
			SET[num][i << 1] = SET[num][i];
			sum[num][i << 1] = (mid - l + 1) * SET[num][i];  //[l, mid]代表左儿子区间
			SET[num][i << 1 | 1] = SET[num][i];
			sum[num][i << 1 | 1] = (r - mid) * SET[num][i];  //[mid + 1, r]代表右儿子区间
			minn[num][i<<1]=maxx[num][i<<1]=maxx[num][i<<1|1]=minn[num][i<<1|1]=SET[num][i];
			add[num][i<<1]=add[num][i<<1|1]=0;
				SET[num][i] = -1;
			
		}
	}
	else
		if (op==1)
		{
			if(add[num][i] != 0)
			{ 
				int mid = (l + r) >> 1;	
				pushDown(i<<1,l,mid,num,2);
				pushDown(i<<1|1,mid+1,r,num,2);//把儿子的SET传下去

				add[num][i << 1 | 1] += add[num][i];
				add[num][i << 1] += add[num][i];
				sum[num][i << 1] += (mid - l + 1) * add[num][i];  //[l, mid]代表左儿子区间
				sum[num][i << 1 | 1] += (r - mid) * add[num][i];  //[mid + 1, r]代表右儿子区间
				minn[num][i<<1]+=add[num][i];
				maxx[num][i<<1]+=add[num][i];
				maxx[num][i<<1|1]+=add[num][i];
				minn[num][i<<1|1]+=add[num][i];
				add[num][i] = 0;
			}
		}
}

void update(int i, int l, int r, int ql, int qr, int val,int num,int op) //更新区间为qlqr,当前区间为l,r,代表当前区间和的节点为i，更新值为val,op为1是add，2是SET
{
    if(l > qr || ql > r)	
		//更新区间不在当前区间内
        return ;
    if(l >= ql && r <= qr)	//要更新的区间把当前区间完全包括，则把当前整个区间+val，然后返回上一层
	{	
		if (op==1)  //add
		{
			if (SET[num][i]!=-1)
				pushDown(i,l,r,num,2);//add之前先把SET传下去
 
			sum[num][i] += (r - l + 1) * val;
			maxx[num][i]+=val;
			minn[num][i]+=val;
			add[num][i]+=val;
        }
		else	
			if (op==2)			//SET
			{
				add[num][i]=0;
				sum[num][i]= (r - l + 1) * val;
				maxx[num][i]=val;
				minn[num][i]=val;
				SET[num][i]=val;
			}
			return ;
    }
				//如果上面没reutrn 表示要往左右儿子区间update，所以把延迟标记放下去
	pushDown(i,l,r,num,2);	//先传SET
		pushDown(i, l, r,num,1);	//再传add
	
	
    int mid = (l + r) >> 1;
    update(i << 1, l, mid, ql, qr, val,num,op);
    update(i << 1 | 1, mid + 1, r, ql, qr, val,num,op);

    sum[num][i] = sum[num][i << 1] + sum[num][i << 1 | 1];
	maxx[num][i] = max(maxx[num][i << 1] , maxx[num][i << 1 | 1]);
	minn[num][i] = min(minn[num][i << 1] , minn[num][i << 1 | 1]);


}

int query(int i, int l, int r, int ql, int qr,int num)	 //查询区间为qlqr,当前区间为l,r,代表当前区间和的节点为i 
{
    if(l > qr || ql > r)
        return 0;
    if(l >= ql && r <= qr)
        return sum[num][i];

	pushDown(i,l,r,num,2);	//先传SET
		pushDown(i, l, r,num,1);	//再传add
    int mid =( l + r) >> 1;
    return query(i << 1, l, mid, ql, qr,num)   + query(i << 1 | 1, mid + 1, r, ql, qr,num);
}
int query_max(int i, int l, int r, int ql, int qr,int num)	 //查询区间为qlqr,当前区间为l,r,代表当前区间和的节点为i 
{
    if(l > qr || ql > r)
        return 0;
    if(l >= ql && r <= qr)
        return maxx[num][i];

	pushDown(i,l,r,num,2);	//先传SET
		pushDown(i, l, r,num,1);	//再传add
    int mid =( l + r) >> 1;
    return max(query_max(i << 1, l, mid, ql, qr,num)   , query_max(i << 1 | 1, mid + 1, r, ql, qr,num));
}
int query_min(int i, int l, int r, int ql, int qr,int num)	 //查询区间为qlqr,当前区间为l,r,代表当前区间和的节点为i 
{
    if(l > qr || ql > r)
        return inf;
    if(l >= ql && r <= qr)
        return minn[num][i];

	pushDown(i,l,r,num,2);	//先传SET
		pushDown(i, l, r,num,1);	//再传add
    int mid =(  l + r) >> 1;
 return min(query_min(i << 1, l, mid, ql, qr,num)   , query_min(i << 1 | 1, mid + 1, r, ql, qr,num));
}
int main()
{
	int op;
    int r,c,m,i,j;
	int x1,x2,y1,y2,v;
	while(    scanf("%d%d%d", &r,&c,&m)!=EOF)
	{
		memset(sum,0,sizeof(sum));
		memset(add,0,sizeof(add));
		memset(maxx,0,sizeof(maxx));
		memset(minn,0,sizeof(minn));
		memset(SET,-1,sizeof (SET));
		for(  i = 1; i <= m; i++)
		{
			scanf("%d",&op); 
			if (op==1)
			{
				scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&v);
				for (j=x1;j<=x2;j++)
				{
					update(1,1,c,y1,y2,v,j,1);
				}
			}
			else
				if (op==2)
				{	
					scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&v); 
					for (j=x1;j<=x2;j++)
					{
						update(1,1,c,y1,y2,v,j,2);
					}
					
				}
				else
				{
					scanf("%d%d%d%d",&x1,&y1,&x2,&y2); 
					int ans1=0;
					int ans2=0;
					int ans3=inf;
					for (j=x1;j<=x2;j++)
					{
						ans1+=query(1,1,c,y1,y2,j); 
						ans2=max(ans2,query_max(1,1,c,y1,y2,j));
						ans3=min(ans3,query_min(1,1,c,y1,y2,j)); 
					}
					printf("%d %d %d\n",ans1,ans3,ans2);
					
				}
				
		}
		
	}
	return 0;
}




yuhong_liu
关注
专栏目录
线段树成段更新----A Simple Problem with Integers
杨鑫newlife的专栏
06-11 987
Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 56841   Accepted: 17230 Case Time Limit: 2000MS Description You have N integers, A1, A2, ... , AN. You need to d
Mayor's posters --- 线段树成段更新，离散化
杨鑫newlife的专栏
06-20 973
Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 39339   Accepted: 11422 Description The citizens of Bytetown, AB, could not stand that the candidates in the mayoral e
E. Sasha and Array——矩阵+线段树
来自Fop_zz的辣鸡题解或教程
06-15 683
E.Sasha and Array 很神奇的一道题 题意大概就是支持区间加，然后求区间和。。。 区间和的求法是 sigema(i,l,r) f[i]  f[i]表示斐波那契序列第i项 例如一个区间的数字是1 2 3 4 5 6  那么其和为1+1+2+3+5+8=20 看到数据范围。。那斐波那契肯定不能离线处理了，只能是用矩阵去算 那我们如果给一个数+x，相当于让他再乘x次，
UVA11992:Fast Matrix Operations(线段树)
ACM！荣耀之路！
06-14 1674
题意： 有一个初始状态全为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 v：查询子矩阵(x1,y1,x2,y2)的元素和，最大值和最小值 思路： 因为总元素葛素不超过10^6，而且更新是对于连续的行进行更新，
Fast Matrix Operations
weixin_30861797的博客
07-10 99
题目：https://cn.vjudge.net/problem/UVA-11992 给定一个r*c(r<=20,r*c<=1e6)的矩阵，其元素都是0，现在对其子矩阵进行操作。 1 x1 y1 x2 y2 val 表示将(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素加val； 2x1 y1 x2 y2 val 表示将(x1,y1,x2,...
Fast Matrix Operations UVA - 11992
海边拾贝，沧海一粟
05-01 133
文章目录 原版刘汝佳的代码有很多maintain语句，个人感觉这样很不符合习惯，所以做出了修改，更符合我的习惯 1： 只保留一个update时的maintain语句 2：时刻记着延迟标记只是延迟往下的标记，并不延迟节点本身，时刻修改，当延迟标记下到某个地方，它的值就应该被修改 // UVa11992 Fast Matrix Operations（更易读、更具一般性的版本） // Rujia Liu...
Fast Matrix Calculation（矩阵快速幂）
Simone chou...
08-19 110
Fast Matrix Calculation Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)Total Submission(s): 28    Accepted Submission(s): 9 Problem Description One day, Alice...
Problem 1007 幸运数 线段树成段更新
05-15
现在有 $m$ 次操作，每次操作为将区间 $[l,r]$ 内的数加上 $x$（即 $a_l\gets a_l+x,a_{l+1}\gets a_{l+1}+x,\cdots,a_r\gets a_r+x$）。 对于每次操作，求出操作后 $[1,n]$ 中有多少个区间的权值为素数。 输入...
《软件方法2024版》公开内容202405更新-epub版
最新发布
05-22
目前指正人有（按指正时间排序）：吴佰钊、王周文、刘学斌、成文华、黄树成、李蜀斌、杨雪鸿、王书伟、高洪江、张志坚、龙燔、陈文飞、郭沼兵、陈自平、张彬、李宏伟、赵志军、孙赛刚、孙军、左科。 2018版《软件...
线段树区间更新
liuyanfeier的专栏
08-12 928
线段树成段更新延迟标记理解区间更新是指每次更新的时候更新的是一个区间里面的所有值，例如将区间[l,r]内的所有点都加或者减去一个数，或者替换成一个数字等等.因为区间更新每次更新的不止一个叶子节点，而叶子节点的值的更新肯定会影响到他的一系列的节点，所以假如每次都按照单点更新的思路将每一个叶子节点以及他的父节点都更新的话，那么工作量太大了并且时间肯定不会是log(n)会超时的。例如线段树总区间为[1,
UVA 11992 Fast Matrix Operations
码代码的猿猿
08-27 1102
每一行都建一个线段树。。。。 Fast Matrix Operations There is a matrix containing at most 106 elements divided into r rows and c columns. Each element has a location (x,y) where 1 1 x1 y1 x2 y2 v
UVA11992 Fast Matrix Operations 快速矩阵操作（线段树）
常胜将军的博客
07-09 189
题意：对矩阵有三种操作，第一是让矩阵的每个元素加上某个数，第二是让矩阵得每个元素变成某个数，第三是输出某个子矩阵的数字总和，最大值和最小值。开始时，矩阵得每个元素为0. 分析：每行建一颗线段树，多维护几个信息。详见《训练指南》P207-208。 代码： #include<bits/stdc++.h> using namespace std; const int maxnode...
UVa11992 - Fast Matrix Operations（线段树）
Coco_T的博客
10-23 257
题目链接简介： 给出一个矩阵，支持以下操作： 子矩阵加，子矩阵赋值，子矩阵查询元素和，最大值和最小值分析： 矩阵不超过20行 所以我们可以每一行都建一棵线段树，这样就转换成了一维问题set和add操作，我们默认set的优先级高 但是我认为，因为我们随时进行push 应该不会有操作叠加或者同时又add和set标记的情况tipUVa上一直RE。。。STO//这里写RE的代码片 #includ
UVA11992 Fast Matrix Operations（线段树）
酱油拌面
08-23 346
题目链接：点击链接题意：给一个矩阵，有三个操作。操作1，将这个矩阵的一个子矩阵的值全加上v。操作2，将这个矩阵的一个子矩阵的值全置为v。操作3，求这个矩阵的一个子矩阵的最大值，最小值，和。思路：看似很简单的set和add操作，我这个线段树菜鸡，卡了很久。前几天写另一道线段树时，突然明白set操作和add操作先后顺序对答案的影响。都明白先add后set，前面的add会被置0，执行Set操作的时候，懒标
uva 11992 - Fast Matrix Operations(线段树)
不慌不忙、不急不躁
08-25 1282
题目链接：uva 11992 - Fast Matrix Operations 题目大意：给定一个R∗C的矩阵，初始状态每个位置均为0， 然后进行Q次操作 1,x1,y1,x2,y2,v：将所有(x,y)满足（x1≤x≤x2,y1≤y≤y2）的点加上v2,x1,y1,x2,y2,v：将所有(x,y)满足（x1≤x≤x2,y1≤y≤y2）的点变成v3,x1,y1,x2,y2：求所
线段树（Fast Matrix Operations，UVA 11992）
xl2015190026的博客
12-12 243
我真的是写线段树写懵逼了。 开始出现成吨低级错误。 代码 #include #define m ((l+r)>>1) #define ls (now<<1) #define rs (ls|1) #define maxn 1000010 #define mem(a) memset(a,0,sizeof(a)) using namespace std; int R,C,M; int o
uva 11992 Fast Matrix Operations 线段树
冷月之殇
08-03 1795
代码：
Uva 11992 - Fast Matrix Operations题解(线段树区间更新+区间Set+Add,查询最大值,最小值,总和)
Twilightuse93的专栏
09-24 669
Problem F Fast Matrix Operations There is a matrix containing at most 106 elements divided into r rows and c columns. Each element has a location (x,y) where 1 1 x1 y1 x2 y2 v Increment each ele
UVA 11992(Fast Matrix Operations-线段树区间加&改)[Template:SegmentTree]
nike0good |Oier&ACMer | 熟能生巧
05-13 946
There is a matrix containing at most 106 elements divided into r rows and c columns. Each element has a location (x; y) where 1  x  r, 1  y  c. Initially, all the elements are zero. You need to
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值