博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
三维扫描 FZU 1063
阅读量:5878 次
发布时间:2019-06-19

本文共 2897 字,大约阅读时间需要 9 分钟。

Description

工业和医学上经常要用到一种诊断技术——核磁共振成像(Magnetic Resonance Imagers)。利用该技术可以对三维物体(例如大脑)进行扫描。扫描的结果用一个三维的数组来保存,数组的每一个元素表示空间的一个象素。数组的元素是0-255的整数,表示该象素的灰度。例如0表示该象素是黑色的,255表示该象素是白色的。

被扫描的物体往往是由若干个部件组合而成的。例如临床医学要对病变的器官进行检查,而器官是由一些不同的组织构成的。在实际问题中,同一个部件内部的色彩变化相对连续,而不同的部件的交界处色彩往往有突变。下面是一个简化的植物细胞的例子。

从细胞的平面图来看,该细胞大致是由四个“部件”构成的,细胞壁、细胞核、液泡和细胞质。为了方便起见,我们对部件的概念做如下的规定: 
1.如果一个象素属于某部件,则或者该象素至少与该部件的一个象素相邻,或者该象素单独组成一个部件。(说明:每一个象素与前后、左右、上下的6个象素相邻) 
2.同一个部件内部,相邻两个象素的灰度差不超过正整数M。M决定了程序识别部件的灵敏度。 
你的任务是对于给定的物体,判断该物体是由几个部件组成的。

Input

输入数据由多组数据组成。每组数据格式如下: 
第一行是三个正整数L,W,H(L,W,H≤50),表示物体的长、宽、高。 
第二行是一个整数M(0≤M≤255),表示识别部件的灵敏度。 
接下来是L×W×H个0-255的非负整数,按照空间坐标从小到大的顺序依次给出每个象素的灰度。 
说明:对于空间两点P1(x1,y1,z1)和P2(x2,y2,z2),P1

Output

对于每组数据,输出仅一行包含一个整数M,即一共识别出几个部件。

Sample Input

2 2 2 0 1 1 1 1 2 2 2 2

Sample Output

2
 
 
 
#include 
#include
#include
using namespace std;int a[52][52][52][2];void dns(int i,int j,int k,int m){ if(a[i][j][k][1]>0&&a[i][j][k+1][1]<0&&a[i][j][k+1][0]>=0&&abs(a[i][j][k][0]-a[i][j][k+1][0])<=m) { a[i][j][k+1][1]=a[i][j][k][1]; dns(i,j,k+1,m); } if(a[i][j][k][1]>0&&a[i][j][k-1][1]<0&&a[i][j][k-1][0]>=0&&abs(a[i][j][k][0]-a[i][j][k-1][0])<=m) { a[i][j][k-1][1]=a[i][j][k][1]; dns(i,j,k-1,m); } if(a[i][j][k][1]>0&&a[i][j-1][k][1]<0&&a[i][j-1][k][0]>=0&&abs(a[i][j][k][0]-a[i][j-1][k][0])<=m) { a[i][j-1][k][1]=a[i][j][k][1]; dns(i,j-1,k,m); } if(a[i][j][k][1]>0&&a[i][j+1][k][1]<0&&a[i][j+1][k][0]>=0&&abs(a[i][j][k][0]-a[i][j+1][k][0])<=m) { a[i][j+1][k][1]=a[i][j][k][1]; dns(i,j+1,k,m); } if(a[i][j][k][1]>0&&a[i-1][j][k][1]<0&&a[i-1][j][k][0]>=0&&abs(a[i][j][k][0]-a[i-1][j][k][0])<=m) { a[i-1][j][k][1]=a[i][j][k][1]; dns(i-1,j,k,m); } if(a[i][j][k][1]>0&&a[i+1][j][k][1]<0&&a[i+1][j][k][0]>=0&&abs(a[i][j][k][0]-a[i+1][j][k][0])<=m) { a[i+1][j][k][1]=a[i][j][k][1]; dns(i+1,j,k,m); }}int main(){ int l,w,h,m; while((scanf("%d%d%d%d",&l,&w,&h,&m))!=EOF) { memset(a,-1,sizeof(a)); for(int i=1;i<=l;i++) { for(int j=1;j<=w;j++) { for(int k=1;k<=h;k++) { scanf("%d",&a[i][j][k][0]); } } } int num=0; for(int i=1;i<=l;i++) { for(int j=1;j<=w;j++) { for(int k=1;k<=h;k++) { if(a[i][j][k][1]<0&&a[i][j][k][0]>=0) { num++; a[i][j][k][1]=num; dns(i,j,k,m); } } } } cout<
<

 

这道题我用的DNS(深度优先搜索 ) 
 

转载于:https://www.cnblogs.com/guofeng1022/p/4213049.html

你可能感兴趣的文章
长春理工大学第十四届程序设计竞赛(重现赛)I.Fate Grand Order
查看>>
好作品地址
查看>>
[翻译]Protocol Buffer 基础: C++
查看>>
runloop与线程的关系
查看>>
[Bzoj2246]迷宫探险(概率+DP)
查看>>
[译] 感受 4px 基线网格带来的便利
查看>>
oracle常用函数
查看>>
MYBATIS
查看>>
详解消息队列的设计与使用
查看>>
iOS 项目优化
查看>>
筛选出sql 查询结果中 不包含某个字符
查看>>
8进制与16进制
查看>>
使用Sqoop从mysql向hdfs或者hive导入数据时出现的一些错误
查看>>
mybatis:Invalid bound statement (not found)
查看>>
电脑中毒的现象
查看>>
django表单操作之django.forms
查看>>
ZipOutputStream出现多层目录问题
查看>>
webSocket vnc rfb
查看>>
列表推导式 生成器表达式
查看>>
控制子窗口的高度
查看>>