POJ1164 The Castle（IOI 1994）

没办法，我就这个点最兴奋，刚两点的时候在楼下看马克思，然后饿了，上来泡面，顺便寂寞一把……

 

题目大意：给出一个 n*m 矩阵，矩阵代表了一个大房间，然后矩阵的每个元素代表了该模块的信息，用 1 表示 WEST 方向是墙壁，2 表示 NORTH 方向是墙， 4 表示 EAST 方向是墙， 8 表示 SOUTH 方向是墙，题目所给矩阵的元素值就是该模块所有墙壁数值总和，现在要求统计这个大房间被分成了几个区域，其中最大区域包含几个模块。

 

显然的，BFS，DFS可解，我用的是并查集（= =！不想写搜索，对搜索的代码有阴影……）。其实这不是重点，关键是怎么把数字矩阵转化成我们想要的图形矩阵，就是从一个模块所有墙壁数值总和看出这个模块都有哪几面墙，这是这道题的亮点

 

本人用的应该是最笨的方法了，直接枚举打表所有可能……就 0 ~ 15 嘛，也不多哈……

 

看了神牛的报告，发现，用位运算的话会风骚很多……

 

1 的2进制 0001,2 的是 0010，4 的是 0100，8 的是 1000，例如 11 取2进制 1011，只有 4 和 11 取 & 运算是 0 ，即该模块在 4 方向（也就是 EAST 方向）是通的。

 

代码：

#include<stdio.h> #include<string.h> #define M 55 int father[M*M],n,m; int findx(int x) { if(father[x]!=x) father[x]=findx(father[x]); return father[x]; } void merge(int a,int b) { int x,y; x=findx(a); y=findx(b); if(x!=y)father[x]=y; //father[a]=findx(a); //father[b]=findx(b); } void fun(int a,int i,int j) { if(a==0){ merge(m*(i-1)+j,m*(i-1)+j+1); merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*(i-2)+j); merge(m*(i-1)+j,m*i+j); } else if(a==1){ merge(m*(i-1)+j,m*(i-1)+j+1); merge(m*(i-1)+j,m*(i-2)+j); merge(m*(i-1)+j,m*i+j); } else if(a==2){ merge(m*(i-1)+j,m*(i-1)+j+1); merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*i+j); } else if(a==3){ merge(m*(i-1)+j,m*(i-1)+j+1); merge(m*(i-1)+j,m*i+j); } else if(a==4){ merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*(i-2)+j); merge(m*(i-1)+j,m*i+j); } else if(a==5){ merge(m*(i-1)+j,m*(i-2)+j); merge(m*(i-1)+j,m*i+j); } else if(a==6){ merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*i+j); } else if(a==7) merge(m*(i-1)+j,m*i+j); else if(a==8){ merge(m*(i-1)+j,m*(i-1)+j+1); merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*(i-2)+j); } else if(a==9){ merge(m*(i-1)+j,m*(i-2)+j); merge(m*(i-1)+j,m*(i-1)+j+1); } else if(a==10){ merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*(i-1)+j+1); } else if(a==11) merge(m*(i-1)+j,m*(i-1)+j+1); else if(a==12){ merge(m*(i-1)+j,m*(i-1)+j-1); merge(m*(i-1)+j,m*(i-2)+j); } else if(a==13) merge(m*(i-1)+j,m*(i-2)+j); else if(a==14) merge(m*(i-1)+j,m*(i-1)+j-1); } void check() { int i,cnt,time[M*M]; memset(time,0,sizeof(time)); for(cnt=0,i=1;i<=n*m;i++) { time[father[i]]++; if(father[i]==i) cnt++; } int max=0; for(i=1;i<=n*m;i++) if(time[i]>max)max=time[i]; printf("%d/n%d/n",cnt,max); } int main() { int i,j,a; while(~scanf("%d%d",&n,&m)) { for(i=1;i<=n;i++) { for(j=1;j<=m;j++) father[m*(i-1)+j]=m*(i-1)+j; } for(i=1;i<=n;i++) { for(j=1;j<=m;j++) { scanf("%d",&a); fun(a,i,j); } } for(i=1;i<=n;i++) { for(j=1;j<=m;j++) father[m*(i-1)+j]=findx(m*(i-1)+j); } check(); } return 0; }