[HDOJ2819]Swap(二分图最大匹配, 匈牙利算法)
                                                            生活随笔
收集整理的這篇文章主要介紹了
                                [HDOJ2819]Swap(二分图最大匹配, 匈牙利算法)
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.                        
                                題目鏈接:http://acm.split.hdu.edu.cn/showproblem.php?pid=2819
題意:給一張n*n的01矩陣,可以任意交換其中的行或者列,問是否可以交換出來一個(gè)對角線上都是1的矩陣。
按行列號建圖,如果(i,j)為1的話,則i和j就有一條邊。匹配出的結(jié)果可以認(rèn)為如何交換使得行列相等,輸出結(jié)果要注意不要忘記更新linker,以保證在正確的位置。
1 /* 2 ━━━━━┒ギリギリ♂ eye! 3 ┓┏┓┏┓┃キリキリ♂ mind! 4 ┛┗┛┗┛┃\○/ 5 ┓┏┓┏┓┃ / 6 ┛┗┛┗┛┃ノ) 7 ┓┏┓┏┓┃ 8 ┛┗┛┗┛┃ 9 ┓┏┓┏┓┃ 10 ┛┗┛┗┛┃ 11 ┓┏┓┏┓┃ 12 ┛┗┛┗┛┃ 13 ┓┏┓┏┓┃ 14 ┃┃┃┃┃┃ 15 ┻┻┻┻┻┻ 16 */ 17 #include <algorithm> 18 #include <iostream> 19 #include <iomanip> 20 #include <cstring> 21 #include <climits> 22 #include <complex> 23 #include <fstream> 24 #include <cassert> 25 #include <cstdio> 26 #include <bitset> 27 #include <vector> 28 #include <deque> 29 #include <queue> 30 #include <stack> 31 #include <ctime> 32 #include <set> 33 #include <map> 34 #include <cmath> 35 using namespace std; 36 #define fr first 37 #define sc second 38 #define cl clear 39 #define BUG puts("here!!!") 40 #define W(a) while(a--) 41 #define pb(a) push_back(a) 42 #define Rint(a) scanf("%d", &a) 43 #define Rll(a) scanf("%I64d", &a) 44 #define Rs(a) scanf("%s", a) 45 #define Cin(a) cin >> a 46 #define FRead() freopen("in", "r", stdin) 47 #define FWrite() freopen("out", "w", stdout) 48 #define Rep(i, len) for(int i = 0; i < (len); i++) 49 #define For(i, a, len) for(int i = (a); i < (len); i++) 50 #define Cls(a) memset((a), 0, sizeof(a)) 51 #define Clr(a, x) memset((a), (x), sizeof(a)) 52 #define Full(a) memset((a), 0x7f7f7f, sizeof(a)) 53 #define lrt rt << 1 54 #define rrt rt << 1 | 1 55 #define pi 3.14159265359 56 #define RT return 57 #define lowbit(x) x & (-x) 58 #define onecnt(x) __builtin_popcount(x) 59 typedef long long LL; 60 typedef long double LD; 61 typedef unsigned long long ULL; 62 typedef pair<int, int> pii; 63 typedef pair<string, int> psi; 64 typedef pair<LL, LL> pll; 65 typedef map<string, int> msi; 66 typedef vector<int> vi; 67 typedef vector<LL> vl; 68 typedef vector<vl> vvl; 69 typedef vector<bool> vb; 70 71 const int maxn = 110; 72 int nu, nv; 73 int G[maxn][maxn]; 74 int linker[maxn]; 75 bool vis[maxn]; 76 77 bool dfs(int u) { 78 For(v, 1, nv+1) { 79 if(G[u][v] && !vis[v]) { 80 vis[v] = 1; 81 if(linker[v] == -1 || dfs(linker[v])) { 82 linker[v] = u; 83 return 1; 84 } 85 } 86 } 87 return 0; 88 } 89 90 int hungary() { 91 int ret = 0; 92 Clr(linker, -1); 93 For(u, 1, nu+1) { 94 Cls(vis); 95 if(dfs(u)) ret++; 96 } 97 return ret; 98 } 99 100 int n, m, k; 101 vector<pii> ans; 102 103 int main() { 104 // FRead(); 105 while(~Rint(n)) { 106 Cls(vis); Cls(G); ans.clear(); 107 nu = nv = n; 108 For(i, 1, n+1) { 109 For(j, 1, n+1) { 110 Rint(G[i][j]); 111 } 112 } 113 int ret = hungary(); 114 if(ret != n) { 115 puts("-1"); 116 continue; 117 } 118 For(i, 1, nu+1) { 119 int j; 120 for(j = i; j <= nu; j++) { 121 if(linker[j] == i) break; 122 } 123 if(i != j) { 124 ans.push_back(pii(i, j)); 125 swap(linker[i], linker[j]); 126 } 127 } 128 printf("%d\n", ans.size()); 129 Rep(i, ans.size()) printf("C %d %d\n", ans[i].first, ans[i].second); 130 } 131 RT 0; 132 }
?
轉(zhuǎn)載于:https://www.cnblogs.com/kirai/p/5804288.html
總結(jié)
以上是生活随笔為你收集整理的[HDOJ2819]Swap(二分图最大匹配, 匈牙利算法)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
 
                            
                        - 上一篇: JZOJ #4722 跳楼机 (最短路模
- 下一篇: 《秋槿》第五句是什么
