[HDOJ2819]Swap（二分图最大匹配, 匈牙利算法）

首页 > [HDOJ2819]Swap（二分图最大匹配, 匈牙利算法）

[HDOJ2819]Swap（二分图最大匹配, 匈牙利算法）

题目链接：http://acm.split.hdu.edu.cn/showproblem.php?pid=2819

题意：给一张n*n的01矩阵，可以任意交换其中的行或者列，问是否可以交换出来一个对角线上都是1的矩阵。

按行列号建图，如果(i,j)为1的话，则i和j就有一条边。匹配出的结果可以认为如何交换使得行列相等，输出结果要注意不要忘记更新linker，以保证在正确的位置。

  1 /*
  2 ━━━━━┒ギリギリ♂ eye！
  3 ┓┏┓┏┓┃キリキリ♂ mind！
  4 ┛┗┛┗┛┃＼○／
  5 ┓┏┓┏┓┃ /
  6 ┛┗┛┗┛┃ノ)
  7 ┓┏┓┏┓┃
  8 ┛┗┛┗┛┃
  9 ┓┏┓┏┓┃
 10 ┛┗┛┗┛┃
 11 ┓┏┓┏┓┃
 12 ┛┗┛┗┛┃
 13 ┓┏┓┏┓┃
 14 ┃┃┃┃┃┃
 15 ┻┻┻┻┻┻
 16 */
 17 #include 
 18 #include 
 19 #include 
 20 #include 
 21 #include 
 22 #include 
 23 #include 
 24 #include 
 25 #include 
 26 #include 
 27 #include 
 28 #include 
 29 #include 
 30 #include 
 31 #include 
 32 #include <set>
 33 #include 
 34 #include 
 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 pll;
 65 typedef map<string, int> msi;
 66 typedef vector<int> vi;
 67 typedef vector vl;
 68 typedef vector 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 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
", ans.size());
129         Rep(i, ans.size()) printf("C %d %d
", ans[i].first, ans[i].second);
130     }
131     RT 0;
132 }

转载于:https://www.cnblogs.com/kirai/p/5804288.html

[HDOJ2819]Swap（二分图最大匹配, 匈牙利算法）

更多相关：

引用-ZIGBEE-ZSTACK网络配置相关问题

BootLoader引导程序制作及移植(一)

基于OHCI的USB主机 —— USB设备常量定义

PCL点云分割（2）

PCL点云配准（2）

PCL点云变换与移除NaN

PCL滤波介绍(3)

PCL滤波介绍（2）

java 读取excel_Java12POI操作Excel

C++ 生成随机数

基于C++的本征图像分解(Intrinsic Image Decomposition)

剑指offer：面试题39. 数组中出现次数超过一半的数字

剑指offer：面试题33. 二叉搜索树的后序遍历序列