HDU 2458 Kindergarten ：二分图用匈牙利求最大独立集

前言

（悄咪咪说一句 好像不是dp， 但是是目前收获挺多的一道题了，一开始hhhh一直列不出来dp[ G ] [ B ] 怎么转移，后面犟太久了顶不住看了题解，发现需要学二分图、最大匹配、匈牙利算法，听都没听过……我两行清泪流下来，为什么当初要犟那么久才看题解）

知识点

这道题我主要是学了下面这篇博客，真的写的超级好，每一步都很细致，通俗易懂又深入浅出 nice ~
夜深人静写算法（八）- 二分图最大匹配：https://blog.csdn.net/WhereIsHeroFrom/article/details/112479644

下面贴一些和这道题有关的概念和结论：

1、二分图

二分图简单来说就是一个图的所有点可以分成两个独立的集合U、V，这里所谓的独立是指集合 U 所有点之间都没有边相连，V 也同理。但是 U 和 V 之间可能会有边相连。

比如下图：

2、匹配

定义：给定一个二分图 G，在 G 的一个子图 M 中，M 的边集 E 中的任意两条边都不依附于同一个顶点，则称 M 是一个匹配。

简单来说就是假如 M 中有 a 到 b 的关系，那么 M 中就不能有任何其它以 a 或者 b 为端点的边。
从匹配的定义可以看出，匹配也是一个图，并且图中的 点数 是 边数 的 2 倍；

3、最大匹配

定义：最大匹配就是找到一个子图，满足是匹配，并且边数（点数）最多。
（一般用匈牙利算法求出）

4、匈牙利算法

（后面的例题是该算法的一个模板）

简单来说，这就是一个讲究时间顺序、对渣男横刀夺爱、对好男人忍痛割爱的爱情故事。

在一个美丽的小镇上，有一群单身男女，他们碍于当时落后的思想，不敢积极主动追求爱情，但是随着新思想不断涌进，他们也日渐开放起来（这里默认是只有男生追女孩子），敢于对平时有好感的异性展开追求（都说了是美好的小镇了，所以我们默认追了就能追到手）。

但我们也知道，爱情是分先来后到的，总会有人比你先勇敢一步，也就是最前面所说的 “ 讲究时间顺序 ”，那么我们假设这个小镇有 4 个男生 5 个女生，对他们分别编号b1、b2、……，g1、g2……（boy 和 girl 的意思），其中 b1、b2……是按照谁先勇敢谁就排在前面的规则排序。

虽然小镇美好，但是人不一定都是好人，一个男生他可能同时对多个女生有好感，现在我们将男女的好感关系列出来（有连线即彼此互相有感觉）：

现在有一个 “ 好朋友 ” 聚会，顾名思义，参加聚会的人不能有超过朋友以上的情感，换句话说，参加聚会的人在上面的关系图中不能存在边相连。现在想问，这个小镇最多有多少人能参加？

问题的答案就是男生数量 + 女生数量 - 该二分图的最大匹配数。因为对二分图去掉最大匹配中每一个匹配的一个端点，最终得到的二分图就是最多顶点且没有任何边相连的子图，而这个子图我们称之为最大独立集。

为什么是“ 男生数量 + 女生数量 - 该二分图的最大匹配数 ” 呢？我们可以反证：假如我们现在已经得到了最大匹配数，根据最大匹配数的定义，此时但凡再多加一条边，都会使得最大匹配不成立，而且增加的边必定有一个端点包含于原先的最大匹配中，那么这条新加的边完全可以靠去除那个相同的端点而一同去除，没有必要再多加一条边进来。

好，那么现在我们再解决一个为什么是只减去该二分图的最大匹配数，而不用 * 2 呢？一条边的两个端点一起去除不是更保证没关系吗？ 这是因为要想消除最大匹配里面的所有边，只需要把一个端点去掉，那么两个端点之间的边也会跟着消除，不必要两个端点同时消失，因为我们要求的是最大独立集，当然是想要满足要求的点最多。

用上面那个故事来说就是，我们把最大数量能参加好朋友聚会的男男女女当成那个最大独立集，因为该集合下，没有任何一个人对另一个有好感。而上面所谓的最大匹配就对应于这个小镇最多能有几对终成眷属，对应那个二分图，就是 b1 - g3、b2 - g5、b3 - g1，此时最划算，因为能凑成最多对情侣，那么我们也由此可以发现，只要去除了 b1、b3、g5 这 3 个人，那么剩下的就是最多人参加聚会的人。（不一定是全去除左端点或者全去除右端点，就像这个例子一样。不过这不影响，因为每条匹配只需要去掉一个端点就行了，至于去掉的是哪个端点，并不影响结论）。

好，罗里吧嗦讲完了，那么在代码上这个集合怎么求呢？
还记得我们一开始说这是一个讲究时间顺序、对渣男横刀夺爱、对好男人忍痛割爱的爱情故事吗？横刀夺爱和忍痛割爱我个人感觉hhhhhh就是匈牙利算法的精髓hhhhhhh。（把上面的图再放一遍）

