【BZOJ5286】[HNOI2018]转盘(线段树)
题面
题解
很妙的一道题目啊。(全世界除了我这题都有40分,就我是一个状压选手
首先来发现一些性质,我们走一圈一定不会更差。
为啥呢?我们反过来看,我们可以钦定一个时间(T),然后从这个时刻出发,每个时刻可以向前走一步或者停留于此,而每个物品有一个消失时间,过了这个时间你还没有到这个位置你就凉了。
那么我们发现我们显然只需要走一圈就可以拿到所有的东西,如果走一圈还有东西拿不到那你走再多圈也拿不到。
那么现在我们要做的就是让钦定的时间(T)最短,那么我们又发现了,你显然不会停,停下来显然不优,因为我们现在要做的就是走一圈,既然所有东西都会消失,那么停下来干啥啊。
行,那么我们再正回来,也就是两个性质:只要我们开始走了,我们就不会停,要停也只会在起点停留若干时间之后一直向前走。另外一个是,我们只会走一圈,走回起点的前一个格子就结束了。
所以,不难发现我们的答案就是停留时间加上(n-1)。
考虑为啥要停留,证明有一个物品出现的时间很晚,我们必须要在所有物品出现以后才能拿。
那么这个出发时间是什么呢?假设物品的位置是(j),出发点是(i),(jge i),物品出现的时间是(t)。假设我们等待的时间是(x)。那么显然(x+j-ige t_j),移项可以得到(xge t_j-j+i),也就是(x=max(j-i+t_j)),写得好看点就是(x=max(t_j-j+i))。
那么我们要求的答案就是((n-1)+min_{i=1}^nmax_{j=i}^{n+n}t_j-j+i)
为什么要到(n+n)我们破环成链之后再原数组后面再接了一份。至于为什么不是到(i+n)是因为(i+n+1)到(n+n)这一段一定不会出现最大值。
现在考虑这个东西怎么维护就好了。
不难发现可能的最大值一定是一个关于(t_j-j)的单调栈。维护单调栈可以参考楼房重建那题。
大致的做法是,令(mx)表示区间的(t_j-j)的最大值。(ans)表示区间的答案。
考虑如何合并两个区间,这个东西是一个后缀区间的单调栈,所有右区间的值是直接拿过来用的。
考虑左区间接过来的答案,我们记录当前的右区间的最值,在左区间上面进行二分,找到最后一个大于最值的位置,那么单调栈我们就一直了,那么答案只需要沿着二分区间一路取(max)就可以了。
时间复杂度两个(log)。
#include<iostream>
#include<cstdio>
using namespace std;
#define MAX 200200
#define lson (now<<1)
#define rson (now<<1|1)
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,m,typ,T[MAX],lans;
int t[MAX<<2],mx[MAX<<2];
int Query(int now,int l,int r,int x)
{
if(l==r)return l+max(x,mx[now]);
int mid=(l+r)>>1;
if(mx[rson]>=x)return min(t[now],Query(rson,mid+1,r,x));
return min(Query(lson,l,mid,x),mid+1+x);
}
void pushup(int now,int l,int r)
{
int mid=(l+r)>>1;
t[now]=Query(lson,l,mid,mx[rson]);
mx[now]=max(mx[lson],mx[rson]);
}
void Build(int now,int l,int r)
{
if(l==r){t[now]=T[l];mx[now]=T[l]-l;return;}
int mid=(l+r)>>1;
Build(lson,l,mid);Build(rson,mid+1,r);
pushup(now,l,r);
}
void Modify(int now,int l,int r,int p)
{
if(l==r){t[now]=T[l];mx[now]=T[l]-l;return;}
int mid=(l+r)>>1;
if(p<=mid)Modify(lson,l,mid,p);
else Modify(rson,mid+1,r,p);
pushup(now,l,r);
}
int main()
{
n=read();m=read();typ=read();
for(int i=1;i<=n;++i)T[i]=T[i+n]=read();
Build(1,1,n+n);printf("%d
",lans=t[1]+n-1);
while(m--)
{
int x=read(),y=read();if(typ)x^=lans,y^=lans;
T[x]=T[x+n]=y;Modify(1,1,n+n,x);Modify(1,1,n+n,x+n);
printf("%d
",lans=t[1]+n-1);
}
return 0;
}