题目
有一个数列,数列中某些项的值可能会变化,但同一个时刻最多只有一个值发生变化。
你已经得知所有可能的变化,请选出一个子序列,使得在任意一种变化中,这个子序列都是不降的。
思路
考虑动态规划。
设 f i f_i fi 表示以 i i i 结尾的不降子序列的长度, c i c_i ci 表示原先的数列, a i a_i ai 表示第 i i i 项最小可能变化到的值, b i b_i bi 表示第 i i i 项最大可能变化到的值。
则有 f i = max f j + 1 f_i=\max{f_j}+1 fi=maxfj+1,其中满足 j < i j<i j<i, b j ≤ c i b_j\le c_i bj≤ci 且 c j ≤ a i c_j\le a_i cj≤ai。
若暴力转移,时间复杂度
O
(
n
2
)
\mathcal{O}(n^2)
O(n2),时间会飞起,于是考虑使用某些毒瘤的数据结构优化。
仔细观察转移方程,我们想到了什么?
三维偏序: a j ≤ a i a_j\le a_i aj≤ai, b j ≤ b i b_j\le b_i bj≤bi 且 c j ≤ c i c_j\le c_i cj≤ci 的 j j j 的数量。
很像三维偏序对吧。
外层树状数组维护 b b b 属性,查询时取 1 ∼ c i 1 \sim c_i 1∼ci,内层权值线段树维护 c c c 属性,查询时取 1 ∼ a i 1 \sim a_i 1∼ai。
树状数组也不再维护对应 f f f 的加,而是维护 f f f 的最大值。
然后就是套板子了。
代码演示
#include <iostream>
#include <algorithm>
using namespace std;
int n, m, f[100010], a[100010], b[100010], c[100010];
int cnt, cnt1, tmp1[100010], root[100010];
struct node {
int lson, rson, val;
} tree[24170412];
int lowbit(int x) { return x & (-x); }
void pushup(int u) { tree[u].val = max(tree[tree[u].lson].val, tree[tree[u].rson].val); }
void update(int &u, int l, int r, int k, int x) {
if (!u) u = ++cnt;
if (l == r) {
tree[u].val = max(tree[u].val, x);
return;
}
int mid = (l + r) >> 1;
if (k <= mid) update(tree[u].lson, l, mid, k, x);
else update(tree[u].rson, mid + 1, r, k, x);
pushup(u);
}
void add(int x, int y, int z) {
for (int i = x; i <= n; i += lowbit(i)) update(root[i], 1, n, y, z);
}
int query(int l, int r, int ci) {
if (l == r) {
int sum = 0;
for (int i = 1; i <= cnt1; i++) sum = max(sum, tree[tmp1[i]].val);
return sum;
}
int mid = (l + r) >> 1;
if (ci <= mid) {
for (int i = 1; i <= cnt1; i++) tmp1[i] = tree[tmp1[i]].lson;
return query(l, mid, ci);
} else {
int sum = 0;
for (int i = 1; i <= cnt1; i++) sum = max(sum, tree[tree[tmp1[i]].lson].val);
for (int i = 1; i <= cnt1; i++) tmp1[i] = tree[tmp1[i]].rson;
return max(sum, query(mid + 1, r, ci));
}
}
int find(int bi, int ci) {
cnt1 = 0;
for (int i = bi; i; i -= lowbit(i))
tmp1[++cnt1] = root[i];
return query(1, n, ci);
}
signed main() {
cin >> n >> m;
for (int i = 1; i <= n; i++) {
cin >> c[i];
a[i] = b[i] = c[i];
}
for (int i = 1; i <= m; i++) {
int x, y;
cin >> x >> y;
a[x] = min(a[x], y);
b[x] = max(b[x], y);
}
for (int i = 1; i <= n; i++) {
f[i] = find(c[i], a[i]) + 1;
add(b[i], c[i], f[i]);
}
int ans = 0;
for (int i = 1; i <= n; i++)
ans = max(ans, f[i]);
cout << ans << endl;
return 0;
}
以下是暴力 dp 50pts 代码:
#include <iostream>
#include <algorithm>
using namespace std;
// f[i] 表示以 i 结尾的不降子序列的长度
int n, m, f[100010], a[100010], b[100010], c[100010];
signed main() {
cin >> n >> m;
for (int i = 1; i <= n; i++) {
cin >> c[i];
a[i] = b[i] = c[i];
}
for (int i = 1; i <= m; i++) {
int x, y;
cin >> x >> y;
a[x] = min(a[x], y);
b[x] = max(b[x], y);
}
for (int i = 1; i <= n; i++) {
f[i] = 1;
for (int j = i - 1; j >= 1; j--)
if (b[j] <= c[i] && c[j] <= a[i])
f[i] = max(f[i], f[j] + 1);
}
int ans = 0;
for (int i = 1; i <= n; i++) ans = max(ans, f[i]);
cout << ans << endl;
return 0;
}