PAT_A1091#Acute Stroke

Source:

PAT A1091 Acute Stroke (30 分)

Description:

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, Land 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 (≤) 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

Keys:

• 广度优先搜索（Breadth First Search）

Attention:

• 4,5测试点用递归会爆栈（段错误）

Code:

 /*
 time: 2019-08-26 17:10:21
 problem: PAT_A1091#Acute Stroke
 AC: 36:38

 题目大意：
 给出脑部各层的切片，患病几率与连续病变区域的数量有关，计算患病区域数量
 输入：
 第一行给出，各层的长度M<=1286和宽度N<=128,层数L<=60,患病阈值T（连续病变区域>=T时患病）
 接下来L行，给出各层的信息，1病变，0正常
 输出：
 打印患病的病变区域数量

 基本思路：
 BFS统计各个块的有效结点数
 */
 #include<cstdio>
 #include<queue>
 using namespace std;
 const int L=,M=,N=;
 struct node
 {
     int x,y,z;
 };
 int l,m,n,t,cnt,sum=,grap[L][M][N],vis[L][M][N];
 int X[]={,-,,,,},Y[]={,,,-,,},Z[]={,,,,,-};

 void BFS(int z,int x, int y)
 {
     grap[z][x][y]=;
     queue<node> q;
     q.push(node{x,y,z});
     while(!q.empty())
     {
         node r = q.front();
         q.pop();cnt++;
         for(int i=; i<; i++)
         {
             int  u=r.x+X[i],v=r.y+Y[i],w=r.z+Z[i];
             if(u>=&&u<m&&v>=&&v<n&&w>=&&w<l)
                 if(grap[w][u][v]==)
                 {
                     grap[w][u][v]=;
                     q.push(node{u,v,w});
                 }
         }
     }
 }

 int main()
 {
 #ifdef ONLINE_JUDGE
 #else
     freopen("Test.txt", "r", stdin);
 #endif // ONLINE_JUDGE

     scanf("%d%d%d%d", &m,&n,&l,&t);
     for(int k=; k<l; k++)
         for(int i=; i<m; i++)
             for(int j=; j<n; j++)
                 scanf("%d", &grap[k][i][j]);
     for(int k=; k<l; k++)
         for(int i=; i<m; i++)
             for(int j=; j<n; j++)
                 if(grap[k][i][j]==)
                 {
                     cnt=;
                     BFS(k,i,j);
                     if(cnt>=t)
                         sum+=cnt;
                 }
     printf("%d", sum);

     return ;
 }