PAT甲级——1091 Acute Stroke (广度优先搜索BFS)

最新推荐文章于 2021-02-04 18:45:13 发布

weixin_30849591

最新推荐文章于 2021-02-04 18:45:13 发布

阅读量62

点赞数

CC 4.0 BY-SA版权

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/yinhao-ing/p/11107933.html

本文介绍了一种用于识别急性脑卒中体积的算法。通过分析MRI切片中核心区域的结果，该算法能够计算出卒中核心的总体积。输入包括切片数量、大小和阈值，输出为卒中核心体积。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

本文同步发布在优快云：https://blog.youkuaiyun.com/weixin_44385565/article/details/94207638

1091 Acute Stroke (30 分)

One important factor to identify acute stroke (急性脑卒中) is the volume of the stroke core. Given the results of image analysis in which the core regions are identified in each MRI slice, your job is to calculate the volume of the stroke core.

Input Specification:

Each input file contains one test case. For each case, the first line contains 4 positive integers:

Then

Figure 1

Output Specification:

For each case, output in a line the total volume of the stroke core.

Sample Input:

Sample Output:

我一开始没理解题目，看完别人的代码才明白，题目本身不难，比较坑的就是题意理解了~

题目大意：在一堆大脑切片里寻找有问题的块里面的值为1的节点数，一共有L片切片，每个切片都由MxN的像素点矩阵组成，形成一个长M宽N高L的长方体。每个像素点的值为0或1, 1就是有问题的像素点，每个节点的邻接点为上下左右前后这6个点。有问题的块为相互连接的值为1且数目超过T的像素点的集合，将有问题的块里面的值为1的像素点的数目累加即为答案。

思路：遍历三维数组，对每个值为1的像素点进行BFS，并且标记已经访问过的节点。

 1 #include <iostream>
 2 #include <queue>
 3 using namespace std;
 4 
 5 struct node {
 6     int x, y, z;
 7 };
 8 
 9 int M, N, L, T,
10     pixel[61][1287][129],
11     X[6] = { -1,1,0,0,0,0 },
12     Y[6] = { 0,0,-1,1,0,0 },
13     Z[6] = { 0,0,0,0,-1,1 };
14 bool visit[61][1287][129] = { false };
15 
16 bool check(int& x, int& y, int& z);
17 int count(const int& x, const int& y, const int& z);
18 
19 int main()
20 {
21     int total = 0;
22     scanf("%d%d%d%d", &M, &N, &L, &T);
23     for (int i = 0; i < L; i++)
24         for (int j = 0; j < M; j++)
25             for (int k = 0; k < N; k++)
26                 scanf("%d", &pixel[i][j][k]);
27     for (int i = 0; i < L; i++) {
28         for (int j = 0; j < M; j++) {
29             for (int k = 0; k < N; k++) {
30                 if (pixel[i][j][k] == 1 && visit[i][j][k] == false) {
31                     total += count(i, j, k);
32                 }
33             }
34         }
35     }
36     printf("%d\n", total);
37     return 0;
38 }
39 
40 int count(const int& x, const int& y, const int& z) {
41     int n = 0;
42     node dot = { x,y,z };
43     queue<node> Q;
44     Q.push(dot);
45     visit[x][y][z] = true;
46     while (!Q.empty()) {
47         dot = Q.front();
48         Q.pop();
49         n++;
50         for (int i = 0; i < 6; i++) {
51             node tmp;
52             tmp.x = dot.x + X[i];
53             tmp.y = dot.y + Y[i];
54             tmp.z = dot.z + Z[i];
55             if (check(tmp.x, tmp.y, tmp.z)) {
56                 Q.push(tmp);
57                 visit[tmp.x][tmp.y][tmp.z] = true;
58             }
59         }
60     }
61     if (n >= T)
62         return n;
63     else
64         return 0;
65 }
66 
67 bool check(int& x, int& y, int& z) {
68     if (x < 0 || x >= L || y < 0 || y >= M || z < 0 || z >= N || 
69         pixel[x][y][z] == 0 || visit[x][y][z] == true) {
70         return false;
71     }
72     return true;
73 }