博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
UVA-11992 Fast Matrix Operations
阅读量:7197 次
发布时间:2019-06-29

本文共 3439 字,大约阅读时间需要 11 分钟。

 

题意:矩阵加,矩阵覆盖,矩阵和查询

 

题目中说了不超过20行,所以不用二维线段树,20颗线段树即可搞定

#include
#include
#include
#define N 1000001using namespace std;int n,m,q,xx1,xx2,yy1,yy2,type;int opl,opr,w,tot;long long key[N<<2],tag_add[N<<2],tag_set[N<<2],tmin[N<<2],tmax[N<<2],sum,minn,maxn;int root[21],lc[N<<2],rc[N<<2];void clear(){ tot=0; memset(root,0,sizeof(root)); memset(lc,0,sizeof(lc)); memset(rc,0,sizeof(rc)); memset(key,0,sizeof(key)); memset(tmin,0,sizeof(tmin)); memset(tmax,0,sizeof(tmax)); memset(tag_add,0,sizeof(tag_add)); memset(tag_set,0,sizeof(tag_set));}void down_add(int k,int l,int r){ if(!lc[k]) lc[k]=++tot; if(!rc[k]) rc[k]=++tot; int mid=l+r>>1; if(tag_set[lc[k]]) tag_set[lc[k]]+=tag_add[k]; else tag_add[lc[k]]+=tag_add[k]; key[lc[k]]+=(mid-l+1)*tag_add[k]; tmin[lc[k]]+=tag_add[k]; tmax[lc[k]]+=tag_add[k]; if(tag_set[rc[k]]) tag_set[rc[k]]+=tag_add[k]; else tag_add[rc[k]]+=tag_add[k]; key[rc[k]]+=(r-mid)*tag_add[k]; tmin[rc[k]]+=tag_add[k]; tmax[rc[k]]+=tag_add[k]; tag_add[k]=0;}void down_set(int k,int l,int r){ if(!lc[k]) lc[k]=++tot; if(!rc[k]) rc[k]=++tot; int mid=l+r>>1; key[lc[k]]=(mid-l+1)*tag_set[k]; key[rc[k]]=(r-mid)*tag_set[k]; tmin[lc[k]]=tmax[lc[k]]=tmin[rc[k]]=tmax[rc[k]]=tag_set[k]; tag_add[lc[k]]=tag_add[rc[k]]=0; tag_set[lc[k]]=tag_set[rc[k]]=tag_set[k]; tag_set[k]=0;}void query(int k,int l,int r){ if(!k) return; if(l>=opl && r<=opr) { sum+=key[k]; minn=min(minn,tmin[k]); maxn=max(maxn,tmax[k]); return ; } if(tag_set[k]) down_set(k,l,r); if(tag_add[k]) down_add(k,l,r); int mid=l+r>>1; if(opl<=mid) query(lc[k],l,mid); if(opr>mid) query(rc[k],mid+1,r);}void add(int &k,int l,int r){ if(!k) k=++tot; if(l>=opl && r<=opr) { key[k]+=w*(r-l+1); tmin[k]+=w; tmax[k]+=w; tag_add[k]+=w; return; } if(tag_set[k]) down_set(k,l,r); if(tag_add[k]) down_add(k,l,r); int mid=l+r>>1; if(opl<=mid) add(lc[k],l,mid); if(opr>mid) add(rc[k],mid+1,r); key[k]=key[lc[k]]+key[rc[k]]; tmin[k]=min(tmin[lc[k]],tmin[rc[k]]); tmax[k]=max(tmax[lc[k]],tmax[rc[k]]);}void set(int &k,int l,int r){ if(!k) k=++tot; if(l>=opl && r<=opr) { key[k]=(r-l+1)*w; tag_add[k]=0; tag_set[k]=w; tmin[k]=tmax[k]=w; return; } if(tag_set[k]) down_set(k,l,r); if(tag_add[k]) down_add(k,l,r); int mid=l+r>>1; if(opl<=mid) set(lc[k],l,mid); if(opr>mid) set(rc[k],mid+1,r); key[k]=key[lc[k]]+key[rc[k]]; tmin[k]=min(tmin[lc[k]],tmin[rc[k]]); tmax[k]=max(tmax[lc[k]],tmax[rc[k]]);}int main(){ while(scanf("%d%d%d",&n,&m,&q)!=EOF) { for(int j=1;j<=q;j++) { scanf("%d%d%d%d%d",&type,&xx1,&yy1,&xx2,&yy2); if(type!=3) scanf("%d",&w); sum=0; minn=2e15; maxn=-2e15; opl=yy1; opr=yy2; for(int i=xx1;i<=xx2;i++) { if(type==1) add(root[i],1,m); else if(type==2) set(root[i],1,m); else query(root[i],1,m); } if(type==3) { if(minn==2e15) minn=0; if(maxn==-2e15) maxn=0; printf("%lld %lld %lld\n",sum,minn,maxn); } } clear(); }}

 

转载于:https://www.cnblogs.com/TheRoadToTheGold/p/7138338.html

你可能感兴趣的文章
matlab从fig文件中提取数据
查看>>
mysql show profiles 使用分析sql 性能
查看>>
POJ 2488 A Knight's Journey (DFS)
查看>>
jvm内存
查看>>
silverlight在XAML资源中实例化CLR对象
查看>>
Java中的字符串不变性
查看>>
Sql 的 RAISERROR用法
查看>>
Css下拉菜单设置
查看>>
Robot Framework学习笔记(八)------ride标签使用
查看>>
BIOS与EC之间关系
查看>>
一款简洁而强大的前端框架—JQuery
查看>>
js中的解绑事件
查看>>
ubuntu16.04下编译安装vim8.1
查看>>
DSSM 深度学习解决 NLP 问题:语义相似度计算
查看>>
真实世界的脉络].(英)戴维.多伊奇.pdf
查看>>
POJ 3710 Christmas Game
查看>>
秒懂神经网络---真·模拟退火算法
查看>>
js进阶 10-9 -of-type型子元素伪类选择器
查看>>
html5--6-14 CSS3中的颜色表示方式
查看>>
thinkphp5空控制器和空操作
查看>>