[Luogu P4425] [BZOJ 5286] [HNOI AHOI2018]转盘

[HNOI/AHOI2018]转盘

洛谷传送门

BZOJ传送门

题目描述

一次小G和小H准备去聚餐，但是由于太麻烦了于是题面简化如下：

一个转盘上有摆成一圈的 $n$ 个物品（编号$1\sim n$ ），其中的 $i$ 个物品会在 $T_i$ 时刻出现。

在$0$时刻时，小G可以任选 $n$ 个物品中的一个，我们将其编号为 $s_0$ 。并且如果 $i$ 时刻选择了物品 $s_i$ ，那么 $i+1$时刻可以继续选择当前物品或选择下一个物品。当 $s_i$ 为 $n$ 时，下一个物品为物品 $1$ ，否则为物品 $s_{i}+1$ 。在每一时刻（包括$0$时刻），如果小G选择的物品已经出现了，那么小G将会标记它。小H想知道，在物品选择的最优策略下，小G什么时候能标记所有物品？

但麻烦的是，物品的出现时间会不时修改。我们将其描述为 $m$ 次修改，每次修改将改变其中一个物品的出现时间。每次修改后，你也需求出当前局面的答案。对于其中部分测试点，小H还追加了强制在线的要求。

输入输出格式

输入格式：

第一行三个非负整数 $n$ 、 $m$ 、 $p$ ，代表一共有 $n$ 个物品， $m$ 次修改。 $p$ 只有$0$或$1$两种取值，强制在线时 $p$ 为$1$，否则 $p$ 为$0$.

接下来一行，有 $n$ 个非负整数，第 $i$ 个数 $T_i$ 代表物品 $i$ 的出现时间。

接下来 $m$ 行，每行两个非负整数 $x$ 、 $y$ ，代表一次修改及询问。修改方式如下：

（1）如果 $p=0$ ，则表示物品 $x$ 的出现时间 $T_x$ 修改为 $y$ 。

（2）如果 $p=1$ ，在先将 $x$ 和 $y$ 分别异或 $LastAns$ ，得到 $x'$ 和 $y'$ ，然后将物品 $x'$的出现时间 $T_{x'}$ 修改为 $y'$ 。其中的 $LastAns$ 是前一个询问的结果；特别的，第一次修改时 $LastAns$ 为初始局面的答案。

保证输入合法。

输出格式：

第一行一个整数，代表初始局面的答案。

接下来 $m$ 行每行一个整数，分别代表每次修改后的答案。

输入输出样例

输入样例#1：

5 3 0
1 2 3 4 5
3 5
5 0
1 4

输出样例#1：

5
7
6
7

说明

【数据范围】

$3≤n≤10^5,0≤m≤10^5,0≤T_i/T_x≤10^5$ 。

解题分析

有个很显然的结论：我们肯定不会停在原地。

似乎不是很显然， 但我们可以倒着想：

一开始你在$t$时刻， 你可以选择停留在原地或向前走一步， 然后$t--$。

每个点上的物品会在$T_i$时消失， 求能得到所有物品的一个最小的$t$。

显然这样考虑我们就不可能停下来了…..

考虑破环为链， 倍长$T$数组。 我们在$[n,2n-1]$位置枚举一个终点$i$， 那么要求是：

$$t-(i-j)\ge T_j(j\in[i-n+1,i]) \\ t_{min}=max_{j=i-n+1}^{i}(T_j-j)+i$$
最终答案为

$$Ans = min_{i=n}^{2n-1}(max_{j=i-n+1}^{i}(T_j-j)+i)$$
这样就有了一个 $O(n^2log(N))$的优秀做法。

把$i$中间的$n$提出来， 设$A_i=T_i-i$，就有：

$$Ans=min_{i=1}^{n}(max_{j=i}^{i + n-1}(A_i)+i)+N-1$$
因为 $A_i=T_i-i>T_i-(i+n)=A_{i+n}$， 所以我们可以减少一些限制：
$$Ans=min_{i=1}^{n}(max_{j=i}^{2n}(A_i)+i)+N-1$$
这样等于我们只需要一个后缀的最大值就可以了， 但是并没有优化复杂度。

