USACO 2.1.1 Castle

Translation

我们憨厚的USACO主人公农夫约翰以无法想象的运气,在他生日那天收到了一份特别的礼物：一张“幸运爱尔兰”（一种彩票）。结果这张彩票让他获得了这次比赛唯一的奖品——坐落于爱尔兰郊外的一座梦幻般的城堡！

喜欢吹嘘的农夫约翰立刻回到有着吹嘘传统的威斯康辛老家开始吹嘘了, 农夫约翰想要告诉他的奶牛们关于他城堡的一切。他需要做一些吹嘘前的准备工作：比如说知道城堡有多少个房间，每个房间有多大。另外，农夫约翰想要把一面单独的墙（指两个单位间的墙）拆掉以形成一个更大的房间。 你的工作就是帮农夫约翰做以上的准备，算出房间数与房间的大小。

城堡的平面图被划分成M*N(1 <=M,N<=50）个正方形的单位，一个这样的单位可以有0到4面墙环绕。城堡周围一定有外墙环绕以遮风挡雨。（就是说平面图的四周一定是墙。）

请仔细研究下面这个有注解的城堡平面图：

     1   2   3   4   5   6   7
   #############################
 1 #   |   #   |   #   |   |   #
   #####---#####---#---#####---#   
 2 #   #   |   #   #   #   #   #
   #---#####---#####---#####---#
 3 #   |   |   #   #   #   #   #   
   #---#########---#####---#---#
 4 # ->#   |   |   |   |   #   #   
   ############################# 


 # =墙壁    -,| = 没有墙壁
 -> =指向一面墙，这面墙推掉的话我们就有一间最大的新房间

友情提示，这个城堡的平面图是7×4个单位的。一个“房间”的是平面图中一个由“#”、“-”、“|”围成的格子（就是图里面的那一个个的格子）。比如说这个样例就有5个房间。（大小分别为9、7、3、1、8个单位）

移去箭头所指的那面墙，可以使2个房间合为一个新房间，且比移去其他墙所形成的房间都大。

城堡保证至少有2个房间，而且一定有一面墙可以被移走。

PROGRAM NAME: castle

INPUT FORMAT: 第一行有两个整数：M和N 城堡的平面图用一个由数字组成的矩阵表示，一个数字表示一个单位，矩阵有N行M列。输入与样例的图一致。

每一个单位的数字告诉我们这个单位的东西南北是否有墙存在。每个数字是由以下四个整数的某个或某几个加起来的(四面都没有墙的话，这个数字应该为0)。

1: 在西面有墙
2: 在北面有墙
4: 在东面有墙
8: 在南面有墙

城堡内部的墙会被规定两次。比如说（1，1）南面的墙，亦会被标记为（2，1）北面的墙。

OUTPUT FORMAT:(file castle.out)

输出包含如下4行:

第 1 行: 城堡的房间数目。

第 2 行: 最大的房间的大小

第 3 行: 移除一面墙能得到的最大的房间的大小

第 4 行: 移除哪面墙可以得到面积最大的新房间。

选择最佳的墙来推倒。有多解时选最靠西的，仍然有多解时选最靠南的。同一格子北边的墙比东边的墙更优先。

用该墙的南邻单位的北墙或西邻单位的东墙来表示这面墙，方法是输出邻近单位的行数、列数和墙的方位（"N"（北）或者"E"（东））。

Solution

PS：好久没写代码了，连广搜都改了好久23333

广搜裸题，先把数字矩阵改造成邻接矩阵，然后广搜得到第一二个问题的答案；随后枚举每一面有效的墙，广搜计算墙被移除后的该房间面积即可。注意枚举墙的顺序

Code

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=50+10;
int map[maxn*2][maxn*2],p[maxn*2][maxn*2],q[maxn*maxn*40][2];
int dir[4][2]={{0,-1},{-1,0},{0,1},{1,0}},mod[4]={1,2,4,8};
int main(){
	int i,j,k,m,n,max_area=0,room_num=0,wall_x,wall_y;
	freopen("castle.in","r",stdin);
	freopen("castle.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(i=2;i<=2*m;i+=2)
		for(j=2;j<=2*n;j+=2){
			scanf("%d",&map[i][j]);
			map[i][j]=15-map[i][j];
			for(k=0;k<=3;k++)
				if(map[i][j]%(2*mod[k])==mod[k]){
					map[i][j]-=mod[k];
					map[i+dir[k][0]][j+dir[k][1]]=2;
				}
			map[i][j]=1;
		}
	for(i=2;i<=2*m;i+=2)
		for(j=2;j<=2*n;j+=2){
			int f=0,l=1;
			int area=0;
			if(!p[i][j]){
				p[i][j]=1;
				q[1][0]=i;q[1][1]=j;
				area++;
				while(f<l){
					f++;
					int x=q[f][0],y=q[f][1];
					for(k=0;k<=3;k++)
						if(map[x+dir[k][0]][y+dir[k][1]]==2 && !p[x+2*dir[k][0]][y+2*dir[k][1]]){
							q[++l][0]=x+2*dir[k][0];
							q[l][1]=y+2*dir[k][1];
							p[x+2*dir[k][0]][y+2*dir[k][1]]=1;
							area++;
						}
				}
				if(max_area<area)max_area=area;
				room_num++;
			}
		}
	printf("%d\n%d\n",room_num,max_area);
	for(j=2;j<=2*n;j++)
		for(i=2*m;i>=2;i--)
			if(!map[i][j] && (i+j)%2){
				memset(p,0,sizeof(p));
				map[i][j]=2;
				int x=i,y=j;
				for(k=0;k<=3;k++)
					if(map[x+dir[k][0]][y+dir[k][1]]){
						x+=dir[k][0];
						y+=dir[k][1];
						break;
					}
				int f=0,l=1;
				int area=0;
				if(!p[x][y]){
					p[x][y]=1;
					q[1][0]=x;q[1][1]=y;
					area++;
					while(f<l){
						f++;
						x=q[f][0],y=q[f][1];
						for(k=0;k<=3;k++)
							if(map[x+dir[k][0]][y+dir[k][1]]==2 && !p[x+2*dir[k][0]][y+2*dir[k][1]] && map[x+2*dir[k][0]][y+2*dir[k][1]]==1){
								q[++l][0]=x+2*dir[k][0];
								q[l][1]=y+2*dir[k][1];
								p[x+2*dir[k][0]][y+2*dir[k][1]]=1;
								area++;
							}
					}
					if(max_area<area){
						max_area=area;
						wall_x=i;wall_y=j;
					}
					map[i][j]=0;
				}
			}
	printf("%d\n",max_area);
	if(map[wall_x+1][wall_y]==1)printf("%d %d N\n",(wall_x+1)/2,wall_y/2);
	else if(map[wall_x][wall_y-1]==1)printf("%d %d E\n",wall_x/2,(wall_y-1)/2);
	return 0;
}