洛谷 P4093: bzoj 4553: [HEOI2016/TJOI2016]序列

题目传送门：洛谷P4093。

题意简述：

给定一个长度为 \(n\) 的序列 \(a\)。

同时这个序列还可能发生变化，每一种变化 \((x_i,y_i)\) 对应着 \(a_{x_i}\) 可能变成 \(y_i\)。

不会同时发生两种变化。

需要找出一个最长的子序列，使得这个子序列在任意一种变化下都是不降的。

只需要求出这个子序列的长度即可。

注意：可以不发生任何变化。

题解：

记 \(f[i]\) 为以第 \(i\) 项结尾的子序列最长长度。

则有转移：\(f[i]=\max_{j<i}(f[j])+1\)，同时还要满足 \(maxval_j\le a_i\) 和 \(a_j\le minval_i\)。

按照项从小到大转移，形成了天然的时间顺序，同时还要满足两个偏序限制。
其中 \(maxval_i\) 表示第 \(i\) 项最大能变成的值，\(minval_i\) 表示第 \(i\) 项最小能变成的值。

算上时间顺序，这是一个三维偏序问题，用 CDQ 分治 + 数据结构（我用了树状数组）就能解决。

 1 #include <cstdio>
 2 #include <algorithm>
 3 using namespace std;
 4 
 5 const int MN = 100005;
 6 const int MC = 100000;
 7 
 8 int N, M;
 9 int A[MN], Mx[MN], Mn[MN];
10 int f[MN], Ans;
11 int p[MN];
12 inline bool cmp1(int i, int j) { return Mx[i] < Mx[j]; }
13 inline bool cmp2(int i, int j) { return A[i] < A[j]; }
14 
15 int B[MN];
16 inline void Ins(int i, int x) { for (; i <= MC; i += i & -i) B[i] = max(B[i], x); }
17 inline void Clr(int i) { for (; i <= MC; i += i & -i) B[i] = 0; }
18 inline int Qur(int i) { int A = 0; for (; i; i -= i & -i) A = max(A, B[i]); return A;}
19 
20 void CDQ(int lb, int rb) {
21     if (lb == rb) {
22         f[lb] = max(f[lb], 1);
23         return;
24     }
25     int mid = lb + rb >> 1;
26     CDQ(lb, mid);
27     for (int i = lb; i <= rb; ++i)
28         p[i] = i;
29     sort(p + lb, p + mid + 1, cmp1);
30     sort(p + mid + 1, p + rb + 1, cmp2);
31     int j = lb;
32     for (int i = mid + 1; i <= rb; ++i) {
33         while (j <= mid && Mx[p[j]] <= A[p[i]]) {
34             Ins(A[p[j]], f[p[j]]);
35             ++j;
36         }
37         f[p[i]] = max(f[p[i]], Qur(Mn[p[i]]) + 1);
38     }
39     for (int i = lb; i <= mid; ++i)
40         Clr(A[i]);
41     CDQ(mid + 1, rb);
42 }
43 
44 int main() {
45     int x, y;
46     scanf("%d%d", &N, &M);
47     for (int i = 1; i <= N; ++i)
48         scanf("%d", &A[i]),
49         Mx[i] = Mn[i] = A[i];
50     for (int i = 1; i <= M; ++i)
51         scanf("%d%d", &x, &y),
52         Mx[x] = max(Mx[x], y),
53         Mn[x] = min(Mn[x], y);
54     CDQ(1, N);
55     for (int i = 1; i <= N; ++i)
56         Ans = max(Ans, f[i]);
57     printf("%d\n", Ans);
58     return 0;
59 }

 

转载于:https://www.cnblogs.com/PinkRabbit/p/10096763.html