[APIO2018] New Home 新家 题解

于 2020-11-24 13:23:36 发布
阅读量209 收藏
点赞数
分类专栏： 题解 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42925924/article/details/110069998
版权
题解 同时被 2 个专栏收录
119 篇文章 1 订阅
订阅专栏
算法
89 篇文章 1 订阅
订阅专栏

[APIO2018] New Home 新家 题解

传送门

首先将时间离散化，然后依据时间建一棵线段树。将每一个店开业歇业时间搞到树上。

这样问题就转化为：有一个数轴每次加入或删除一个点。每一个点有一个颜色。在任意时间询问某一个点距离某一种颜色的最远距离。到某一个颜色的距离为到这个颜色的点的最小距离。

显然可以二分答案。每次询问区间 [ l , r ] [l,r] [l,r]是否包含了所有颜色的点。

对于每一个颜色我们搞一个set。记录每一个点到前面的那个同色的点，也就是前驱，记为pre[i]。

这个玩意的充要条件是在区间 [ l , ∞ ] [l,\infin] [l,∞] 的pre[i]都满足>=r。

对于这个我们也可以将坐标离散化搞到一个树上。

时间复杂度是 O ( N × log ⁡ 2 N ) O(N\times \log ^2 N) O(N×log2N)

这题也有单log的做法,就是在线段树上二分的地方的优化，不想多讲。

code:

(这个代码被卡常了我也懒得改了)

#include<bits/stdc++.h>
#define rb(a,b,c) for(int a=b;a<=c;++a)
#define rl(a,b,c) for(int a=b;a>=c;--a)
#define LL long long
#define IT iterator
#define PB push_back
#define II(a,b) make_pair(a,b)
#define FIR first
#define SEC second
#define FREO freopen("check.out","w",stdout)
#define rep(a,b) for(int a=0;a<b;++a)
#define SRAND mt19937 rng(chrono::steady_clock::now().time_since_epoch().count())
#define random(a) rng()%a
#define ALL(a) a.begin(),a.end()
#define POB pop_back
#define ff fflush(stdout)
#define fastio ios::sync_with_stdio(false)
#define check_min(a,b) a=min(a,b)
#define check_max(a,b) a=max(a,b)
using namespace std;
inline int read()
{
	int num = 0, f = 1;
	char ch = getchar();
	while( !isdigit( ch ) ) { if(ch == '-') f = -1; ch = getchar(); }
	while( isdigit( ch ) ) num = (num << 3) + (num << 1) + (ch ^ 48), ch = getchar();
	return num * f;
}
const int INF=0x3f3f3f3f;
typedef pair<int,int> mp;
/*}
*/ 
const int MAXN=6e5+20;
int n,k,q;
vector<mp> sor; 
vector<int> pos_;
int x[MAXN],t[MAXN],a[MAXN],b[MAXN],l[MAXN],y[MAXN];
map<int,int> M;
int find(int f){
	//第一个>=f的位置
	return lower_bound(ALL(pos_),f)-pos_.begin()+1; 
}
const int N=1<<20;
int rest[MAXN];

