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. 树状数组的应用
我们需要一个数据结构,它能够:
- 在某个位置
pos添加一个元素。 - 查询区间
[L, R]中元素的个数。
树状数组 (Fenwick Tree) 非常适合这个任务。我们可以用一个大小为 n 的树状数组来维护砖块的位置。
4. 算法流程
-
数据结构化:
- 将
n个砖块存为一个结构体数组blocks,每个元素包含(height, position)。 - 将
m个查询存为一个结构体数组queries,每个元素包含(H, L, R, original_index)。
- 将
-
排序:
- 对
blocks数组按height升序排序。 - 对
queries数组按H升序排序。
- 对
-
处理查询:
- 初始化一个大小为
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]。
- 对于当前查询
- 初始化一个大小为
-
输出:
- 所有查询处理完毕后,按
0到m-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;
}