在诊断肿瘤疾病时,计算肿瘤体积是很重要的一环。给定病灶扫描切片中标注出的疑似肿瘤区域,请你计算肿瘤的体积。
输入格式:
输入第一行给出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
测试点0 | 同样例 |
测试点1 | 只有一层切片(L=1) |
测试点2 | 无连通体积超过阈值的肿瘤块,输出为0 |
测试点3 | 有且只有一条环状肿瘤;阈值T为1 |
测试点4 | 最大数据,dfs算法会发生段错误 |
测试点5 | 最大数据,dfs算法会发送段错误 |
dfs算法(25分)
虽然代码简单但是最后两个测试点会因为递归层数过多发生段错误
#include<bits/stdc++.h>
using namespace std;
int n,m,l,coord[130][1288][62];
vector<vector<int>>cube={{1,0,0},{-1,0,0},{0,1,0},{0,-1,0},{0,0,1},{0,0,-1}};
int dfs(int a,int b,int c){
if(coord[a][b][c]!=1)return 0;
coord[a][b][c]=0;
int sum=1;
for(auto j:cube)sum+=dfs(a+j[0],b+j[1],c+j[2]);
return sum;
}
int main()
{
int T,p,aa=0;
cin>>m>>n>>l>>T;
for(int z=1;z<=l;z++){
for(int y=1;y<=m;y++){
for(int x=1;x<=n;x++)cin>>coord[x][y][z];
}
}
for(int z=1;z<=l;z++){
for(int y=1;y<=m;y++){
for(int x=1;x<=n;x++){
int k=dfs(x,y,z);
if(k>=T)aa+=k;
}
}
}
cout<<aa;
}
bfs算法(30分)
#include<bits/stdc++.h>
using namespace std;
int n,m,l,coord[130][1288][62];
vector<vector<int>>cube={{1,0,0},{-1,0,0},{0,1,0},{0,-1,0},{0,0,1},{0,0,-1}};
//6个方位
int bfs(int a,int b,int c){//统计相连的肿瘤块数
if(!coord[a][b][c])return 0;
queue<vector<int>>q;
int sum=1;
coord[a][b][c]=0;//清除已被统计过的肿瘤
q.push({a,b,c});
while(!q.empty()){//直到相连的肿瘤被统计完
vector<int>s=q.front();
q.pop();
for(auto t:cube){
int x=s[0]+t[0],y=s[1]+t[1],z=s[2]+t[2];
if(coord[x][y][z]){
coord[x][y][z]=0;//清除已被统计过的肿瘤
sum++;
q.push({x,y,z});
}
}
}
return sum;
}
int main()
{
int T,p,aa=0;
cin>>m>>n>>l>>T;
for(int z=1;z<=l;z++){
for(int y=1;y<=m;y++){
for(int x=1;x<=n;x++)cin>>coord[x][y][z];//建立空间直角坐标系
}
}
for(int z=1;z<=l;z++){
for(int y=1;y<=m;y++){
for(int x=1;x<=n;x++){
int k=bfs(x,y,z);
if(k>=T)aa+=k;//只计连通体积不小于阈值的肿瘤
}
}
}
cout<<aa;
}