struct SEGMENTTREE{
	int tree[N+N];//最小 
	SEGMENTTREE(){
		memset(tree,63,sizeof(tree));
	} 
	void modify(int index,int val){
		index+=N-1;
		tree[index]=val;
		index>>=1;
		while(index){
			tree[index]=min(tree[index<<1],tree[index<<1|1]);
			index>>=1;
		}
	}
	int query(int a,int b,int now=1,int l=1,int r=N+1){
		if(r<=a||l>=b){
			return INF;
		}
		if(r<=b&&l>=a){
			return tree[now];
		}
		int mid=(l+r)>>1;
		return min(query(a,b,now<<1,l,mid),query(a,b,now<<1|1,mid,r));
	}
}sgt;
set<int> pos[MAXN];
void add(int x_,int t_){
	pos[t_].insert(x_);
	auto ite=pos[t_].upper_bound(x_);
	if(ite!=pos[t_].end()){
		sgt.modify(*ite,x_);
	}
	ite=pos[t_].lower_bound(x_);
	if(ite!=pos[t_].begin()){
		ite--;
		sgt.modify(x_,*ite);
	}
	else {
		sgt.modify(x_,-INF);
	}
}
void del(int x_,int t_){
	sgt.modify(x_,INF);
	pos[t_].erase(x_);
	auto ite=pos[t_].upper_bound(x_);
	if(ite==pos[t_].end()) return;
	if(ite==pos[t_].begin()){
		sgt.modify(*ite,-INF);
	}
	else{
		int z=*ite;
		ite--;
		sgt.modify(z,*ite);
	}
}
bool check(int x_,int dis){
	int st=find(x_+dis+1);
	if(sgt.query(st,N+1)<find(x_-dis)) return 0;
	return 1;
}
struct EVENTSEGTREE{
	vector<int> tree[N+N];
	EVENTSEGTREE(){
		
	}
	void add_event(int st,int ed,int val,int now=1,int l=1,int r=N+1){
		if(r<=st||l>=ed){
			return ;
		}
		if(r<=ed&&l>=st){
			tree[now].PB(val);
			return ;
		}
		int mid=(l+r)>>1;
		add_event(st,ed,val,now<<1,l,mid);
		add_event(st,ed,val,now<<1|1,mid,r);
	}
	
	void run(int now=1,int l_=1,int r_=N+1){
		for(auto it:tree[now]){
			if(it>0){
				add(x[it],t[it]);
			}	
		}
		if(l_==r_-1){
			for(auto it:tree[now]){
				if(it<0){
					int lb=0,rb=1e8+1;
					if(!check(l[-it],rb)){
						rest[-it]=-1;
					}
					else{
						while(lb<rb){
							int mid=(lb+rb)>>1;
							if(check(l[-it],mid)) rb=mid;
							else lb=mid+1;
						}
						rest[-it]=lb;
					}
				}
			}
		}
		else{
			int mid=(l_+r_)>>1;
			run(now<<1,l_,mid);
			run(now<<1|1,mid,r_);
		}
		for(auto it:tree[now]){
			if(it>0){
				del(x[it],t[it]);
			}
		}
	}
}esgt;
int main(){
//	scanf("%d%d%d",&n,&k,&q);
	n=read();
	k=read();
	q=read();
	rb(i,1,n){
		int xi,ti,ai,bi;
//		scanf("%d%d%d%d",&xi,&ti,&ai,&bi);
		xi=read();
		ti=read();
		ai=read();
		bi=read();
		a[i]=ai;
		b[i]=bi;
		M[ai]=M[bi]=1;
		t[i]=ti;
		sor.PB(II(xi,i));
	}	
	rb(i,1,q){
		l[i]=read();
		y[i]=read();
	}
	rb(i,1,q){
//		scanf("%d%d",&l[i],&y[i]);
		M[y[i]]=1;
	}
	rb(i,1,k){
		n++;
		a[n]=-INF;
		b[n]=INF;
		M[a[n]]=M[b[n]]=1;
		x[n]=INF;
		t[n]=i;
		sor.PB(II(x[n],n));
	}
	sort(ALL(sor));
	rep(i,n){
		x[sor[i].SEC]=i+1;
		pos_.PB(sor[i].FIR);
	}
	int cnt=0;
	for(auto ite=M.begin();ite!=M.end();ite++){
		ite->SEC=++cnt;
	}
	rb(i,1,n){
		a[i]=M[a[i]];
		b[i]=M[b[i]];
		esgt.add_event(a[i],b[i]+1,i);
	}
	rb(i,1,q){
		y[i]=M[y[i]];
		esgt.add_event(y[i],y[i]+1,-i);
	}
	esgt.run();
	rb(i,1,q){
		printf("%d\n",rest[i]);
	}
	return 0;
}

UPDATE:

