zoj 1654 Place the Robots

纠结了N久，不知不觉可都11点半了！！赶紧写完。。明早还有课呢。。

 

这题刚做最大二分匹配就瞅见了，因为看过黑书，觉得应该是最大独立集神马的知识，就放着没做。

 

今天下决心给它A掉，正好下的资料中有解说这题的，其实这题很像 1002 的扩大版。。。

 

给你地图，机器人可以发射上下左右的激光，如果俩机器人之间有墙隔着就可以放同一行，同一列。。。

 

真的和1002很像。。。这里面还有玻璃。。玻璃不能放机器人，但是可以穿过激光。

 

问最大放机器人的数量，就是最大独立集。

 

建图很是纠结，我觉得我建得也没错。。。但是一直WA，又换了种建图方法，依旧WA = =。。。搜题解，人家建图跟我差不多嘛。。。 = =。。

 

背疼。。。睡觉。。。呃。。

 

---

呃，昨晚想了想，中间有点过程还是有点纠结的，昨晚睡不着，想了下，终于明白了撒~

 

等于说分块后，把每个为 o 的点标记两次（行，列），如果相冲突，即一个 o 被标了两次，就建一条边，等于说是，这个点只能利用一次哈。要么第一次被标记，要么第二次标记的。呵呵。利用最大匹配的性质，每个点只能用一次。

 

 

#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <algorithm> #define MAX 55 #define MMAX 2505 using namespace std; char chmap[MAX][MAX]; int map[MAX][MAX]; int rmap[MAX][MAX]; bool g[MMAX][MMAX]; int used[MMAX],mat[MMAX]; int n,m; int nn,mm; int Augement(int x) { int i; for(i=1; i<=mm; i++) if( !used[i] && g[x][i] ) { used[i] = 1; if( mat[i] == 0 || Augement(mat[i]) ) { mat[i] = x; return 1; } } return 0; } int Hungary() { int i,sum = 0; memset(mat,0,sizeof(mat)); for(i=1; i<=nn; i++) { memset(used,0,sizeof(used)); if( Augement(i) ) sum++; } return sum; } int main() { int ncases,i,k; int mmax,cases = 1; scanf("%d",&ncases); while( ncases-- ) { scanf("%d%d",&n,&m); memset(map,0,sizeof(map)); memset(rmap,0,sizeof(map)); memset(g,0,sizeof(g)); for(i=0; i<n; i++) scanf("%s",chmap[i]); nn = mm = 0; for(i=0; i<n; i++) { k = 0; while( k < m ) if( chmap[i][k] == 'o' ) { nn++; while( k < m && chmap[i][k] != '#' ) { map[i][k] = nn; k++; } } else k++; } for(i=0; i<m; i++) { k = 0; while( k < n ) if( chmap[k][i] == 'o' ) { mm++; while( k < n && chmap[k][i] != '#' ) { rmap[k][i] = mm; k++; } } else k++; } for(i=0; i<n; i++) for(k=0; k<m; k++) if( chmap[i][k] == 'o' ) g[map[i][k]][rmap[i][k]] = 1; int ans = Hungary(); printf("Case :%d/n%d/n",cases++,ans); } return 0; }