「CEOI2010」 MP3Player 题解

「CEOI2010」 MP3Player

题意

\(~~~~\) 给出 \(n\) 次对音量的操作,每次操作为使音量 \(+1\) 或 \(-1\) ,同时若音量操作后超出 \([0,V_{\max}]\) 或该操作与上一个操作的时间间隔 \(>t\) 则该操作无效(不论上一个操作是否有效)。已知若干次操作后的最终音量,求 \(t\) 可能的最大值以及该最大值下可能的最大的初始音量。

\(~~~~\) \(1\leq n\leq 2\times 10^5,2\leq V_{\max}\leq 5000\)

题解

\(~~~~\) 线段树可以维护的东西增加了,线段树yyds

\(~~~~\) 先给出一个暴力通用思路:枚举所有在操作中出现的 \(t\) ,则对于任意一个初始音量 \(V_1\),我们都可以 \(\mathcal{O(n)}\)求出其对应的结束音量 \(V_2\) 。并且输出出来我们可以发现若将 \(V_2\) 视作 \(V_1\) 的函数,则该函数不降,故可以二分求出是否存在 \(V_1\) 在该 \(t\) 下使得结束音量为 \(V_2\) 。时间复杂度:\(\mathcal{O(n^2\log n)}\) 似乎一个subtask都过不了(

\(~~~~\) 考虑继续发掘一下该函数的性质,可以发现该函数存在 \(A,B\) ,使得 \([0,A)\) 内所有 \(V_1\) 的 \(V_2\) 相同,\([A,B)\) 内所有 \(V_1\) 的 \(V_2\) 每个相对前一个多 \(1\) ,\([B,V_{\max}]\) 内所有 \(V_1\) 的 \(V_2\) 也相同。简单来说该函数依次是:常函数、一次函数、常函数。理性理解一下, \(V_1\) 在从 \(0\) 开始增长时中途一定会有一段时间卡在下界、增长到 \(V_{max}\) 前有一段时间卡在上界,当然这段时间可以为空。那么显然,该函数的值域也一定是连续的整数区间。

\(~~~~\) 由此,我们可以考虑维护两个常函数的值 \(a,b\) 和一次函数的 \(y=x+c\) 的 \(c\) 。则当我们知道两段相连区间的 \(a,b,c\) 时,我们就可以考虑复合两个函数得到这个大区间的 \(a,b,c\) 。稍微手玩一下就可以推出 \(\mathcal{O(1)}\) 转移,这点可以见代码。自然而然地,我们想到用线段树维护区间的函数值。

\(~~~~\) 因此,我们仍然保留枚举所有可能的 \(t\) ,然后每次在线段树上删除某个点对应的贡献,并检查该函数值域中是否有 \(V_2\) 。该过程中每个点只会被删去一次,所以时间复杂度 \(\mathcal{O(n\log n)}\) 。

代码

查看代码
#include <cstdio>
#include <algorithm>
using namespace std;
int n,Vmax,V2;
int val[200005],Time[200005];
struct node{
	int val,Time,id;
	node(){}
	node(int V,int T,int I){val=V,Time=T,id=I;}
}P[200005];
struct Seg{
	#define ls p<<1
	#define rs p<<1|1
	#define lson p<<1,l,mid
	#define rson p<<1|1,mid+1,r
	struct Func{
		int a,b,c;
		Func(){}
		Func(int A,int B,int C){a=A,b=B,c=C;}
		int F(int x){return min(b,max(a,x+c));}
	}tr[800005];
	//a一平的值 b二平的值 c一次函数 y=x+c 
	void pushUp(int p)
	{
		tr[p].a=min(tr[rs].b,max(tr[rs].a,tr[ls].a+tr[rs].c));
		tr[p].b=max(tr[rs].a,min(tr[rs].b,tr[ls].b+tr[rs].c));
		tr[p].c=tr[ls].c+tr[rs].c;
	}
	void Build(int p,int l,int r)
	{
		if(l==r)
		{
			tr[p]=Func(0,Vmax,val[l]);
			return;
		}
		int mid=(l+r)>>1;
		Build(lson);Build(rson);
		pushUp(p);
	}
	void Modify(int p,int l,int r,int aim)
	{
		if(l==r)
		{
			tr[p]=Func(0,Vmax,0);
			return;
		} 
		int mid=(l+r)>>1;
		if(aim<=mid) Modify(lson,aim);
		if(mid<aim)  Modify(rson,aim);
		pushUp(p);
	} 
}Seg;
int Calc()
{
	int l=0,r=Vmax,mid,Ans=0;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(Seg.tr[1].F(mid)<=V2) l=mid+1,Ans=mid;
		else r=mid-1;
	}
	return Ans;
}
bool cmp(node x,node y){return x.Time>y.Time;}
int main() {
	scanf("%d %d %d",&n,&Vmax,&V2);
	char S[10];
	for(int i=1;i<=n;i++)
	{
		scanf("%s %d",S+1,&Time[i]);
		val[i]=(S[1]=='+')?1:-1;
	}n--;
	for(int i=1;i<=n;i++)
	{
		val[i]=val[i+1];
		Time[i]=Time[i+1]-Time[i];
		P[i]=node(val[i],Time[i],i);
	}
	Seg.Build(1,1,n);
	sort(P+1,P+1+n,cmp);
	if(Seg.tr[1].F(0)<=V2&&V2<=Seg.tr[1].F(Vmax)) return puts("infinity")&0;
	int j;
	for(int i=1;i<=n;i=j)
	{
		int Tmp=P[i].Time;
		for(j=i;P[j].Time==P[i].Time;j++) Seg.Modify(1,1,n,P[j].id);
		if(Seg.tr[1].F(0)<=V2&&V2<=Seg.tr[1].F(Vmax)) return printf("%d %d\n",P[i].Time-1,Calc())&0;
	}
	return 0;
}
上一篇:「题解」Codeforces 1372E Omkar and Last Floor


下一篇:「NOI2020」 美食家 【矩阵快速幂】