还是补充一下 O ( N × log ⁡ N ) O(N\times \log N) O(N×logN)的做法吧，(UOJ 已ac！)

#include<bits/stdc++.h>
#define rb(a,b,c) for(int a=b;a<=c;++a)
#define rl(a,b,c) for(int a=b;a>=c;--a)
#define LL long long
#define IT iterator
#define PB push_back
#define II(a,b) make_pair(a,b)
#define FIR first
#define SEC second
#define FREO freopen("check.out","w",stdout)
#define rep(a,b) for(int a=0;a<b;++a)
#define SRAND mt19937 rng(chrono::steady_clock::now().time_since_epoch().count())
#define random(a) rng()%a
#define ALL(a) a.begin(),a.end()
#define POB pop_back
#define ff fflush(stdout)
#define fastio ios::sync_with_stdio(false)
#define check_min(a,b) a=min(a,b)
#define check_max(a,b) a=max(a,b)
using namespace std;
inline int read()
{
	int num = 0, f = 1;
	char ch = getchar();
	while( !isdigit( ch ) ) { if(ch == '-') f = -1; ch = getchar(); }
	while( isdigit( ch ) ) num = (num << 3) + (num << 1) + (ch ^ 48), ch = getchar();
	return num * f;
}
const int INF=0x3f3f3f3f;
typedef pair<int,int> mp;
/*}
*/ 
const int MAXN=6e5+20;
int n,k,q;
vector<mp> sor; 
vector<int> pos_;
int x[MAXN],t[MAXN],a[MAXN],b[MAXN],l[MAXN],y[MAXN];
map<int,int> M;
int find(int f){
	//第一个>=f的位置
	return lower_bound(ALL(pos_),f)-pos_.begin()+1; 
}
const int N=1<<20;
int rest[MAXN];

