HDU2819Swap(二分图最大匹配)
题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=2819
题目大意很明确,交换图的某些行或者是某些列(可以都换),使得这个N*N的图对角线上全部都是1.
这里有一点需要说明,就是说题目的交换,其实是将原来图的某一行移到最后图的某一行,而不是指先交换两行,得到一个新图,再交换新图的两行。感觉这里比较坑。
这里先说明的一点就是,如果通过交换某些行没有办法的到解的话,那么只交换列 或者 既交换行又交换列 那也没办法得到解。其实个人感觉这个可以用矩阵的秩来解释,所有的对角线都是1,所以也就是矩阵的秩就是N,所以秩小于N就无解。另外,根据矩阵的性质,任意交换矩阵的两行 或者 两列,矩阵的秩不变,也就保证了如果通过 只交换行 或 只交换列 无法得到解的话,那么其他交换形式也必然无解。
既然说是用二分图的最大匹配,那怎么构建二分图呢,我们构建的二分图,第一部分X表示的是横坐标,第二部分Y表示纵坐标,所以范围都是1~N,然后如果a[i][j]是1,那我们就从X的i向Y的j引一条边,那么这条边的含义就可以解释为可以将Y的第j列(因为Y表示的是列的集合)移到第i列,使得a[i][i]变成1,这样就相当于是第i行第i列就变成了1,也就是说对角线多了一个1。
因此我们求这个二分图的最大匹配(目的是为了让每一列只与X中的某一行匹配),这样来就形成了N条边,那我们只需要将所有匹配的边的右边(列) 和 左边(行)所在的列 交换,这样一来对角线上这一行就成了1.
上面也也正好提示了如果最大匹配是N,那就存在解,否则无解。
1 #include <map> 2 #include <set> 3 #include <stack> 4 #include <queue> 5 #include <cmath> 6 #include <ctime> 7 #include <vector> 8 #include <cstdio> 9 #include <cctype> 10 #include <cstring> 11 #include <cstdlib> 12 #include <iostream> 13 #include <algorithm> 14 using namespace std; 15 #define eps 1e-15 16 #define MAXN 105 17 #define INF 1000000007 18 #define MAX(a,b) (a > b ? a : b) 19 #define MIN(a,b) (a < b ? a : b) 20 #define mem(a) memset(a,0,sizeof(a)) 21 22 bool G[MAXN][MAXN],vis[MAXN]; 23 int Left[MAXN],N,M,T,a[MAXN],b[MAXN]; 24 25 bool DFS(int u) 26 { 27 for(int v=0;v<=N;v++) if(G[u][v] && !vis[v]) 28 { 29 vis[v] = true; 30 if(!Left[v] || DFS(Left[v])) 31 { 32 Left[v] = u; 33 return true; 34 } 35 } 36 return false; 37 } 38 39 int main() 40 { 41 while(~scanf("%d", &N)) 42 { 43 mem(G); mem(Left); 44 int x,ans = 0; 45 for(int i=1;i<=N;i++) for(int j=1;j<=N;j++) 46 { 47 scanf("%d", &x); 48 if(x)G[i][j] = true; 49 } 50 for(int i=1;i<=N;i++)//求最大匹配 51 { 52 mem(vis); 53 if(DFS(i)) ans ++; 54 } 55 if(ans < N){printf("-1\n");continue;}//小于N无解 56 int tot = 0,j; 57 for(int i=1;i<=N;i++) 58 { 59 for(j=1;j<=N && Left[j]!=i ;j++); 60 if(i != j)//交换第i列和第j列 61 { 62 a[tot] = i; b[tot] = j; tot ++;//记录结果 63 int t = Left[i]; Left[i] = Left[j]; Left[j] = t; 64 } 65 } 66 printf("%d\n",tot); 67 for(int i=0;i<tot;i++) printf("C %d %d\n", a[i],b[i]); 68 } 69 return 0; 70 }