火神最爱的就是吃鱼了,所以某一天他来到了一个池塘边捕鱼。池塘可以看成一个二维的平面,而他的渔网可以看成一个与坐标轴平行的矩形。
池塘里的鱼不停地在水中游动,可以看成一些点。有的时候会有鱼游进渔网,有的时候也会有鱼游出渔网。所以火神不知道什么时候收网才可以抓住最多的鱼,现在他寻求你的帮助。
他对池塘里的每条鱼都给予了一个标号,分别从$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]火神的鱼的更多相关文章

  1. AlloyTouch 0.2.0发布--鱼和熊掌兼得

    原文链接:https://github.com/AlloyTeam/AlloyTouch/wiki/AlloyTouch-0.2.0 背景 公司师姐昨日在KM发了篇长文,主要结论RAF+transfo ...

  2. 鱼搜_鱼搜官网_鱼搜搜索_http://www.7yusou.com

    收集了N多视频小站,然后花了3天时间弄了一个鱼搜搜索网站.欢迎大家访问哟. http://www.7yusou.com

  3. 洛谷 P1736 创意吃鱼法 Label:dp || 前缀和

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

  4. [课程设计]Scrum 3.8 多鱼点餐系统开发进度(留言反馈系统设计)

    Scrum 3.8 多鱼点餐系统开发进度(留言反馈系统设计) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点餐系统 ...

  5. [课程设计]Scrum 3.7 多鱼点餐系统开发进度(留言板选择方案)

    Scrum 3.7 多鱼点餐系统开发进度(留言板选择方案) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点餐系统W ...

  6. [课程设计]Scrum 3.6 多鱼点餐系统开发进度(用户测试反馈页面构思&留言板设计)

    Scrum 3.6 多鱼点餐系统开发进度(用户测试反馈页面构思&留言板设计) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团 ...

  7. [课程设计]Scrum 3.5 多鱼点餐系统开发进度(修复Bug&美化页面)

    Scrum 3.5 多鱼点餐系统开发进度(修复Bug&美化页面) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅 ...

  8. [课程设计]Scrum 3.4 多鱼点餐系统开发进度(下单详细信息页面&会员信息页面)

    Scrum 3.4 多鱼点餐系统开发进度(下单详细信息页面&会员信息页面) 1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队 ...

  9. [课程设计]Scrum 3.3 多鱼点餐系统开发进度(下单详细信息页面设计)

    Scrum 3.3 多鱼点餐系统开发进度(下单详细信息页面设计)  1.团队名称:重案组 2.团队目标:长期经营,积累客户充分准备,伺机而行 3.团队口号:矢志不渝,追求完美 4.团队选题:餐厅到店点 ...

随机推荐

  1. mootools框架里如何使用ajax

    ajax可通过直接写源码实现,但有点繁琐,现在流行的ajax框架都集成了ajax的功能,而且写起来非常简单方便.当然mootools也不例外.mootools是一个非常优秀的javascript的库, ...

  2. 一个JavaScript日期格式化扩展函数

    我们都知道在Java和PHP语言中,有专门用于格式化日期对象的类和函数,例如Java中的DateFormat等等,通过这些类和函数,我们可以方便的将一个日期对象按照格式的要求输出为字符串,例如对于同一 ...

  3. 关于 WizTools.org RESTClient的使用

    今天分享一个很好用的测试service的工具,很好用 提供两种方法使用这个东东. 第一种方法 通过cmd命令窗口. (1)cd C:\Users\li_weifeng\Desktop\c4d2(文件存 ...

  4. bzoj 2525 [Poi2011]Dynamite 二分+树形dp

    [Poi2011]Dynamite Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 270  Solved: 138[Submit][Status][D ...

  5. cloudera manager配置

    cloudera manager的数据库配置文件位置:    /etc/cloudera-scm-server/db.properties

  6. Linux设置虚拟内存-创建和启用Swap交换区

    如果你的服务器的总是报告内存不足,并且时常因为内存不足而引发服务被强制kill的话,在不增加物理内存的情况下,启用swap交换区作为虚拟内存是一个不错的选择,如果是SSD硬盘,正常读写速度都在300M ...

  7. iebackground+icon图标兼容

    <!DOCTYPE > <html> <head> <title>zepto</title> <meta name="nam ...

  8. MyBatis查询结果resultType返回值类型详细介绍

    一.返回一般数据类型 比如要根据 id 属性获得数据库中的某个字段值. mapper 接口: // 根据 id 获得数据库中的 username 字段的值 String getEmpNameById( ...

  9. bzoj4900 [CTSC2017]密钥

    传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4900 [题解] 恭喜bzoj达到40页 考场由于傻逼基数排序写挂了而gg. 竟然忘了考试前一 ...

  10. bzoj 2190 线性生成欧拉函数表

    首先我们知道,正方形内个是对称的,关于y=x对称,所以只需要算出来一半的人数 然后乘2+1就行了,+1是(1,1)这个点 开始我先想的递推 那么我们对于一半的三角形,一列一列的看,假设已经求好了第I- ...