【ACWing】592. 雨

题目地址：

https://www.acwing.com/problem/content/594/

海里有一个小岛。岛可以描述为一个具有$R$行$C$列的矩阵，其中包含的数字$H[i][j]$表示每个单元格的高度。以下是一个$3\times 3$岛屿的示例：

3 5 5
5 4 5
5 5 5

有时，大雨会在这个岛上的每个单元格上均匀地降落。你可以假设降水量可以是无限大。在这样的大雨之后，岛上的一些区域（由沿着边缘连接的一个或多个单元格形成）可能会存在积水。一片区域可以存在积水的前提是与该区域内单元格沿着边缘连接的区域外的单元格的高度均高于区域内单元格。周围的海洋被视为高度为$0$的无限单元格。否则，水将始终流入其他区域并最终出海。你可以假设海的高度永远不会改变。在大雨过后，我们将使用$W[i][j]$来表示岛屿上每个单元格的高度。以下是对大雨过后的示例岛屿的高度的分析。最中央的初始高度为$4$的单元格被初始高度为$5$的单元格包围着，因此在该单元格将产生积水，在积水高度到达$5$以后，没有更多区域被高度更高的单元格包围，因此将不再继续积水。请注意，水不会直接在两个共享顶点的单元格之间流动；水必须沿着共同的边缘流动，所以中央单元格的积水不会从左上角处流走。以下是雨后示例岛屿的高度：

3 5 5
5 5 5
5 5 5

给定岛的高度矩阵，你能计算出大雨后的每个单元格的增加高度$(W[i][j]-H[i][j])$的总和吗？

输入格式：
第一行包含整数$T$，表示共有$T$组测试数据。每组数据第一行包含两个整数$R$和$C$。接下来$R$行每行包含$C$个整数$H[i][j]$，表示矩阵内单元格的高度。

输出格式：
每组数据输出一个结果，每个结果占一行。结果表示为Case #x: y，其中$x$是组别编号（从$1$开始），$y$是增加高度的总和。

数据范围：
$1\le T\le 100$
$1\le H[i][j]\le 1000$
$1\le R,C\le 50$

每个位置的水量是由其外围的最矮的柱子的高度减去其自己的柱子高度决定的（当然如果减出来是负数则视为$0$），所以可以用最小堆，先将最外围的柱子入堆，然后每出堆一个柱子，就看一下其四周哪个柱子的水量没算出过，累加其水量后，再将注水后的高度和坐标入堆（好像注水之后柱子长高了一样）。参考https://blog.csdn.net/qq_46105170/article/details/108575830。代码如下：

#include <iostream>
#include <cstring>
#include <queue>
using namespace std;

// 第一维存柱子高度，后面两维存坐标
typedef pair<int, pair<int, int> > PII;

const int R = 55, C = 55, d[] = {1, 0, -1, 0, 1};
int r, c;
int h[R][C];
bool st[R][C];

int get_vol() {
    memset(st, 0, sizeof st);
    
    priority_queue<PII, vector<PII>, greater<PII> > heap;
    for (int i = 1; i <= r; i++) {
        heap.push({h[i][1], {i, 1}}), heap.push({h[i][c], {i, c}});
        st[i][1] = st[i][c] = true;
    }
    for (int i = 1; i <= c; i++) {
        heap.push({h[1][i], {1, i}}), heap.push({h[r][i], {r, i}});
        st[1][i] = st[r][i] = true;
    }

    int res = 0;
    while (heap.size()) {
        auto t = heap.top(); heap.pop();
        int hei = t.first, x = t.second.first, y = t.second.second;
        
        for (int i = 0; i < 4; i++) {
            int nx = x + d[i], ny = y + d[i + 1];
            if (!st[nx][ny] && 1 <= nx && nx <= r && 1 <= ny && ny <= c) {
                res += max(0, hei - h[nx][ny]);
                heap.push({max(hei, h[nx][ny]), {nx, ny}});
                // 算出(nx, ny)这个位置的水量之后就标记一下
                st[nx][ny] = true;
            }
        }
    }
    
    return res;
}

int main() {
    int T;
    cin >> T;
    for (int t = 1; t <= T; t++) {
        cin >> r >> c;
        for (int i = 1; i <= r; i++)
            for (int j = 1; j <= c; j++)
                cin >> h[i][j];

        printf("Case #%d: %d\n", t, get_vol());
    }

    return 0;
}

时间复杂度$O(RC\log (R+C))$，空间$O(RC)$。