[Contest20171028]火神的鱼
火神最爱的就是吃鱼了,所以某一天他来到了一个池塘边捕鱼。池塘可以看成一个二维的平面,而他的渔网可以看成一个与坐标轴平行的矩形。
池塘里的鱼不停地在水中游动,可以看成一些点。有的时候会有鱼游进渔网,有的时候也会有鱼游出渔网。所以火神不知道什么时候收网才可以抓住最多的鱼,现在他寻求你的帮助。
他对池塘里的每条鱼都给予了一个标号,分别从$1$到$n$标号,$n$表示池塘里鱼的总数。鱼的游动可以概括为两个动作:
$1\ l\ r\ d$:表示标号在$[l,r]$这个区间内的鱼向$x$轴正方向游动了$d$个单位长度。
$2\ l\ r\ d$:表示标号在$[l,r]$这个区间内的鱼向$y$轴正方向游动了$d$个单位长度。
在某些时刻,火神会询问你现在有多少条鱼在渔网内(边界上的也算),请你来帮助他吧。
渔网的左下角坐标和右上角坐标为$(x_1,y_1)$和$(x_2,y_2)$
这大概是我近几个月来写的代码量最大的题了,A掉的时候感觉世界都明亮了起来23333
直接上线段树,每个节点维护$7$个标记:
$mx_1$表示当前节点的子树中最大的$x\lt x_1$
$mx_2$表示当前节点的子树中最大的$x_1\leq x\leq x_2$
$my_1$表示当前节点的子树中最大的$y\lt y_1$
$my_2$表示当前节点的子树中最大的$y_1\leq y\leq y_2$
$sum$表示当前节点所代表的区间中(入网的鱼的数量)
$laz_x$表示$x$轴增值的懒标记
$laz_y$表示$y$轴增值的懒标记
对于修改操作:
①直接区间修改的同时维护前四个标记
②当$mx_1\geq x_1$时,不停地把最大值移至$mx_2$并维护$sum$
③当$mx_2\gt x_2$时,不停地把最大值赋为$-\infty$并维护$sum$
因为鱼只会往正方向游,所以每条鱼最多转移标记一次,修改的总时间复杂度是$O\left(nlog_2n\right)$
对于查询:
因为修改的时候已经维护好了$sum$,所以我们查询的时候直接求$sum$的和即可
#include<stdio.h> #define inf 2147483647 int mx1[120010],mx2[120010],my1[120010],my2[120010],sum[120010],lazx[120010],lazy[120010],x[30010],y[30010],x1,y1,x2,y2; int max(int a,int b){return a>b?a:b;} void pushup(int x){ mx1[x]=max(mx1[x<<1],mx1[x<<1|1]); mx2[x]=max(mx2[x<<1],mx2[x<<1|1]); my1[x]=max(my1[x<<1],my1[x<<1|1]); my2[x]=max(my2[x<<1],my2[x<<1|1]); sum[x]=sum[x<<1]+sum[x<<1|1]; } void build(int l,int r,int u){ lazx[u]=lazy[u]=0; if(l==r){ if(x[l]<x1){ mx1[u]=x[l]; mx2[u]=-inf; sum[u]=0; }else if(x[l]<=x2){ mx1[u]=-inf; mx2[u]=x[l]; sum[u]=1; }else{ mx1[u]=mx2[u]=-inf; sum[u]=0; } if(y[l]<y1){ my1[u]=y[l]; my2[u]=-inf; sum[u]=0; }else if(y[l]<=y2){ my1[u]=-inf; my2[u]=y[l]; sum[u]&=1; }else{ my1[u]=my2[u]=-inf; sum[u]=0; } return; } int mid=(l+r)>>1; build(l,mid,u<<1); build(mid+1,r,u<<1|1); pushup(u); } void pushdown(int x){ if(lazx[x]){ lazx[x<<1]+=lazx[x]; mx1[x<<1]+=lazx[x]; mx2[x<<1]+=lazx[x]; lazx[x<<1|1]+=lazx[x]; mx1[x<<1|1]+=lazx[x]; mx2[x<<1|1]+=lazx[x]; lazx[x]=0; } if(lazy[x]){ lazy[x<<1]+=lazy[x]; my1[x<<1]+=lazy[x]; my2[x<<1]+=lazy[x]; lazy[x<<1|1]+=lazy[x]; my1[x<<1|1]+=lazy[x]; my2[x<<1|1]+=lazy[x]; lazy[x]=0; } } void modifyx(int L,int R,int v,int l,int r,int x){ if(L<=l&&r<=R){ mx1[x]+=v; mx2[x]+=v; lazx[x]+=v; return; } int mid=(l+r)>>1; pushdown(x); if(L<=mid)modifyx(L,R,v,l,mid,x<<1); if(mid<R)modifyx(L,R,v,mid+1,r,x<<1|1); pushup(x); } int querymx1(int L,int R,int l,int r,int x){ if(L<=l&&r<=R)return mx1[x]; pushdown(x); int mid=(l+r)>>1,ans=-inf; if(L<=mid)ans=max(ans,querymx1(L,R,l,mid,x<<1)); if(mid<R)ans=max(ans,querymx1(L,R,mid+1,r,x<<1|1)); return ans; } int querymx2(int L,int R,int l,int r,int x){ if(L<=l&&r<=R)return mx2[x]; pushdown(x); int mid=(l+r)>>1,ans=-inf; if(L<=mid)ans=max(ans,querymx2(L,R,l,mid,x<<1)); if(mid<R)ans=max(ans,querymx2(L,R,mid+1,r,x<<1|1)); return ans; } void modifyx12(int L,int R,int l,int r,int x){ if(l==r){ if(mx1[x]>=x1){ mx2[x]=mx1[x]; mx1[x]=-inf; } if(mx2[x]>=x1&&my2[x]>=y1) sum[x]=1; else sum[x]=0; return; } pushdown(x); int mid=(l+r)>>1; if(mx1[x<<1]>mx1[x<<1|1]) modifyx12(L,R,l,mid,x<<1); else modifyx12(L,R,mid+1,r,x<<1|1); pushup(x); } void modifyx20(int L,int R,int l,int r,int x){ if(l==r){ if(mx2[x]>x2)mx1[x]=mx2[x]=-inf; if(mx2[x]>=x1&&my2[x]>=y1) sum[x]=1; else sum[x]=0; return; } pushdown(x); int mid=(l+r)>>1; if(mx2[x<<1]>mx2[x<<1|1]) modifyx20(L,R,l,mid,x<<1); else modifyx20(L,R,mid+1,r,x<<1|1); pushup(x); } void modifyy(int L,int R,int v,int l,int r,int x){ if(L<=l&&r<=R){ my1[x]+=v; my2[x]+=v; lazy[x]+=v; return; } int mid=(l+r)>>1; pushdown(x); if(L<=mid)modifyy(L,R,v,l,mid,x<<1); if(mid<R)modifyy(L,R,v,mid+1,r,x<<1|1); pushup(x); } int querymy1(int L,int R,int l,int r,int x){ if(L<=l&&r<=R)return my1[x]; pushdown(x); int mid=(l+r)>>1,ans=-inf; if(L<=mid)ans=max(ans,querymy1(L,R,l,mid,x<<1)); if(mid<R)ans=max(ans,querymy1(L,R,mid+1,r,x<<1|1)); return ans; } int querymy2(int L,int R,int l,int r,int x){ if(L<=l&&r<=R)return my2[x]; pushdown(x); int mid=(l+r)>>1,ans=-inf; if(L<=mid)ans=max(ans,querymy2(L,R,l,mid,x<<1)); if(mid<R)ans=max(ans,querymy2(L,R,mid+1,r,x<<1|1)); return ans; } void modifyy12(int L,int R,int l,int r,int x){ if(l==r){ if(my1[x]>=y1){ my2[x]=my1[x]; my1[x]=-inf; } if(mx2[x]>=x1&&my2[x]>=y1) sum[x]=1; else sum[x]=0; return; } pushdown(x); int mid=(l+r)>>1; if(my1[x<<1]>my1[x<<1|1]) modifyy12(L,R,l,mid,x<<1); else modifyy12(L,R,mid+1,r,x<<1|1); pushup(x); } void modifyy20(int L,int R,int l,int r,int x){ if(l==r){ if(my2[x]>y2)my1[x]=my2[x]=-inf; if(mx2[x]>=x1&&my2[x]>=y1) sum[x]=1; else sum[x]=0; return; } pushdown(x); int mid=(l+r)>>1; if(my2[x<<1]>my2[x<<1|1]) modifyy20(L,R,l,mid,x<<1); else modifyy20(L,R,mid+1,r,x<<1|1); pushup(x); } int query(int L,int R,int l,int r,int x){ if(L<=l&&r<=R)return sum[x]; int mid=(l+r)>>1,ans=0; if(L<=mid)ans+=query(L,R,l,mid,x<<1); if(mid<R)ans+=query(L,R,mid+1,r,x<<1|1); return ans; } int main(){ int T,n,m,i,op,l,r,d; scanf("%d",&T); while(T--){ scanf("%d%d%d%d%d",&n,&x1,&y1,&x2,&y2); for(i=1;i<=n;i++)scanf("%d%d",x+i,y+i); build(1,n,1); scanf("%d",&m); while(m--){ scanf("%d%d%d",&op,&l,&r); if(op==1){ scanf("%d",&d); modifyx(l,r,d,1,n,1); while(querymx1(l,r,1,n,1)>=x1)modifyx12(l,r,1,n,1); while(querymx2(l,r,1,n,1)>x2)modifyx20(l,r,1,n,1); } if(op==2){ scanf("%d",&d); modifyy(l,r,d,1,n,1); while(querymy1(l,r,1,n,1)>=y1)modifyy12(l,r,1,n,1); while(querymy2(l,r,1,n,1)>y2)modifyy20(l,r,1,n,1); } if(op==3)printf("%d\n",query(l,r,1,n,1)); } } }
[Contest20171028]火神的鱼的更多相关文章
- AlloyTouch 0.2.0发布--鱼和熊掌兼得
原文链接:https://github.com/AlloyTeam/AlloyTouch/wiki/AlloyTouch-0.2.0 背景 公司师姐昨日在KM发了篇长文,主要结论RAF+transfo ...
- 鱼搜_鱼搜官网_鱼搜搜索_http://www.7yusou.com
收集了N多视频小站,然后花了3天时间弄了一个鱼搜搜索网站.欢迎大家访问哟. http://www.7yusou.com
- 洛谷 P1736 创意吃鱼法 Label:dp || 前缀和
题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...
- [课程设计]Scrum 3.8 多鱼点餐系统开发进度(留言反馈系统设计)
Scrum 3.8 多鱼点餐系统开发进度(留言反馈系统设计) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点餐系统 ...
- [课程设计]Scrum 3.7 多鱼点餐系统开发进度(留言板选择方案)
Scrum 3.7 多鱼点餐系统开发进度(留言板选择方案) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点餐系统W ...
- [课程设计]Scrum 3.6 多鱼点餐系统开发进度(用户测试反馈页面构思&留言板设计)
Scrum 3.6 多鱼点餐系统开发进度(用户测试反馈页面构思&留言板设计) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团 ...
- [课程设计]Scrum 3.5 多鱼点餐系统开发进度(修复Bug&美化页面)
Scrum 3.5 多鱼点餐系统开发进度(修复Bug&美化页面) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅 ...
- [课程设计]Scrum 3.4 多鱼点餐系统开发进度(下单详细信息页面&会员信息页面)
Scrum 3.4 多鱼点餐系统开发进度(下单详细信息页面&会员信息页面) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队 ...
- [课程设计]Scrum 3.3 多鱼点餐系统开发进度(下单详细信息页面设计)
Scrum 3.3 多鱼点餐系统开发进度(下单详细信息页面设计) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点 ...
随机推荐
- mootools框架里如何使用ajax
ajax可通过直接写源码实现,但有点繁琐,现在流行的ajax框架都集成了ajax的功能,而且写起来非常简单方便.当然mootools也不例外.mootools是一个非常优秀的javascript的库, ...
- 一个JavaScript日期格式化扩展函数
我们都知道在Java和PHP语言中,有专门用于格式化日期对象的类和函数,例如Java中的DateFormat等等,通过这些类和函数,我们可以方便的将一个日期对象按照格式的要求输出为字符串,例如对于同一 ...
- 关于 WizTools.org RESTClient的使用
今天分享一个很好用的测试service的工具,很好用 提供两种方法使用这个东东. 第一种方法 通过cmd命令窗口. (1)cd C:\Users\li_weifeng\Desktop\c4d2(文件存 ...
- bzoj 2525 [Poi2011]Dynamite 二分+树形dp
[Poi2011]Dynamite Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 270 Solved: 138[Submit][Status][D ...
- cloudera manager配置
cloudera manager的数据库配置文件位置: /etc/cloudera-scm-server/db.properties
- Linux设置虚拟内存-创建和启用Swap交换区
如果你的服务器的总是报告内存不足,并且时常因为内存不足而引发服务被强制kill的话,在不增加物理内存的情况下,启用swap交换区作为虚拟内存是一个不错的选择,如果是SSD硬盘,正常读写速度都在300M ...
- iebackground+icon图标兼容
<!DOCTYPE > <html> <head> <title>zepto</title> <meta name="nam ...
- MyBatis查询结果resultType返回值类型详细介绍
一.返回一般数据类型 比如要根据 id 属性获得数据库中的某个字段值. mapper 接口: // 根据 id 获得数据库中的 username 字段的值 String getEmpNameById( ...
- bzoj4900 [CTSC2017]密钥
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4900 [题解] 恭喜bzoj达到40页 考场由于傻逼基数排序写挂了而gg. 竟然忘了考试前一 ...
- bzoj 2190 线性生成欧拉函数表
首先我们知道,正方形内个是对称的,关于y=x对称,所以只需要算出来一半的人数 然后乘2+1就行了,+1是(1,1)这个点 开始我先想的递推 那么我们对于一半的三角形,一列一列的看,假设已经求好了第I- ...