CodeForces_1348EF - Phoenix and Memory 贪心+线段树找区间最小值

题目要求在满足一定限制条件下找出解,通过贪心策略将限制按R值排序,从最小的R开始选择,并使用线段树查找区间最小值。当存在两个位置可交换时,满足特定条件,通过线段树检查交换可行性。若找不到可行交换则输出YES,否则输出包含两个解的序列。
摘要由CSDN通过智能技术生成

这题如果没有输出2个解就很简单。

是个之前做过的类型:

把所有限制按R排序,

然后每次取出R最小的,然后从其L开始选,尽量选能选的中最小的。

这样选如果能选完,就说明有解。

贪心正确性显然:R大的至少可以选则R做为点来用。所以按R升序遍历,每次优先选左边的,能让后边的可选的更多。

用set维护可选的数即可。

这题加了个输出2个方案。

我们考虑最简单的情况:即确定一个序列后,是否有2个位置,posi,posj,使得i,j可以交换。

显然必须满足    Lj<=posi<posj<Ri

我们枚举i,只需要找到任意一个j满足上式,就可以交换i,j然后输出。

我们可以先把所有的Lj存到以 posj为区间的线段树里,每次i,查询区间posi+1  -   Ri   之间的Lj最小值是否小于等于posi。

如果是则可以交换输出NO。加上2个序列

如果没找到,就输出YES。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
const double PI= acos(-1.0);
const int M = 2e5+7;
/*
int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt,val;}ee[M*2];
void add(int x,int y){ee[++cnt].nxt=head[x],ee[cnt].to=y,head[x]=cnt;}
*/
struct node{
	int l,r,id,pos;
	bool operator <(const node &x)const
	{
		return r<x.r;
	}
}p[M];
struct no{
	int x,id;
};
set<int>s;
int id[M],mp[M];
no mn[M<<2];
void bd(int o,int l,int r)
{
	if(l==r)
	{
		mn[o].x=1e6;
		mn[o].id=l;
		return ;
	}
	int m=(l+r)/2;
	bd(ls,l,m);
	bd(rs,m+1,r);
	if(mn[ls].x<mn[rs].x)mn[o]=mn[ls];
	else mn[o]=mn[rs];
}
void up(int o,int l,int r,int x,int d)
{
	if(l==r)
	{
		mn[o].x=d;
		mn[o].id=l;
		return ;
	}
	int m=(l+r)/2;
	if(x<=m)up(ls,l,m,x,d);
	else up(rs,m+1,r,x,d);
	if(mn[ls].x<mn[rs].x)mn[o]=mn[ls];
	else mn[o]=mn[rs];
}
no qu(int o,int l,int r,int x,int y)
{
	if(x<=l&&r<=y)
	{
		return mn[o];
	}
	int m=(l+r)/2;
	no t1={1000000+7,0},t2={1000000+7,0};
	if(x<=m)t1=qu(ls,l,m,x,y);
	if(y>m)t2=qu(rs,m+1,r,x,y);
//	cout<<"----     "<<t1.x<<"-"<<t1.id<<"  "<<t2.x<<"-"<<t2.id<<endl;
	if(t1.x<t2.x)return t1;
	else return t2;
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
  	int n;
  	cin>>n;
  	for(int i=1;i<=n;i++)
  	{
  		cin>>p[i].l>>p[i].r;
  		p[i].id=i;
  		s.insert(i);
	}
	bool f=true;
	sort(p+1,p+1+n);
	bd(1,1,n);
	for(int i=1;i<=n;i++)
	{
		auto it=s.lower_bound(p[i].l);//题目保证了合法 
		id[p[i].id]=*it;
		mp[*it]=p[i].id;
		p[i].pos=*it;
	//	cout<<" ---------------------------- "<<*it<<"   "<<p[i].l<<"  = "<<endl;
		up(1,1,n,*it,p[i].l);//单点更新 
		s.erase(it); 
	}
	
	for(int i=1;i<=n;i++)
	{
		int t=p[i].pos;
		no lj=qu(1,1,n,t+1,p[i].r);//区间最小值 ,+其下标 
		if(lj.x<=t)
		{
			cout<<"NO"<<endl;
	//		cout<<p[i].id<<"--  "<<p[i].pos<<"  "<<lj.id<<"-- "<<mp[lj.id]<<" "<<lj.x<<endl;
			for(int j=1;j<=n;j++)cout<<id[j]<<" ";cout<<endl;
			swap(id[p[i].id],id[mp[lj.id]]);
			for(int j=1;j<=n;j++)cout<<id[j]<<" ";cout<<endl;
			return 0;
		}
	}
	cout<<"YES"<<endl;
	for(int j=1;j<=n;j++)cout<<id[j]<<" ";cout<<endl;
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值