struct SEGMENTTREE{
	int tree[N+N];//最小 
	SEGMENTTREE(){
		memset(tree,63,sizeof(tree));
	} 
	void modify(int index,int val){
		index+=N-1;
		tree[index]=val;
		index>>=1;
		while(index){
			tree[index]=min(tree[index<<1],tree[index<<1|1]);
			index>>=1;
		}
	}
	int query(int a,int b,int now=1,int l=1,int r=N+1){
		if(r<=a||l>=b){
			return INF;
		}
		if(r<=b&&l>=a){
			return tree[now];
		}
		int mid=(l+r)>>1;
		return min(query(a,b,now<<1,l,mid),query(a,b,now<<1|1,mid,r));
	}
	int low_bound(int pos,int now=1,int l=1,int r=N+1,int MIN=INF){
		if(l==r-1) {
			int ret=INF;
			ret=min(ret,pos_[l-1]-pos);
			if(min(tree[now],MIN)>0)
			ret=min(ret,pos-pos_[min(tree[now],MIN)-1]);
			return ret;	
		}
		int mid=(l+r)>>1;
		if(mid>pos_.size()) return low_bound(pos,now<<1,l,mid,min(MIN,tree[now<<1|1]));
		if(pos_[mid-1]<=pos) return low_bound(pos,now<<1|1,mid,r,MIN);
		if(min(tree[now<<1|1],MIN)!=-INF&&(pos_[mid-2]-pos>=pos-pos_[min(tree[now<<1|1],MIN)-1])){
			return low_bound(pos,now<<1,l,mid,min(MIN,tree[now<<1|1]));
		}
		return low_bound(pos,now<<1|1,mid,r,MIN);
	}
}sgt;
set<int> pos[MAXN];
void add(int x_,int t_){
	pos[t_].insert(x_);
	auto ite=pos[t_].upper_bound(x_);
	if(ite!=pos[t_].end()){
		sgt.modify(*ite,x_);
	}
	ite=pos[t_].lower_bound(x_);
	if(ite!=pos[t_].begin()){
		ite--;
		sgt.modify(x_,*ite);
	}
	else {
		sgt.modify(x_,-INF);
	}
}
void del(int x_,int t_){
	sgt.modify(x_,INF);
	pos[t_].erase(x_);
	auto ite=pos[t_].upper_bound(x_);
	if(ite==pos[t_].end()) return;
	if(ite==pos[t_].begin()){
		sgt.modify(*ite,-INF);
	}
	else{
		int z=*ite;
		ite--;
		sgt.modify(z,*ite);
	}
}
bool check(int x_,int dis){
	int st=find(x_+dis+1);
	if(sgt.query(st,N+1)<find(x_-dis)) return 0;
	return 1;
}
int main(){
//	scanf("%d%d%d",&n,&k,&q);
	n=read();
	k=read();
	q=read();
	rb(i,1,n){
		int xi,ti,ai,bi;
//		scanf("%d%d%d%d",&xi,&ti,&ai,&bi);
		xi=read();
		ti=read();
		ai=read();
		bi=read();
		a[i]=ai;
		b[i]=bi;
		M[ai]=M[bi]=1;
		t[i]=ti;
		sor.PB(II(xi,i));
	}	
	rb(i,1,q){
		l[i]=read();
		y[i]=read();
	}
	rb(i,1,q){
//		scanf("%d%d",&l[i],&y[i]);
		M[y[i]]=1;
	}
	rb(i,1,k){
		n++;
		a[n]=-INF;
		b[n]=INF;
		M[a[n]]=M[b[n]]=1;
		x[n]=INF;
		t[n]=i;
		sor.PB(II(x[n],n));
	}
	sort(ALL(sor));
	rep(i,n){
		x[sor[i].SEC]=i+1;
		pos_.PB(sor[i].FIR);
	}
	int cnt=0;
	for(auto ite=M.begin();ite!=M.end();ite++){
		ite->SEC=++cnt;
	}
	vector<mp> eve;
	rb(i,1,n){
		a[i]=M[a[i]];
		b[i]=M[b[i]];
		eve.PB(II(a[i],i-INF));
		eve.PB(II(b[i]+1,i));
	}
	rb(i,1,q){
		y[i]=M[y[i]];
		eve.PB(II(y[i],i+INF));
	}
	sort(ALL(eve));
	for(auto &it:eve){
		if(it.SEC<0){
			add(x[it.SEC+INF],t[it.SEC+INF]);
			continue;
		}
		if(it.SEC<INF){
			del(x[it.SEC],t[it.SEC]);
			continue;
		}
		it.SEC-=INF;
		int lb=0,rb=1e8+1;
		if(!check(l[it.SEC],rb)){
			rest[it.SEC]=-1;
		}
		else{
				rest[it.SEC]=sgt.low_bound(l[it.SEC]);
			}
	}
	rb(i,1,q){
		printf("%d\n",rest[i]);
	}
	return 0;
}
Gary_2005
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
APIO2018T1 New Home 题解 吕欣
07-26
题目大意 在一个数轴上： 给定 N 个商店，每个商店有一个开业时间、关门时间、坐标和销售物品的种类 同时有 M 个询问，每个询问给你一个时间 t[i] 和地点 d[i]，试求在 t[i] 时刻，一个住在 d[i] 的人，为了买某种商品最远需要跑多长距离。
APIO2018解题报告
weixin_34238642的博客
12-19 173
今年的APIO好邪啊。 T1铁人两项 题目大意 给一个无向图，问有多少三元组(三个元素两两不同)使得它们构成一条简单路径 。 题解 无向图这种东西不好直接处理，考虑点双缩点建圆方树。 然后就出现了一个比较神的做法，把方点点权设为点双大小，圆点点权设为-1，那么我们枚举三元组中的起点和终点，那么这条路径的权值就是可能的中间点个数。 为什么？考虑两种圆点。 1、起点和终点，因为它们已经...
APIO2018 题解
weixin_44666356的博客
08-12 23
铁人两项 铁人两项 看到不经过重复点，想到网络流 / 圆方树。 前者复杂度太大不考虑，后者，缩点后成圆方树，似乎比较好统计？ 简单个 p 我写了一页的理论，结果实现调了一年都没有调出来，我裂开了。 对于满足条件的 \(s, c, f\)，一定是 \(c\) 在 \(s \Rightarrow f\) 的一条简单路径上，考虑 \(s\) 与 \(f\) 的所有简单路径的点集之并是 \(s\) 与 ...
【LOJ2585】「APIO2018」新家
cz_xuyixuan的博客
05-23 734
【题目链接】点击打开链接【思路要点】对时间轴进行扫描线，在时间轴上的一段区间可以表示为一次插入操作和一次删除操作。问题被转化为：维护一个序列，支持在某处插入/删除一个数，以及询问以某个位置为中心，包含所有种类的数的区间的最小长度。对于询问，不难想到二分答案，二分答案后我们需要支持的是询问区间内是否出现了所有种类的数。对于每一种数，我们维护一个支持查询前驱后继的平衡树（std::set）。在序列上维...
APIO2018 新家
yzyyylx的博客
12-29 498
题面 题意 在一条路上有n家商店，每个商店有一个类型，在某段时间内在某个位置上存在，每次询问给出一个时间和一个地点，求该时间从该地点出发最少走多少路才可以到达任意一种类型的商店。 做法 首先可以离线处理，根据时间进行排序，那么没家商店就都可以看作是在某个时间加入，再从某个时间删除，并询问某个点到所有类型点的最小距离的最大值。 这个可以进行二分，这样问题就转化为了，询问此时某一段区间内是否存在所有类...
APIO2018 T2T3 题解 陈俊锟 吕欣 张瑞喆
07-26
第二题 重复下列算法： 每次选出半径最大的圆 删除和它有交的所有圆 这里有交指圆盘公共区域不为 0 求每个圆是哪一次被删除的 第三题 先从简单的题开始 给出一张无向图 求有多少 (
APIO2018资料
06-30
APIO2018的题面与部分课件与pdf和ppt，.在此分享一下。
APIO2018游记
include / Hdhd
05-14 290
Day 0:碗了一套APIO2015的题。由于OJ没有打包测试的能力，加上试题暴力分分高，居然在没有AC的情况下上了200。信心满满。 Day 1：到北京已经是中午。午饭异常拥挤，开始怀念南开食堂。下午去清华膜拜nodgd和Wo_Ai_Wang_Yuan，共进晚餐。 Day 2：开始讲课了。第一位老师讲的是计算机与哲学。啥都没听懂，只是知道了有一系列的折纸问题。第二位老师较为有趣，讲了图像处理...
CTSC/APIO2018咕咕记&&THUPC2018打铁记
时光真疯狂, 我一路执迷于匆忙.
05-15 1280
CTSC Day -1 为了赶明早的飞机在广州住了一晚上 水果竟然是免费的好评 本能寺真好玩，我塔喵肝爆 Day 0 报道日，发现房间是按NOIP成绩分的。。。 于是我有幸和beginend dalao一间房，一起战斗很资瓷啊 八十中饭堂好评，就是为什么一张票只能换一种东西，或者6瓶优酸乳_ (:з」∠) _ 房间海星，至少是两张床的，而且麦老大和beginend的儿子同学也过来...
APIO2018小记
szh_0808的博客
05-12 473
一个不敢参加CTSC只来了APIO却依然没有什么好下场的蒟蒻的小记。 T1 New Home 一开始敲了Subtask1，n²暴力。 觉得只拿五分不甘心啊，去看Subtask2，想了一个神奇的做法：离线处理，将询问按时间排序，再将所有商店也按时间排序，枚举k，对于每一个k分别开一个set和一个priority_queue，set里面放坐标，priority_queue是pair&lt;结束时...
(APIO2018)从DFA到后缀自动机_张云帆.pdf
07-29
(APIO2018)从DFA到后缀自动机_张云帆.pdf
APIO2008.pdf
07-21
"APIO2008.pdf" 本文档是亚洲信息学奥林匹克竞赛（APIO）2008年的竞赛题目，包含三个问题，分别是珠链交换器、换器交和DNA。本文档详细介绍了每个问题的描述、输入格式、输出格式、时间限制和空间限制等信息。 ...
[笔记] APIO 2018 Day1
weixin_34075551的博客
05-11 164
计算折纸 computaional origami 全息算法(???) margulis napkin problem 素数里有任意长的等差数列 xor gate Σxi or gate(exact version) 1-π(1-xi) or gate(approximate versiton) :Randomly pick r1,r2,…,rn...
[APIO2018] New Home 新家
weixin_34111790的博客
02-23 123
扫描线+线段树+二分答案+set+STL 就是把区间数颜色做得很好 时间看成线段，扫描线 对于某一个询问位置x 二分答案转化，看区间内有没有k种颜色。。 一个区间数颜色的套路是，prei上一个该颜色出现位置 查[x-mid,x+mid]pre小于x-mid的有几个。 树套树！！（TLE飞起） 其实并不关心是哪些个，只关心是否有k种 对于(prei,i)开区间，一定没有i这个...
APIO2018凉凉记
Luckfort
05-14 301
总起 这次能够去APIO其实已经很不容易了，毕竟考虑刚考完GDOI文化课后我都成文盲了。 Day0 由于教练漏帮我报吃特色菜，所以我只能够来APIO了。一个人坐飞机，这是第一次。 辗转了很久终于来到了丽都，住的是广州二中大佬住过的房间，酒店的服务还挺周到的。 Day1 又久经辗转跑到海淀区？领完东西后，很艰难地抢到一间房。 下午练习赛，题水，所以去学习其他的东西了。 Day2 结...
APIO2018酱油记
TimeTraveller-VictoryCzt
05-12 553
对于我这种各类比赛都围观吃瓜的群众(都不能叫选手)，居然误打误撞的来参加APIO2018，当然是来打酱油的。 和我一起的还有两位前辈大佬(DraZxlNDt和Jeff),无限Orz%%%%%%。 就是打酱油，也算是一场人生经历吧，简要记述一下吧。 Day0 报到日。（2018.5.10） 早上SCMY还有点下雨，从学校寝室起来，偷偷溜出校门，坐上车回家收拾了东西，8:20就来到了机场...
CTSC/APIO 2018 咕咕记
LYD729
05-15 521
Preface 省选结束后还能继续逃几天文化课的学习，还是很资瓷的 日期从5.6一直到5.14号，时间很长就不用Day x的形式了 本文大概是记录每天日程的流水账 5.5 坐车到广州机场附近住了一晚，一个人一间房体验极佳（除了经常有飞机飞过） 由于忘带电脑晚上基本上是看番度过的 5.6 早上的飞机， 5.7 5.8 5.9 5.10 5.11 5...
HDU 5405 Sometimes Naive 题解(HLD+BIT)
qq_42925924的博客
08-25 1202
一题非常巧妙的DS题。 假设1为根。 首先可以发现，把这条连从图中抠掉，那么图中会剩下几个联通快:S1,S2,S3..SkS_1,S_2,S_3..S_kS1​,S2​,S3​..Sk​则Ans=(∑i=1nai)2−∑i=1k(si2)Ans=(\sum_{i=1}^n a_i)^2-\sum_{i=1}^k(s_i^2)Ans=(∑i=1n​ai​)2−∑i=1k​(si2​)。 以下所将的答案为∑i=1k(si2)\sum_{i=1}^k(s_i^2)∑i=1k​(si2​). 这题的难点就在于怎样维
APIO2_VDD怎么设置GPIO
最新发布
04-24
根据提供的引文，APIO2_VDD的GPIO设置如下： ```shell + vcc5v0_otg: vcc5v0-otg-regulator { + ... + gpio = ; + ... + }; ``` 其中，`gpio1 3`表示GPIO引脚编号为1的GPIO控制器上的第3个引脚，`GPIO_ACTIVE_HIGH`...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值