1 条题解
-
0
翻译自 ROI 2023 D1T3。
如果对于所有 ,都有 ,则称 为峰值。
如果对于所有 ,都有 ,则称 为反峰值。
给定大小为 的排列 。需要将其分为两个非空子序列 和 。每个元素 必须恰好被分到一个子序列中。你需要最大化 中的峰值数量和 中的反峰值数量之和。
。
我们按顺序考虑每一个 。由于我们要把序列划分成两个子序列,分别是前缀最大和前缀最小,所以我们定义 表示当前考虑到第 个数的前缀最大序列中最后一个值, 表示当前考虑到第 个数的前缀最小序列中最后一个值。显然,我们可以通过储存这两个值来进行 DP。设 表示当前考虑到第 个数的 和 值,则转移是朴素的。
我们想要更快的算法,于是我们考虑研究 DP 转移的过程(也是按顺序对 分配的过程)。对于每一个数 ,当前的 我们已经知道了:
- 对于 来说:
- 如果 ,则我们将 加到 数列里会更新,加到 数列里没影响。
- 如果 ,则我们将 加到 数列里没影响,加到 数列里会更新。
- 如果 ,则加到任何一个数列都会更新。
- 对于 来说:
- 如果 ,则我们将 加到 数列里会更新,加到 数列里没影响。
- 如果 ,则我们将 加到 数列里没影响,加到 数列里会更新。
- 如果 ,则加到任何一个数列都没有影响。
对于两种情况,显然只有 在中间的时候情况不同。似乎如果更新少一点更简单,所以我们考虑 。
我们考虑 。对于这种情况,我们发现,对于每个数,如果它可以但是不用在 或者 数列中的话,它可以选择不影响答案。因此,对于所有后面 的数,我可以选择一些来构成最长上升子序列,其他的划分到 去,不对答案造成影响。 的数同理。容易发现,这其实是 和 问题,而这个问题我们可以用 轻松计算出每个 开头的 值。
但是对于 ,由于 的存在,并不是每个数都能选择不影响序列。如果我们直接考虑研究 的转移过程,那是没有任何性质的。但是,如果我们观察答案整体的结构,我们发现 和 分别是这个序列中的前缀和后缀部分,所以不难想到枚举分界点。
(为了方便区分,我们后文将 的部分称为前面的部分, 的部分称为后面的部分)
我们考虑枚举分界点 。如果我们要求后面的部分,它的答案,那么由 ,我们需要知道 和 的值来划分后面的数分别可以分到哪一个数列。因此,我们考虑枚举在分界点 选择完之后, 的值的。(此时,)
接着,我们考虑前面的部分怎么求。首先,由 ,前面的操作过程没什么救,但是我们知道前面 啊!由于所有 序列的所有数都 , 序列的所有数都 ,而 ,所以 序列的所有值都小于 序列的所有值!因此,我们再考虑枚举分界点 还没选择前, 的值。这样,我们可以预处理出前面所有 的数的前缀最大值个数,和所有 的前缀最小值的个数。
当然,实际上你在枚举 的时候,有可能存在 ,使得 。这种情况显然不合法!因此,我们枚举的 和 需要保证不存在 使得 。
但这样复杂度直接起飞,非常不好。考虑我们都枚举了什么东西。
- 分界点
- 分界点选择之前 值
- 分界点选择之后 值
我们发现,分界点选择之前的 和之后的 肯定有很多重复啊!假设 选择之前是 ,那 选择,肯定要么更新 ,要么更新 ,然后都是变成 ,所以实际上我们只需要枚举 ,然后分类讨论一下即可。
我们明明已经找到了这么多性质,时间复杂度还是 。怎么回事呢?但其实不然。由 ,我们知道 肯定是前面的 的值,而我们不能存在 在 之间,所以实际上 是 的后继。因此,你只需要枚举 两个值即可,复杂度为 。
我们考虑这个算法的细节。首先,枚举 ,则 是很容易通过维护 算出来的。然后,我们发现对于一个 ,肯定在前面出现过,而且这个 后面的数是没用的,所以我们可以对于每个 ,预处理出所有 的数的前缀最大值个数, 同理。对于后面的部分,我们可以通过 ,求出以每个 为开头的 。然后,由于我们要求以 为开头的最长子序列,和以 或者 开头的最长子序列,因此我们按顺序扫一下 就可以快速求出。
我们考虑优化这个算法。
我们大概的思路是:枚举一个值,另一个值数据结构维护,所以我们可以枚举 。结果我们发现,枚举 完全不可做!!!???这怎么办?
我们考虑枚举 为什么不优秀。枚举 ,给我们提供了 或 的信息,让我们知道从 开始进入后面的部分,但是我们不知道 ;不知道 后面应该找什么数(这里就很离谱了:我们不知道 ,还要对于 来二维数点!),但是知道 后面一定是 。如果枚举 的话,如果 在 的右边,那我们不能知道 是什么。
我们对于 维护的信息是不平衡的,但是我们又要通过 求出来 。于是,我们应该要枚举一个点,这个点左边涵盖了 ,右边是我们要找的转折点 。如果 在 左边,这个 可以是 ,否则可以是 ……对!我们可以枚举 和 最靠右的那个点!
假设我们枚举了 ,其中 都在 左边,且其中有一个是 。由 ,这两个值在 左边的值域上肯定是相邻的,所以我们可以通过 求出来另一个 。那么,我们分类讨论 是 还是 ,对于每种情况,我们还要考虑是 还是 的更新造成了 的大小关系变化。因此,我们枚举 ,然后分类讨论四种情况。不妨设 是 ,我们算出来了 , 跳完后超过了 ,则跳的下标 要满足 ,然后求 的最大值,这是经典的二位数点。
这份代码是 @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
- 上传者