hdu1507最大匹配
题意:给出n*m的矩阵方块,某些方块是ponds,这些方块不能出售,然后能出售的是1*2的方块
要求求出有多少块1*2的方块能够出售,并输出这些方块的坐标(specia judge)
关键是怎么建出二分图,然后找最大匹配
可以把i+j==奇数的点提取出来,成为X集合,那么与(i,j)相邻的点就是Y集合了,然后就是二分图最大匹配
1 //关键在于建图 2 #include <stdio.h> 3 #include <string.h> 4 #include <vector> 5 using namespace std; 6 const int N = 10000+10; 7 int Map[111][111]; 8 int n,m,k; 9 vector<int> G[N]; 10 int dx[] = {0,0,-1,1}; 11 int dy[] = {-1,1,0,0}; 12 bool vis[N]; 13 int cy[N]; 14 int cx[N]; 15 bool dfs(int u) 16 { 17 for(int i=0; i<G[u].size(); ++i) 18 { 19 int v = G[u][i]; 20 if(!vis[v]) 21 { 22 vis[v] = true; 23 if(cy[v] == -1 || dfs(cy[v])) 24 { 25 cy[v] = u; 26 cx[u] = v; 27 return true; 28 } 29 } 30 } 31 return false; 32 } 33 int MaxMatch() 34 { 35 memset(cy, -1, sizeof(cy)); 36 memset(cx, -1, sizeof(cx)); 37 int t = m * n; 38 int cnt = 0; 39 for(int i=0; i<t; ++i) 40 { 41 if(cx[i] == -1) 42 { 43 memset(vis,0,sizeof(vis)); 44 cnt += dfs(i); 45 } 46 } 47 return cnt; 48 } 49 int main() 50 { 51 int i,j,x,y; 52 while(true) 53 { 54 scanf("%d%d",&n,&m); 55 if(n==0 && m==0) 56 break; 57 scanf("%d",&k); 58 for(i=0; i<N; ++i) 59 G[i].clear(); 60 memset(Map, 0, sizeof(Map)); 61 for(i=0; i<k; ++i) 62 { 63 scanf("%d%d",&x,&y); 64 x-=1; y-=1; 65 Map[x][y] = 1; 66 } 67 for(i=0; i<n; ++i) 68 for(j=0; j<m; ++j) 69 if(Map[i][j] == 0 && (i+j)&1==1) 70 { 71 x = i * n + j; 72 for(k=0; k<4; ++k) 73 { 74 int a = dx[k] + i; 75 int b = dy[k] + j; 76 if(a>=0 && b>=0 && a<n && b<m && Map[a][b] == 0) 77 { 78 y = a * n + b; 79 G[x].push_back(y); 80 } 81 } 82 } 83 int ans = MaxMatch(); 84 printf("%d\n",ans); 85 int t = n * m; 86 for(i=0; i<t; ++i) 87 if(cx[i] != -1) 88 printf("(%d,%d)--(%d,%d)\n",i/n+1,i%n+1,cx[i]/n+1,cx[i]%n+1); 89 90 } 91 return 0; 92 }