匈牙利算法模板（无算法讲解）

二分图最大匹配和匈牙利算法想学习知识点的误入

二分图及匹配的定义不介绍（还是那句话，网上有很多相关资料，说白了就是我懒）

匈牙利算法，总的来说就是对于每个点出发找增广路，找到一个非匹配点就连，如果找不到就去在增广路上抢别的点的匹配点，然后让别的点换地方。（网上还是有很多资料）

其实这篇文章我就是想来贴两道陈芝麻烂谷子的模板题。。。

模板题1

https://www.luogu.org/problem/show?pid=3386 【模板】二分图匹配

这道题目不会的，应该去看所谓的“很多资料”。

题解：匈牙利算法（您略长、复杂度玄学的Dinic算法等也可）

代码：

#include <bits/stdc++.h>
using namespace std;
const int N = 1005;
int cnt, head[N<<1], mth
; bool vis
;
struct edge {int to, next;} e[N*N];
void addedge (int u, int v) {
e[++cnt] = (edge){v, head[u]}; head[u] = cnt;
}
bool hgy (int x) {
for (int i = head[x]; i; i = e[i].next)
if (!vis[e[i].to]) {
vis[e[i].to] = 1;
if (!mth[e[i].to] || hgy (mth[e[i].to])) {
mth[e[i].to] = x; return 1;
}
}
return 0;
}
int main () {
int n, m, ec, u, v, ans = 0;
scanf ("%d%d%d", &n, &m, &ec);
for (int i = 1; i <= ec; ++ i) {
scanf ("%d%d", &u, &v);
if (u<=n && v<=m) addedge (u, v);
}
for (int i = 1; i <= n; ++ i) {
memset (vis, 0, sizeof (vis));
if (hgy (i)) ++ ans;
} printf ("%d\n", ans);
return 0;
}

来一道升级一点的模板

模板题2

http://www.lydsy.com/JudgeOnline/problem.php?id=2663 [Beijing wc2012]灵魂宝石

R对K有单调性，所以可以二分。rmin直接跑最大匹配就行；rmax可以理解成求一个最小的极大匹配（不知道对不对），大概就是把分属于两个集合在原图中没有边的点在新图连边，有边的新图中不连边，然后用(n-最大匹配数)与k的大小关系作为判断的条件。大概是因为如果这样的匹配出来是k，说明对于当前的R，不可能有小于k的匹配满足条件，这样就能算出r的最大值。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef double DB;
const int N = 105;
const DB eps = 1e-8;
int head
, cnt, n, k, mth
;
bool vis
;
struct edge {int to, next;} e[N*N];
struct point {int x, y;} a
, b
;
int cmp (DB a, DB b) {
if (a-b>=0 && a-b<=eps) return 0;
if (b-a>=0 && b-a<=eps) return 0;
if (a-b>0) return 1; return -1;
}
int judis (point A, point B, DB d) {
int dx = A.x-B.x, dy = A.y-B.y;
DB dis = 1.0*dx*dx+dy*dy;
return cmp (dis, d*d);
}
void addedge (int u, int v) {
e[++cnt] = (edge){v, head[u]}; head[u] = cnt;
}
void init () {
cnt = 0;
memset (head, 0, sizeof (head));
memset (mth, 0, sizeof (mth));
for (int i = 1; i <= n*n; ++ i) {
e[i].to = e[i].next = 0;
}
}
bool hgy (int x) {
for (int i = head[x]; i; i = e[i].next)
if (!vis[e[i].to]) {
vis[e[i].to] = 1;
if (!mth[e[i].to] || hgy (mth[e[i].to])) {
mth[e[i].to] = x; return 1;
}
}
return 0;
}
bool Match (bool fl, DB r) {
int ans = 0; init ();
for (int i = 1; i <= n; ++ i)
for (int j = 1; j <= n; ++ j) {
if (fl && judis (a[i], b[j], r)<=0) addedge (i, j);
if (!fl && judis (a[i], b[j], r)>=0) addedge (i, j);
}
for (int i = 1; i <= n; ++ i) {
memset (vis, 0, sizeof (vis));
if (hgy (i)) ++ ans;
}
if (fl && ans>=k) return 1;
if (!fl && ans>=n-k) return 1;
return 0;
}
int main () {
DB l, r, mid; scanf ("%d%d", &n, &k);
for (int i = 1; i <= n; ++ i)
scanf ("%d%d", &a[i].x, &a[i].y);
for (int i = 1; i <= n; ++ i)
scanf ("%d%d", &b[i].x, &b[i].y);
l = 0, r = 1e6;
for (int i = 1; i <= 45; ++ i) {
mid = (l+r)/2;
if (Match (1, mid)) r = mid;
else l = mid;
}
printf ("%.2lf ", l);
if (n==k) {puts ("+INF"); return 0;}
l = 0, r = 1e6;
for (int i = 1; i <= 45; ++ i) {
mid = (l+r)/2;
if (Match (0, mid)) l = mid;
else r = mid;
}
printf ("%.2lf\n", l);
return 0;
}