时间限制:2秒
空间限制:65536K
研究地球空间科学的永强想研究海岸线的长度和海岸线面积之间的关系,为此他找来了很多航拍图像。在航拍图像上利用图像分割的方法,把图像的每个像素标记成陆地(1)和水面(0)。
示例图片:
现在永强想知道每张图中陆地部分的面积。
已知每张图最底部的一条边都是陆地,并且在一张图上陆地都是四邻域联通的。
但是永强发现分割的结果有很多的噪声,于是他定义了如下规则试图去除噪声:
a) 如果一个水面区域被陆地包围,则将这个区域记为陆地;
b) 在a的基础上如果一个陆地区域不和底边的陆地相连,那么这是一个岛屿,不计入陆地的面积。
输入描述:
第一行为两个整数m和n, 接下来m行每行会有n个数,0表示水面,1表示陆地。
输出描述:
去噪后陆地部分的面积。
输入例子1:
5 6 1 0 0 1 0 0 0 0 1 0 1 0 1 1 1 1 1 0 1 0 0 1 1 1 1 1 1 1 1 1
输出例子1:
21
#include<iostream>
using namespace std;
int m,n;
int *p;
int sum;
/**让一个值为key的节点和他周围值为key的节点连通
*key!=set
*/
void link(int key,int sx,int sy,int set)
{
sum++;
p[sy*n+sx]=set;
if(sx-1>=0&&p[sy*n+sx-1]==key) link(key,sx-1,sy,set);
if(sx+1<n&&p[sy*n+sx+1]==key) link(key,sx+1,sy,set);
if(sy-1>=0&&p[(sy-1)*n+sx]==key) link(key,sx,sy-1,set);
if(sy+1<m&&p[(sy+1)*n+sx]==key) link(key,sx,sy+1,set);
}
int main()
{
int i,j;
cin>>m>>n;
p=new int[m*n];
for(i=0;i<m*n;i++)
{
cin>>p[i];
}
/*连通所有海洋且标记为-1*/
for(i=0;i<n;i++)
{
if(p[i]==0) link(0,i,0,-1);
}
for(i=0;i<m;i++)
{
if(p[i*n]==0) link(0,0,i,-1);
if(p[i*n+n-1]==0) link(0,n-1,i,-1);
}
/*其余的都为陆地,将地中海变为陆地*/
for(i=1;i<m-1;i++)
{
for(j=1;j<n-1;j++)
{
if(p[i*n+j]==0) p[i*n+j]=1;
}
}
sum=0;
/*计算陆地面积*/
link(1,0,m-1,2);
cout<<sum<<endl;
}