是这样子的，男生 b1，是男生中最先主动的人，那么现在他追到 g1 了，两个人谈起了甜甜的恋爱；
第二个主动的是 b2 ，刚刚好，他主动的也还算早，g5 还单身，所以他也追到了，b2 也和 g5 谈起了恋爱；
但是 b3 就太闹心了，喜欢也不主动点，他发现 g1 已经和 b1谈起了恋爱，但是他又只喜欢 g1， 所以没有办法劝自己下一个更乖悄然离去，而且很重要的一点是，他发现 b1 很渣，虽然和 g1 谈恋爱，但是同时和 g3、g4 搞暧昧，g1 已经因为这些事和他吵了很多次了，于是 b3 就觉得，哎呀机会来了！趁着两个人冷战 b1 去找 g3 诉苦，b3 也去听 g1 诉苦，一来二往，b3 就横刀夺爱和 g1 在一起了，而渣男 b1 也和 g3 在一起了，至此男生 1 2 3号都有了对象啦；
接下来就到最后一颗独苗了， 可是 b4 刚想迈出那一步，结果发现已经被 b2 捷足先登了！而且 b2 是个专一的人，他并没有机会下手，于是就忍痛割爱，为自己的慢一步买单去了。

到这里我们可以发现，为什么 b3 能得手，b4 不能？因为假如 b3 和 g1 在一起了，b1 可以找别人，换句话说，这样做能增加情侣的对数。而就算 b4 最终从 b2 手中抢下 g5，也没有其他女生能和 b2 匹配，总的情侣对数并不会增加，倒不如退一步成人之美是吧。

而这里得到的最多情侣对数就对应于我们想求的最大匹配数，有了最大匹配数，最小顶点覆盖、最大独立集就套公式可以求出。

5、最小顶点覆盖

定义：选了一个点就相当于覆盖了以它为端点的所有边。最小顶点覆盖就是选择最少的点来覆盖整张图所有的边。

结论：二分图最小顶点覆盖 = 二分图最大匹配

6、最大独立集

定义：选取最多的点，使得选中的点在图中任意两点都没有关系。

结论：二分图的最大独立集 = 顶点总数 - 最小顶点覆盖 = 顶点总数 - 二分图最大匹配。

例题 - HDU 2458 Kindergarten

原题链接：http://acm.hdu.edu.cn/showproblem.php?pid=2458

题目大意

幼儿园有一个班，班上有 G 个女生、B 个男生，其中女生两两之间互相认识，男生也是两两认识。男生编号 1 、2、……，女生同理。现在存在 M 对关系，关系 1 - 2 指女生1 和男生2认识。问该班上最多有几个人两两之间互为认识。
G、B、M 皆不超过 200。

输入格式
G B M
M行……
0 0 0 输入结束

输出格式
Case 1: 答案
……

输入
2 3 3
1 1
1 2
2 3
2 3 5
1 1
1 2
2 1
2 2
2 3
0 0 0

输出
Case 1: 3
Case 2: 4

思路

显然，a 和 b 认识，b 和 c 认识并不能推出 a 和 c 认识，所以这道题不能用并查集。但是 女女、男男之间互相认识，如果把女生男生各自当成一个点集，显然不能构成二分图。但是因为我们是要求两两认识的最大数，那么我们就是要把不认识的最大匹配数求出来，那么 男 + 女 - 不认识的最大匹配数 = 认识的最大独立集。

所以基于上面对最大独立集求法的分析，我们按照建图 - 匈牙利算法求最大匹配 - 输出顺序完成。

代码

#include <iostream>
#include <algorithm>
#include <iomanip>
#include <sstream>
#include <string>
#include <stack>
#include <queue>
#include <deque>
#include <vector>
#include <map>
#include <set>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <climits>

using namespace std;

#define getlen(array) {return (sizeof(array) / sizeof(array[0]));}
#define ll long long 
#define MEM(x, y) memset(x, y, sizeof x)
#define rin scanf("%d", &n)
#define rln scanf("%lld", &n)
#define rit scanf("%d", &t)
#define sc scanf
#define pr printf

const int INF = 0x3f3f3f3f;
const int N = 10000; 

//(val & 1) == 0偶， == 1奇。

bool graph[210][210];  //图
int pre[210];    //匹配
bool book[210];  //判断当前的男生是否遍历过
int g, b, m;

//找看看能不能增加匹配数，能 - true，不能 - false
bool findMatch(int k) {
	for (int i = 1; i <= b; ++ i) {
		// if(两个人认识 && 当次还没遍历过的男生)
		if (graph[k][i] == false && book[i] == false) {
			book[i] = true;    //true表示当次不必再在i号男生查询
			//if(直接匹配 || 让原先和i号男生匹配的女生找别的男生匹配)
			if (pre[i] == 0 || findMatch(pre[i])) {
				pre[i] = k;  //更新和i男生匹配的女生
				return true;  //能增加匹配数
			}
		}
	}
	return false;  //不能增加匹配数
}

int main() {
	//freopen("D:\\in.txt", "r", stdin); 
	//freopen("D:\\out.txt", "w", stdout);
	int t = 1;
	while (sc("%d %d %d", &g, &b, &m) != -1 && g + b + m > 0) {
		//置零
		MEM(graph, false);
		MEM(pre, 0);
		//建图
		int x, y;
		for (int i = 0; i < m; ++ i) {
			sc("%d %d", &x, &y);
			graph[x][y] = true;
		}
		///计算最大匹配
		int cnt = 0;
		for (int i = 1; i <= g; ++ i) {
			//每一次匹配，都应当对男生重置，因为从头到尾把男生走一遍才能确保不会漏
			MEM(book, false);  
			if (findMatch(i)) 
				++ cnt;
		}
		//输出
		pr("Case %d: %d\n", t ++, g + b - cnt);
	}
	
	return 0;
}