[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.团队选题:餐厅到店点 ...
随机推荐
- bzoj4589 FWT xor版本
4589: Hard Nim Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 865 Solved: 484[Submit][Status][Disc ...
- PowerShell官方文档
PowerShell PowerShell 在 .NET Framework 基础之上构建,是一种基于任务的命令行 Shell 脚本语言:专门面向系统管理员和高级用户,可快速自动化多个操作系统(Lin ...
- Linux引导过程
早期时,启动一台计算机意味着要给计算机喂一条包含引导程序的纸带,或者手工使用前端面板地址/数据/控制开关来加载引导程序.尽管目前的计算机已经装备了很多工具来简化引导过程,但是这一切并没有对整个过程进行 ...
- 【BZOJ】5028: 小Z的加油店
[算法]数学+线段树/树状数组 [题解] 首先三个操作可以理解为更相减损术或者辗转相除法(待证明),所以就是求区间gcd. 这题的问题在线段树维护gcd只能支持修改成一个数,不支持加一个数. 套路:g ...
- set .net principle
var ticket = new FormsAuthenticationTicket(1, username, DateTime.Now, DateTime.Now.AddMinutes(FormsA ...
- 【洛谷 P1070】道路游戏 (DP)
题目链接 这题还是很好想的,看到\(90%\)的数据点时,我就知道要用\(n^3\)的算法(最后10分就算了吧) 然后,数据水,直接暴力\(n^3\)卡过了. 显然是道DP. 设\(f[i]\)表示第 ...
- [bzoj1031][JSOI2007]字符加密Cipher——后缀数组
Brief Description 给定一个长度为n的字符串,你需要对其进行加密. 把字符串围成一个环 显然从任意一个位置开始都可以有一个长度为n的串 把产生的n个串按字典序排序,把这n个串的最后一个 ...
- CentOS 7 主机加固手册-下
CentOS 7 主机加固手册-上 CentOS 7 主机加固手册-中 CentOS 7 主机加固手册-下 0x1f 删除禁用非必要的服务 删除非必要的服务 # Remove yum remove ...
- 内置函数补充,__str__方法、__del__方法 和 __call__方法和元祖
一 .内置函数补充 1.isinstance函数: isinstance(obj,cls)检查obj是否是类 cls 的对象 使用该函数来判断一个函数的类型 2. issubclass(sub, su ...
- JavaScript秘密花园
译文地址 bonsaiden.github.io/JavaScript-Garden/zh/#intro.authors 之前被人问到JS一些概念性的东西,感觉很模糊,可能层次比较浅,偏理论的东西实践 ...