建立线段树， 设节点$P$管辖范围为$[L,R]$， $ans_P=min_{i=L}^{mid}(max_{j=i}^{R}(A_i)+i)$，$mx_P=max_{i=L}^{R}(A_i)$。

当$L=R$时， $ans_P=T_i,mx_P=A_i$。

现在考虑如何合并两个区间。

显然靠右的一个区间的最大值一定会被计算入左边的这一部分， 设与$P$合并的节点为$H$，因此实际上我们求的是：

$$ans_{fat[P]}=min_{i=L}^{R}(max(max_{j=i}^{R}(A_i), mx_H)+i)$$
用 $query(L,R,mx_H)$表示这一过程， 那么：

• 当$L=R$时，$ans_{fat[P]}=max(A_L,mx_H)+L$
• 当$mx_{rs[P]}\ge mx_H$时， 左半部分无关， 所以$ans_{fat[P]}=min(ans_P, query(mid+1,R,mx_H))$
• 否则当$i\in[mid+1,R]$时答案为$mid+1+mx_H$，另一部分递归到左区间处理， $ans_{fat[P]}=min(mid+1+mx_H,query(L,mid,mx_H))$。

代码如下：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <algorithm>
#include <assert.h>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 200050
#define ls (now << 1)
#define rs (now << 1 | 1)
template <class T>
IN void in(T &x)
{
    x = 0; R char c = gc;
    for (; !isdigit(c); c = gc);
    for (;  isdigit(c); c = gc)
    x = (x << 1) + (x << 3) + c - 48;
}
int T[MX], A[MX], dot, q, typ, lastans, bd;
struct Node {int mx, ans;} tree[MX << 2];
namespace SGT
{
    int query(R int now, R int lef, R int rig, R int lim)
    {
        if(lef == rig) return lef + std::max(lim, A[lef]);
        int mid = lef + rig >> 1;
        if(tree[rs].mx >= lim) return std::min(tree[now].ans, query(rs, mid + 1, rig, lim));
        return std::min(query(ls, lef, mid, lim), mid + 1 + lim);
    }
    IN void pushup(R int now, R int lef, R int rig)
    {
        tree[now].mx = std::max(tree[ls].mx, tree[rs].mx);
        tree[now].ans = query(ls, lef, lef + rig >> 1, tree[rs].mx);
    }
    void build(R int now, R int lef, R int rig)
    {
        if(lef == rig) return tree[now].ans = T[lef], tree[now].mx = A[lef], void();
        int mid = lef + rig >> 1;
        build(ls, lef, mid), build(rs, mid + 1, rig);
        pushup(now, lef, rig);
    }
    void modify(R int now, R int lef, R int rig, R int tar)
    {
        if(lef == rig) return tree[now].ans = T[lef], tree[now].mx = A[lef], void();
        int mid = lef + rig >> 1;
        if(tar > mid) modify(rs, mid + 1, rig, tar);
        else modify(ls, lef, mid, tar);
        pushup(now, lef, rig);
    }
}
int main(void)
{
    int a, b;
    in(dot), in(q), in(typ); bd = dot << 1;
    for (R int i = 1; i <= dot; ++i) in(T[i]), A[i] = T[i] - i;
    for (R int i = dot + 1; i <= bd; ++i) T[i] = T[i - dot], A[i] = A[i - dot] - dot;
    SGT::build(1, 1, bd); printf("%d\n", lastans = tree[1].ans + dot - 1);
    W (q--)
    {
        in(a), in(b); if(typ) a ^= lastans, b ^= lastans;
        T[a] = T[a + dot] = b;
        A[a] = T[a] - a; A[a + dot] = A[a] - dot;
        SGT::modify(1, 1, bd, a), SGT::modify(1, 1, bd, a + dot);
        printf("%d\n", lastans = tree[1].ans + dot - 1);
    }
}