PAT 1091 Acute Stroke [难][bfs]

本文介绍了一种基于MRI图像分析的急性脑卒中核心体积计算方法,通过三维矩阵存储和BFS遍历算法,准确识别并计算超过阈值的中风区域体积。

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

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: M, N, L and T, where M and N are the sizes of each slice (i.e. pixels of a slice are in an M×N matrix, and the maximum resolution is 1286 by 128); L (60) is the number of slices of a brain; and T is the integer threshold (i.e. if the volume of a connected core is less than T, then that core must not be counted).

Then L slices are given. Each slice is represented by an M×N matrix of 0's and 1's, where 1 represents a pixel of stroke, and 0 means normal. Since the thickness of a slice is a constant, we only have to count the number of 1's to obtain the volume. However, there might be several separated core regions in a brain, and only those with their volumes no less than T are counted. Two pixels are connected and hence belong to the same region if they share a common side, as shown by Figure 1 where all the 6 red pixels are connected to the blue one.

 

Figure 1

Output Specification:

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

Sample Input:

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

Sample Output:

26

 题目大意:先给出4个数,M行,N列,L共有几块,T阈值。1表示中风,0表正常,如果两块有相同的边,那么它们就是相连的。肿瘤块大小大于T才是肿瘤。

 //这个真看不太懂,以前看过的,都不理解,什么意思,什么是块连在一起了呢?真不理解啊。

 代码转自:https://www.liuchuo.net/archives/2307

#include <cstdio>
#include <queue>
using namespace std;
struct node {
    int x, y, z;
};
int m, n, l, t;
int X[6] = {1, 0, 0, -1, 0, 0};//方向依次是:同一个平面上:下,右,往下一个平面正下,
int Y[6] = {0, 1, 0, 0, -1, 0};//同一个平面,上,左,往上一个平面上
int Z[6] = {0, 0, 1, 0, 0, -1};
int arr[1300][130][80];
bool visit[1300][130][80];
bool judge(int x, int y, int z) {
    if(x < 0 || x >= m || y < 0 || y >= n || z < 0 || z >= l) return false;
    if(arr[x][y][z] == 0 || visit[x][y][z] == true) return false;
    return true;
}
int bfs(int x, int y, int z) {
    int cnt = 0;
    node temp;
    temp.x = x, temp.y = y, temp.z = z;
    queue<node> q;
    q.push(temp);
    visit[x][y][z] = true;
    while(!q.empty()) {
        node top = q.front();
        q.pop();
        cnt++;//这里对块进行计数。
        for(int i = 0; i < 6; i++) {
            int tx = top.x + X[i];//6个方向进行判断。
            int ty = top.y + Y[i];
            int tz = top.z + Z[i];
            if(judge(tx, ty, tz)) {
        //这里需要判断下标是否是合法的,并且数组是否是1或者未被访问过。
                visit[tx][ty][tz] = true;
                temp.x = tx, temp.y = ty, temp.z = tz;
                q.push(temp);
            }
        }
    }
    if(cnt >= t)
        return cnt;
    else
        return 0;//否则返回0.

}
int main() {
    scanf("%d %d %d %d", &m, &n, &l, &t);
    for(int i = 0; i < l; i++)
        for(int j = 0; j < m; j++)
            for(int k = 0; k < n; k++)//使用三维数组存储。
                scanf("%d", &arr[j][k][i]);//先读入进来。
    int ans = 0;
    for(int i = 0; i < l; i++) {
        for(int j = 0; j < m; j++) {
            for(int k = 0; k < n; k++) {//对每一个点都开始遍历,
                if(arr[j][k][i] == 1 && visit[j][k][i] == false)
                    ans += bfs(j, k, i);//i=0时,也就是从第0块开始。
            }
        }
    }
    printf("%d", ans);
    return 0;
}

 

//代码是看懂了,但是有点疑问,就是如果l不是5呢?为什么它们就有公共的边呢?这里也没有进行判断啊,还不不明白。

1.使用三维数据进行输入,并且标记是否被访问过,使用bfs进行遍历。

2.对六个方向都进行了判断,并且使用了一个judge函数判断下标是否越界,并且是否可以访问该点,十分厉害,值得学习。

3.使用bfs也就是需要用到队列,判断条件是队列是否为空,学习了。

转载于:https://www.cnblogs.com/BlueBlueSea/p/9693094.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值