三维空间bfs（L3-004 肿瘤诊断 （30 分））

原题链接
L3-004 肿瘤诊断 （30 分）
在诊断肿瘤疾病时，计算肿瘤体积是很重要的一环。给定病灶扫描切片中标注出的疑似肿瘤区域，请你计算肿瘤的体积。

输入格式：
输入第一行给出4个正整数：M、N、L、T，其中M和N是每张切片的尺寸（即每张切片是一个M×N的像素矩阵。最大分辨率是1286×128）；L（≤60）是切片的张数；T是一个整数阈值（若疑似肿瘤的连通体体积小于T，则该小块忽略不计）。

最后给出L张切片。每张用一个由0和1组成的M×N的矩阵表示，其中1表示疑似肿瘤的像素，0表示正常像素。由于切片厚度可以认为是一个常数，于是我们只要数连通体中1的个数就可以得到体积了。麻烦的是，可能存在多个肿瘤，这时我们只统计那些体积不小于T的。两个像素被认为是“连通的”，如果它们有一个共同的切面，如下图所示，所有6个红色的像素都与蓝色的像素连通。

输出格式：
在一行中输出肿瘤的总体积。

输入样例：
3 4 5 2
1 1 1 1
1 1 1 1
1 1 1 1
0 0 1 1
0 0 1 1
0 0 1 1
1 0 1 1
0 1 0 0
0 0 0 0
1 0 1 1
0 0 0 0
0 0 0 0
0 0 0 1
0 0 0 1
1 0 0 0
输出样例：
26

---

分析：
1.题目数据：题目数据是按照z x y给出的，这是一个长方体结构，按照案例的长方体的长是3宽是4高是5。
2.题目要求：题目要统计总的肿瘤的体积，而且肿瘤可以有多个。不难想到这题是一个三维的体连通量统计，只要在二维的基础上修改方向数组即可。

int xx[6] = {0, 0, 0, 0, 1, -1}; //上 下 左 右 前 后
int yy[6] = {0, 0, -1, 1, 0, 0};
int zz[6] = {1, -1, 0, 0, 0, 0};

3.答案要求：如果ans小于阈值，则归为0，最后统计一下即可。

---

类似的题目 请看我之前的博客
图的连通性
抠图

---

#include<iostream>
#include<queue>

using namespace std;
struct node {
    int x, y, z;
};
bool visit[61][1287][129] = {false};  // z x y
int val[61][1287][129]; // z x y
int xx[6] = {0, 0, 0, 0, 1, -1}; //上 下 左 右 前 后
int yy[6] = {0, 0, -1, 1, 0, 0};
int zz[6] = {1, -1, 0, 0, 0, 0};
int m, n, l, t; // x y z
int judge(int x, int y ,int z) { //是否越界 是否被访问 是否是1
    if(x < 1 || x > m || y < 1 || y > n || z < 1 || z > l) return 0;
    if(val[z][x][y] == 0 || visit[z][x][y] == true) return 0;
    return 1;
}
int bfs(int x,int y, int z) {
    int ans = 0;
    node tem;
    tem.x = x;
    tem.y = y;
    tem.z = z;
    queue<node >que;
    que.push(tem);
    visit[z][x][y] = true;
    while(!que.empty()) {
        node t = que.front();
        que.pop();      //每次出列 则表示一个点的完成
        ans++;
        int tx,ty,tz;
        for(int i = 0; i < 6; i ++) {
            tx = t.x + xx[i];
            ty = t.y + yy[i];
            tz = t.z + zz[i];
            if(judge(tx, ty, tz)) {
                visit[tz][tx][ty] = true;
                tem.x = tx; tem.y = ty; tem.z = tz;
                que.push(tem);
            }
        }
    }
    if(ans < t) {
        ans = 0;
    }
    return ans;
}
int main() {
    int ans = 0;
    cin >> m >> n >> l >> t;
    for (int c = 1; c <= l; c++) {
        for (int a = 1; a <= m; a++) {
            for (int b = 1; b <= n; b++) {
                cin >> val[c][a][b];
            }
        }
    }
    for (int c = 1; c <= l; c++) {
        for (int a = 1; a <= m; a ++) {
            for (int b = 1; b <= n; b ++) {
                if (val[c][a][b] == 1 && visit[c][a][b] == false) {  // 枚举起点
                    ans += bfs(a, b, c); // x y z
                }
            }
        }
    }
    cout << ans;
}