Super Mario

题目链接

HDU-4417 Super Mario

题目大意

给定一个长度为 n 的序列,代表每个位置的砖块高度。有 m 次询问,每次询问一个区间 [L, R] 和一个最大跳跃高度 H,要求回答在 [L, R] 这个区间内,有多少砖块的高度小于或等于 H

暴力代码

cpp
copy
        
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
#include <iostream> #include <vector> #include <cstdio> using namespace std; void solve(int caseNum) { int n, m; // 读取 n (道路长度) 和 m (查询次数) if (!(cin >> n >> m)) return; vector<int> h(n); // 读取每个砖块的高度 for (int i = 0; i < n; ++i) { cin >> h[i]; } printf("Case %d:\n", caseNum); // 处理 m 次查询 for (int i = 0; i < m; ++i) { int L, R, H; cin >> L >> R >> H; int count = 0; // 暴力核心:直接遍历 L 到 R 区间 // 题目保证 0 <= L <= R < n,所以不需要做边界检查 for (int j = L; j <= R; ++j) { if (h[j] <= H) { count++; } } printf("%d\n", count); } } int main() { int t; cin >> t; for (int i = 1; i <= t; ++i) { solve(i); } return 0; }

解题思路

这是一个典型的二维数点问题,查询涉及位置和数值两个维度。对于这类问题,如果可以离线处理,一个常见的优化技巧是将查询排序,从而将二维问题降为一维问题

1. 离线处理

我们可以将所有查询读入并存储起来,而不是立即回答。然后,我们按照查询的高度 H 进行升序排序。同时,我们也把所有的砖块按照高度进行升序排序。

2. 排序带来的好处

当我们按高度 H 从小到大的顺序处理查询时,我们会发现一个关键特性:对于一个查询 q_i,所有满足其高度限制 H_i 的砖块,对于后续任何查询 q_j (因为 H_j >= H_i) 也同样满足。

这意味着,我们可以维护一个数据结构,随着我们处理的查询 H 越来越大,我们将越来越多满足高度条件的砖块加入到这个数据结构中。

3. 树状数组的应用

我们需要一个数据结构,它能够:

  1. 在某个位置 pos 添加一个元素。
  2. 查询区间 [L, R] 中元素的个数。

树状数组 (Fenwick Tree) 非常适合这个任务。我们可以用一个大小为 n 的树状数组来维护砖块的位置。

4. 算法流程

  1. 数据结构化

    • n 个砖块存为一个结构体数组 blocks,每个元素包含 (height, position)
    • m 个查询存为一个结构体数组 queries,每个元素包含 (H, L, R, original_index)
  2. 排序

    • blocks 数组按 height 升序排序。
    • queries 数组按 H 升序排序。
  3. 处理查询

    • 初始化一个大小为 n 的树状数组 bit,所有元素为0。
    • 初始化一个 block_ptr = 0,用于指向 blocks 数组。
    • 遍历排序后的 queries 数组:
      • 对于当前查询 q,其高度为 H。我们先将所有高度小于等于 H 的砖块加入树状数组。具体操作是:while (block_ptr < n && blocks[block_ptr].height <= H),执行 bit.add(blocks[block_ptr].position, 1),然后 block_ptr++
      • 此时,树状数组中值为1的位置,就代表了所有高度 ≤ H 的砖块。
      • 查询 [L, R] 区间内的砖块数,即为 bit.query(R) - bit.query(L-1)
      • 将结果存入 ans[q.original_index]
  4. 输出

    • 所有查询处理完毕后,按 0m-1 的顺序输出 ans 数组中的结果。

通过这种方式,每个砖块和每个查询都只被处理常数次,树状数组的单次操作复杂度为 O(log n),排序的复杂度为 O(n log n + m log m)。因此,算法总的时间复杂度为 O(n log n + m log m),足以通过本题。

代码实现

cpp
copy
        
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
#include <algorithm> #include <bits/stdc++.h> using namespace std; const int maxn = 2e5+5; int n,m; int unique_cnt; int a[maxn],b[maxn]; int root[maxn]; // root[i] 第i颗树的根 struct Node { int l,r,sum; // sum 区间内权值(数字的数量) } tree[maxn << 5]; // 2^5 倍 int node_cnt = 0; int get_node() {return ++node_cnt;} // 建立线段树 void build() { } int update(int u,int l,int r,int id) { int rt = get_node(); tree[rt] = tree[u]; // 复制原来的结点 tree[rt].sum ++; // 多了一个值,增加1 if( l==r) return rt; // 叶子结点,返回 int mid = (l+r) >>1; if( id <= mid ) tree[rt].l = update(tree[u].l,l,mid,id); if( id > mid) tree[rt].r = update(tree[u].r,mid+1,r,id); return rt; } // u,v 对以两个树上对应的结点 // 本质是求区间和 int query(int u,int v,int l,int r,int k) { // 结点对应的区间,完全被包含 if( r <= k ) { return tree[v].sum - tree[u].sum; } if( l == r) { return tree[v].sum - tree[u].sum; } // 当前结点的左子树sum int mid = (l+r) >>1; int lsum = 0, rsum = 0; if( mid >= 1 ) lsum = query(tree[u].l,tree[v].l,l,mid,k); if( mid+1 <= k) rsum = query(tree[u].r,tree[v].r, mid+1, r, k); return lsum + rsum; } void init() { std::cin >> n >> m; for(int i = 1;i <= n ;++i ) // i: 1->n { std::cin >> a[i]; b[i] = a[i]; } sort(b+1,b+1+n); unique_cnt= std::unique(b+1,b+1+n)- (b+1);// 离散化 } int main (int argc, char *argv[]) { int T; std::cin >> T; int case_num = 1; while (T--) { printf("Case %d:\n", case_num++); node_cnt = 0; init(); // 建立n颗线段树 for(int i = 1;i <= n ;++i ) // i: 1->n { // 找到a[i] 对应的离散化后的值 int id = lower_bound(b+1, b+1+unique_cnt, a[i]) - b; // cout << id << " "; root[i] = update(root[i-1],1,unique_cnt,id); } // cout << endl; while (m--) { int x,y,k; std::cin >> x >> y >> k; x++;y++; // 得到k对应的离散化后的值 int kk = lower_bound(b+1, b+1+unique_cnt,k) - b; // 细节!!: 如果查询到的值就是k // 表明k 本身 离散化后存在 // 如果k 不在离散化后的值, if( b[kk] != k ) { kk--; } if( kk == 0) { // 特判 // 这里是一个细节 // 线段树里面的区间 [1,unique_cnt] // 如果要 查询 <=0,那么 就会出现错误 printf("0\n"); } else { int t = query(root[x-1], root[y],1,unique_cnt,kk); std::cout << t << "\n"; } } } return 0; }