L3-004. 肿瘤诊断

L3-004. 肿瘤诊断

时间限制

400 ms

内存限制

65536 kB

代码长度限制

8000 B

判题程序

Standard

作者

陈越

在诊断肿瘤疾病时，计算肿瘤体积是很重要的一环。给定病灶扫描切片中标注出的疑似肿瘤区域，请你计算肿瘤的体积。

输入格式：

输入第一行给出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个红色的像素都与蓝色的像素连通。

Figure 1

输出格式：

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

输入样例：

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

思路分析：

本题目可抽象为一个三维数组的遍历，使用BFS搜索即可。建立一个bool型的三维judge数组用来判别各个结点是否已经访问过。在进行BFS搜索时，可以利用增量数组来遍历与当前结点相邻的结点。

主要步骤：

1. 遍历所有结点，如果结点未被访问且值为1，BFS遍历结点，计算相邻块中的有效结点数量并返回ans。

2. 如果ans满足要讲，则将ans累加进结果

3. 输出结果

BFS思路：

1. 将当前结点压入队列q，且令结点对应的judge = true；队不空时循环执行2

2. 出队队首元素，计算加1， 利用增量数组确定相邻结点，若相邻结点的坐标合法且未访问过且值为1，入队，否则跳过

题解：

#include <stdio.h>
#include <queue>
using namespace std;

struct Node{
	int x, y, z;
	Node(int a, int b, int c){
		x = a, y = b, z = c;
	}
	Node(){}
};
int g[1286][128][60] = {0}, sum = 0;
bool judge[1286][128][60] = {false};
int addx[6] = {1, -1, 0, 0, 0, 0};
int addy[6] = {0, 0, 1, -1, 0, 0};
int addz[6] = {0, 0, 0, 0, 1, -1};

int BFS(int a, int b, int c){
	int num = 0;
	Node s = Node(a, b, c);
	queue<Node> q;
	q.push(s);
	judge[a][b][c] = true;
	while(!q.empty()){
		Node front = q.front();
		q.pop();
		num++;
		for(int i = 0; i < 6; i++){
			s.x = front.x + addx[i];
			s.y = front.y + addy[i];
			s.z = front.z + addz[i];
			if(s.x >= 0 && s.x < 1286 && s.y >= 0 && s.y < 128 && s.z >= 0 && s.z < 60){
				if(!judge[s.x][s.y][s.z] && g[s.x][s.y][s.z]){
					q.push(s);
					judge[s.x][s.y][s.z] = true;
				} 
			}
		}
	}
	return num;
}

int main(){
	int m, n, l, t;
	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", &g[j][k][i]);
			}
		}
	}
	for(int i = 0; i < l; i++){
		for(int j = 0; j < m; j++){
			for(int k = 0; k < n; k++){
				if(!judge[j][k][i] && g[j][k][i]){
					int ans = BFS(j, k, i);
					if(ans >= t) sum += ans;
				} 
			}
		}
	}
	printf("%d", sum);
	return 0;
}