ZJU 1002 Fire Net - 二分图最大匹配

题目大意：

一个N*N的网格中(N<=4)，有些格子是墙。问在这个网格中最多能放置多少个士兵，使得他们在四个方向上不能相互攻击。

分析：

这是一个古老的题目了，我最开始学ACM的时候一眼就能看出来这是一个搜索题，并且也能毫不费力地写出一段搜索代码把它过掉。

后来听罗老师讲了这个题之后才知道，原来还有更高效的匹配算法来解决这个问题。（按照罗老师的原话说，这个题目出题的本意就是让大家用匹配来做……>_<）

可以这样来构造匹配模型。假设没有墙阻挡，那么把地图的行和列作为二部图的两侧顶点，i行和j列有交叉，那么二部图中g[i][j]=1。现在有墙隔开，如果墙把第i行分割成了2段，那么二部图左边节点i就拆成两个点i和i'，若i'和j有交叉，那么二部图中g[i'][j]=1。列的拆分同理。

这样得到的二部图，求出一个最大匹配即得到了最多的放置士兵数。

代码也容易写出，感觉自己写得不是很精炼。

 

1. /*
2. ZJU1002 Pro
3. */
5. #include <stdio.h>
6. #include <memory.h>
8. #define clr(a) memset(a,0,sizeof(a))
9. #define N 8
11. int find(int i,int m,int g[][N],int mat[],int tmat[]){
12.     int v,j;
13.     for(j=0;j<m;j++) if(g[i][j]&&tmat[j]==0){
14.         tmat[j]=1;v=mat[j];mat[j]=i;
15.         if(v==-1||find(v,m,g,mat,tmat)) return 1;
16.         mat[j]=v;
17.     }
18.     return 0;
19. }
20. int match(int g[][N],int n,int m){
21.     int i,k=0,tmat[N],mat[N];
22.     for(i=0;i<m;i++) mat[i]=-1;
23.     for(i=0;i<n;i++){clr(tmat); k+=find(i,m,g,mat,tmat);}
24.     return k;
25. }
27. int main()
28. {
29.     int i,j,k,w,n,m;
30.     char map[N][N];
31.     char a[N][N],fa;
32.     char b[N][N],fb;
33.     int g[N][N];
34.     
35.     while(scanf("%d",&w),w){
36.         //init
37.         clr(g); n=m=0;
38.         //input
39.         for(i=0;i<w;i++){
40.             scanf("%s",map[i]);
41.         }
42.         //create a b map
43.         for(i=0;i<w;i++){
44.             fa=0; fb=0;
45.             for(j=0;j<w;j++){
46.                 if(map[i][j]=='.'){
47.                     if(fa==0) n++;
48.                     fa=1; a[i][j]=n;
49.                 }
50.                 else fa=0;
51.                 if(map[j][i]=='.'){
52.                     if(fb==0) m++;
53.                     fb=1; b[j][i]=m;
54.                 }
55.                 else fb=0;
56.             }
57.         }
58.         //create g
59.         for(i=0;i<w;i++) for(j=0;j<w;j++){
60.             if(map[i][j]=='.'){
61.                 g[a[i][j]-1][b[i][j]-1]=1;
62.             }
63.         }
64.         //matching
65.         k=match(g,n,m);
66.         printf("%d/n",k);
67.     }
68.     return 0;
69. }