【BZOJ5286】[HNOI2018]转盘(线段树)
题面
题解
很妙的一道题目啊。(全世界除了我这题都有40分,就我是一个状压选手
首先来发现一些性质,我们走一圈一定不会更差。 为啥呢?我们反过来看,我们可以钦定一个时间\(T\),然后从这个时刻出发,每个时刻可以向前走一步或者停留于此,而每个物品有一个消失时间,过了这个时间你还没有到这个位置你就凉了。 那么我们发现我们显然只需要走一圈就可以拿到所有的东西,如果走一圈还有东西拿不到那你走再多圈也拿不到。 那么现在我们要做的就是让钦定的时间\(T\)最短,那么我们又发现了,你显然不会停,停下来显然不优,因为我们现在要做的就是走一圈,既然所有东西都会消失,那么停下来干啥啊。 行,那么我们再正回来,也就是两个性质:只要我们开始走了,我们就不会停,要停也只会在起点停留若干时间之后一直向前走。另外一个是,我们只会走一圈,走回起点的前一个格子就结束了。 所以,不难发现我们的答案就是停留时间加上\(n-1\)。 考虑为啥要停留,证明有一个物品出现的时间很晚,我们必须要在所有物品出现以后才能拿。 那么这个出发时间是什么呢?假设物品的位置是\(j\),出发点是\(i\),\(j\ge i\),物品出现的时间是\(t\)。假设我们等待的时间是\(x\)。那么显然\(x+j-i\ge t_j\),移项可以得到\(x\ge 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#include 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\n",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\n",lans=t[1]+n-1); } return 0;}