1 条题解

  • 0
    @ 2026-5-3 7:52:44

    Part.1\mathbf{Part. -1}

    翻译自 ROI 2023 D1T3

    如果对于所有 1j<i1 \le j < i,都有 aj<aia_j < a_i,则称 aia_i 为峰值。

    如果对于所有 1j<i1 \le j < i,都有 aj>aia_j > a_i,则称 aia_i 为反峰值。

    给定大小为 nn 的排列 p1,p2,,pnp_1,p_2,\dots,p_n。需要将其分为两个非空子序列 qqrr。每个元素 pp 必须恰好被分到一个子序列中。你需要最大化 qq 中的峰值数量和 rr 中的反峰值数量之和。

    1n2×1051 \leq n \leq 2 \times 10^5

    Part.1\mathbf{Part. 1}

    我们按顺序考虑每一个 aia_i。由于我们要把序列划分成两个子序列,分别是前缀最大和前缀最小,所以我们定义 mxmx 表示当前考虑到第 ii 个数的前缀最大序列中最后一个值,mnmn 表示当前考虑到第 ii 个数的前缀最小序列中最后一个值。显然,我们可以通过储存这两个值来进行 DP。设 fi,mx,mnf_{i, mx, mn} 表示当前考虑到第 ii 个数的 mxmxmnmn 值,则转移是朴素的。

    我们想要更快的算法,于是我们考虑研究 DP 转移的过程(也是按顺序对 aia_i 分配的过程)。对于每一个数 aia_i,当前的 mx,mnmx,mn 我们已经知道了:

    • 对于 mx<mnmx < mn 来说:
      • 如果 mx<mn<aimx < mn < a_i,则我们将 aia_i 加到 mxmx 数列里会更新,加到 mnmn 数列里没影响。
      • 如果 ai<mx<mna_i < mx < mn,则我们将 aia_i 加到 mxmx 数列里没影响,加到 mnmn 数列里会更新。
      • 如果 mx<ai<mnmx < a_i < mn,则加到任何一个数列都会更新。(1)^{{\color{red}(1)}}
    • 对于 mn<mxmn < mx 来说:
      • 如果 mn<mx<aimn < mx < a_i,则我们将 aia_i 加到 mxmx 数列里会更新,加到 mnmn 数列里没影响。
      • 如果 ai<mn<mxa_i < mn < mx,则我们将 aia_i 加到 mxmx 数列里没影响,加到 mnmn 数列里会更新。
      • 如果 mn<ai<mxmn < a_i < mx,则加到任何一个数列都没有影响。

    对于两种情况,显然只有 aia_i 在中间的时候情况不同。似乎如果更新少一点更简单,所以我们考虑 mn<mxmn < mx

    Part.2\mathbf{Part. 2}

    我们考虑 mn<mxmn < mx。对于这种情况,我们发现,对于每个数,如果它可以但是不用在 mxmx 或者 mnmn 数列中的话,它可以选择不影响答案。因此,对于所有后面 >mx> mx 的数,我可以选择一些来构成最长上升子序列,其他的划分到 mnmn 去,不对答案造成影响。<mn< mn 的数同理(2)^{\color{red}(2)}。容易发现,这其实是 LIS\mathbf{LIS}LCS\mathbf{LCS} 问题,而这个问题我们可以用 BIT\mathbf{BIT} 轻松计算出每个 aia_i 开头的 LIS,LCS\mathbf{LIS,LCS} 值。

    但是对于 mx<mnmx < mn,由于 (1)\color{red}(1) 的存在,并不是每个数都能选择不影响序列。如果我们直接考虑研究 mx<mnmx < mn 的转移过程,那是没有任何性质的(3)^{\color{red}(3)}。但是,如果我们观察答案整体的结构,我们发现 mx>mnmx > mnmn>mxmn > mx 分别是这个序列中的前缀和后缀部分,所以不难想到枚举分界点

    (为了方便区分,我们后文将 mx<mnmx < mn 的部分称为前面的部分,mx<mxmx < mx 的部分称为后面的部分)

    我们考虑枚举分界点 xx。如果我们要求后面的部分,它的答案,那么由 (2)\color{red}(2),我们需要知道 mnmnmxmx 的值来划分后面的数分别可以分到哪一个数列。因此,我们考虑枚举在分界点 xx 选择完之后,mn,mxmn,mx 的值的。(此时,mn<mxmn < mx

    接着,我们考虑前面的部分怎么求。首先,由 (3)\color{red}(3),前面的操作过程没什么救,但是我们知道前面 mx<mnmx < mn 啊!由于所有 mxmx 序列的所有数都 mx\leq mxmnmn 序列的所有数都 mn\geq mn,而 mx<mnmx <mn,所以 mxmx 序列的所有值都小于 mnmn 序列的所有值!因此,我们再考虑枚举分界点 xx 还没选择前,mn,mxmn,mx 的值。这样,我们可以预处理出前面所有 mx\leq mx 的数的前缀最大值个数,和所有 mn\geq mn 的前缀最小值的个数。

    当然,实际上你在枚举 mn,mxmn,mx 的时候,有可能存在 aia_i,使得 mx<ai<mnmx < a_i < mn。这种情况显然不合法!因此,我们枚举的 mxmxmnmn 需要保证不存在 aia_i 使得 mx<ai<mn(4)mx < a_i < mn^{\color{red}{(4)}}

    但这样复杂度直接起飞,非常不好。考虑我们都枚举了什么东西。

    • 分界点 xx
    • 分界点选择之前 mx,mnmx,mn
    • 分界点选择之后 mx,mnmx,mn

    我们发现,分界点选择之前的 mx,mnmx,mn 和之后的 mx,mnmx,mn 肯定有很多重复啊!假设 xx 选择之前是 mx,mnmx,mn,那 xx 选择,肯定要么更新 mxmx,要么更新 mnmn,然后都是变成 axa_x,所以实际上我们只需要枚举 mx,mnmx,mn,然后分类讨论一下即可。

    我们明明已经找到了这么多性质,时间复杂度还是 O(n3)\mathcal{O}(n^3)。怎么回事呢?但其实不然。由 (4)\color{red}(4),我们知道 mx,mnmx,mn 肯定是前面的 aia_i 的值,而我们不能存在 aia_imx,mnmx,mn 之间,所以实际上 mxmxmnmn 的后继。因此,你只需要枚举 mx,xmx,x 两个值即可,复杂度为 O(n2)\mathcal{O}(n^2)

    我们考虑这个算法的细节。首先,枚举 x,mxx,mx,则 mnmn 是很容易通过维护 set\mathbf{set} 算出来的。然后,我们发现对于一个 mxmx,肯定在前面出现过,而且这个 mxmx 后面的数是没用的,所以我们可以对于每个 aimxa_i \leq mx,预处理出所有 mx\leq mx 的数的前缀最大值个数,mnmn 同理。对于后面的部分,我们可以通过 BIT\mathbf{BIT},求出以每个 aia_i 为开头的 LIS,LCS\mathbf{LIS,LCS}。然后,由于我们要求以 xx 为开头的最长子序列,和以 i>x,ai<mni > x,a_i < mn 或者 i>x,ai>mxi > x, a_i > mx 开头的最长子序列,因此我们按顺序扫一下 mxmx 就可以快速求出。

    Part.3\mathbf{Part. 3}

    我们考虑优化这个算法。

    我们大概的思路是:枚举一个值,另一个值数据结构维护,所以我们可以枚举 xx。结果我们发现,枚举 xx 完全不可做!!!???这怎么办?

    我们考虑枚举 xx 为什么不优秀。枚举 xx,给我们提供了 mx<mn<axmx < mn < a_xax<mx<mna_x < mx < mn 的信息,让我们知道从 ii 开始进入后面的部分,但是我们不知道 mx,mnmx,mn;不知道 mnmn 后面应该找什么数(这里就很离谱了:我们不知道 mnmn,还要对于 x>i,ax<mnx > i,a_x < mn 来二维数点!),但是知道 mxmx 后面一定是 ii。如果枚举 mxmx 的话,如果 mnmnmxmx 的右边,那我们不能知道 mnmn 是什么。

    我们对于 mx,mnmx,mn 维护的信息是不平衡的,但是我们又要通过 mxmx 求出来 mnmn。于是,我们应该要枚举一个点,这个点左边涵盖了 mn,mxmn,mx,右边是我们要找的转折点 xx。如果 mxmxmnmn 左边,这个 xx 可以是 mxmx,否则可以是 mnmn……对!我们可以枚举 mnmnmxmx 最靠右的那个点!

    假设我们枚举了 xx,其中 mn,mxmn,mx 都在 xx 左边,且其中有一个是 axa_x。由 (4)\color{red}(4),这两个值在 xx 左边的值域上肯定是相邻的,所以我们可以通过 mxmx 求出来另一个 mnmn。那么,我们分类讨论 xxmxmx 还是 mnmn,对于每种情况,我们还要考虑是 mnmn 还是 mxmx 的更新造成了 mn,mxmn,mx 的大小关系变化。因此,我们枚举 xx,然后分类讨论四种情况。不妨设 axa_xmxmx,我们算出来了 mnmnmxmx 跳完后超过了 mnmn,则跳的下标 pp 要满足 p>x,ap>mnp > x, a_p > mn,然后求 LIS(p)\mathbf{LIS}(p) 的最大值,这是经典的二位数点。

    Part.+\mathbf{Part. +\infin}

    这份代码是 @Albert_Wei 写的,本人还没写,写了 update。

    int n, p[N], res[N][2], l[N], r[N];
    
    struct BIT {
        int tr[N];
        void init() { F(i, 1, n + 1) tr[i] = 0; }
        void modify(int x, int y) { for (; x <= n + 1; x += (x & -x)) tr[x] += y; }
        int query(int x) { int ans = 0; for (; x; x -= (x & -x)) ans += tr[x]; return ans; }
    } S;
    struct BIT2 {
        int tr[N];
        void init() { F(i, 1, n) tr[i] = 0; }
        void modify(int x, int y) { for (; x <= n; x += (x & -x)) Fmax(tr[x], y); }
        int query(int x) { int ans = 0; for (; x; x -= (x & -x)) Fmax(ans, tr[x]); return ans; }
    } U, V;
    
    void solve() {
        n = read();
        F(i, 1, n) p[i] = read();
        S.init();
        set<int> st;
        st.insert(0), st.insert(n + 1);
        F(i, 1, n) {
            l[i] = *(-- st.lower_bound(p[i])), r[i] = *st.upper_bound(p[i]);
            S.modify(l[i] + 1, 1), S.modify(r[i] + 1, -1);
            res[i][0] = S.query(p[i] + 1), res[i][1] = S.query(p[i]), st.insert(p[i]);
        }
        int ans = 0;
        U.init(), V.init();
        G(i, n, 1) {
            if (r[i] != n + 1) Fmax(ans, res[i][0] + U.query(n - r[i]) + V.query(r[i]));
            Fmax(ans, res[i][0] + U.query(n - p[i]) + V.query(p[i]));
            Fmax(ans, res[i][1] + U.query(n - p[i]) + V.query(p[i]));
            if (l[i]) Fmax(ans, res[i][1] + U.query(n - l[i]) + V.query(l[i]));
            U.modify(n - p[i] + 1, U.query(n - p[i]) + 1);
            V.modify(p[i], V.query(p[i]) + 1);
        }
        F(i, 0, n) Fmax(ans, U.query(n - i) + V.query(i));
        cout << ans << '\n';
    }
    
    • 1

    信息

    ID
    7336
    时间
    3000ms
    内存
    512MiB
    难度
    10
    标签
    递交数
    1
    已通过
    1
    上传者