KM(Kuhn-Munkres)算法求带权二分图的最佳匹配
相关概念
这个算法个人觉得一开始时有点难以理解它的一些概念，特别是新定义出来的，因为不知道是干嘛用的。但是，在了解了算法的执行过程和原理后，这些概念的意义和背后的作用就渐渐的显示出来了。因此，先暂时把相关概念列出来，看看，有个大概印象就好，等到了解了算法的流程后，在看原理中会有这些概念，那个时候回来细看就好了。
完备匹配：定义 设G＝
二分图最优匹配：对于二分图的每条边都有一个权（非负），要求一种完备匹配方案，使得所有匹配边的权和最大，记做最优完备匹配。（特殊的，当所有边的权为1时，就是最大完备匹配问题）
二分图带权匹配与最优匹配：什么是二分图的带权匹配？二分图的带权匹配就是求出一个匹配集合，使得集合中边的权值之和最大或最小。而二分图的最优匹配则一定为完备匹配，在此基础上，才要求匹配的边权值之和最大或最小。二分图的带权匹配与最优匹配不等价，也不互相包含。
以下是一些转换的思路：
KM算法是求最大权完备匹配，如果要求最小权完备匹配怎么办？方法很简单，只需将所有的边权值取其相反数，求最大权完备匹配，匹配的值再取相反数即可。
KM算法的运行要求是必须存在一个完备匹配，如果求一个最大权匹配(不一定完备)该如何办？依然很简单，把不存在的边权值赋为0。
KM算法求得的最大权匹配是边权值和最大，如果我想要边权之积最大，又怎样转化？还是不难办到，每条边权取自然对数，然后求最大和权匹配，求得的结果a再算出e^a就是最大积匹配。至于精度问题则没有更好的办法了。
下面是算法需要而建立的概念：
顶标：每个节点与另一个集合中节点之间的最大权值
可行顶标：对于原图中的任意一个结点，给定一个函数L(node)求出结点的顶标值。我们用数组lx(x)记录集合X中的结点顶标值，用数组ly(y)记录集合Y中的结点顶标值。 并且，对于原图中任意一条边\(edge(x,y)\)，都满足\(lx(x)+ly(x)>=weight(x,y)\)。
相等子图：设 G(V,E) 为二部图， G'(V,E') 为二部图的子图。如果对于 G' 中的任何边
算法流程
这里有一篇比较好的男女找对象指南博客讲解了KM算法的执行过程，生动而又形象。主要理解过程中是怎么满足配对条件的，不满足后怎么办，过程中降低期望值d是怎么求得的，在配对失败期望值进行更改以后发生了那些变化，这些变化为新的一轮的匹配带来了什么？
此处有链接
算法原理：
搞清楚算法流程后对算法的原理就会很好奇（可能吧~），为什么要这样做，这样做后为什么会产生这些改变，为什么利用这些改变就可以得出最终结果？
在详细介绍原理前，可以先看一下这篇博客（有图），根据算法流程介绍原理，可以直观感受一下，写的还是很清晰的，就是字有点小
此处有链接
算法的正确性基于以下定理：(相关概念就可以倒回去看了)
若由二分图中所有满足\(lx(x)+ly(y)==weight(x,y)\)的边(x,y)构成的子图（称做相等子图）有完备匹配，那么这个完备匹配就是二分图的最大权匹配。
因为对于二分图的任意一个匹配，如果它包含于相等子图，那么它的边权和等于所有顶点的顶标和；如果它有的边不包含于相等子图，那么它的边权和小于所有顶点的顶标和（即不是最优匹配）。所以相等子图的完备匹配一定是二分图的最大权匹配。
初始时为了使\(lx(x)+ly(y)>=weight(x,y)\)恒成立，令lx(x)为所有与顶点Xi关联的边的最大权，ly(y)=0。如果当前的相等子图没有完备匹配，就按下面的方法修改顶标以使扩大相等子图，直到相等子图具有完备匹配为止。
我们求当前相等子图的完备匹配失败了，是因为对于某个X顶点，我们找不到一条从它出发的交错路。这时我们获得了一棵交错树，它的叶子结点全部是X顶点（因为没法跟Y匹配下去了，终止在X）。现在我们把交错树中X顶点的顶标全都减小某个值d，Y顶点的顶标全都增加同一个值d（注意是交错树中的顶点，也就是匹配的女生和被访问的男生），那么我们会发现：（这个过程就对应于男女生配对中女生配对失败后的情况）
1）两端都在交错树中的边(x,y)，lx(x)+ly(y)的值没有变化。也就是说，它原来属于相等子图，现在仍属于相等子图。
2）两端都不在交错树中的边(x,y)，lx(x),ly(y)都没有变化。也就是说，它原来属于（或不属于）相等子图，现在仍属于（或不属于）相等子图。
3）X端不在交错树中，Y端在交错树中的边(x,y)，它的lx(x)+ly(y)的值有所增大。它原来不属于相等子图，现在仍不属于相等子图。
4）X端在交错树中，Y端不在交错树中的边(x,y)，它的lx(x)+ly(y)的值有所减小。也就说，它原来不属于相等子图，现在可能进入了相等子图，因而使相等子图得到了扩大。
现在的问题就是求d值了。为了使\(lx(x)+ly(y)>=weight(x,y)\) 始终成立，且至少有一条边进入相等子图，d应该等于：{ \(Minlx(x)+ly(y)-weight(x,y)|x在交错树中，y不在交错树中\)}。
以上是KM算法的基本思想。但是朴素的实现方法，时间复杂度为O(n4)——需要找\(O(n)\)次增广路，每次增广最多需要修改\(O(n)\)次顶标，每次修改顶标时由于要枚举边来求d值，复杂度为\(O(n^2)\)。实际上KM算法的复杂度是可以做到\(O(n^3)\)的。我们给每个Y顶点一个“松弛量”函数slack，每次开始找增广路时初始化为无穷大。在寻找增广路的过程中，检查边(x,y)时，如果它不在相等子图中，则让slack[y]变成原值与weight(x,y)-lx(x)-ly(y)的较小值（就是男生女生配对的求d的过程）。这样，在修改顶标时，取所有不在交错树中的Y顶点的slack值中的最小值作为d值即可。但还要注意一点：修改顶标后，要把所有的不在交错树中的Y顶点的slack值都减去d（对应就是未被访问的男生离女生更近了一步）。
代码
以下代码基于HDU2255奔小康赚大钱
#include
#include
#include
#define max_n 305
#define INF 0x3f3f3f3f
using namespace std;
int n;//一个集合中的顶点数
int love[max_n][max_n];//两个不同集合点之间的权值，即男女好感度
int ex_girl[max_n];//女生顶标
int ex_boy[max_n];//男生顶标
int vis_girl[max_n];//记录女生是否参与匹配
int vis_boy[max_n];//记录男生是否被尝试匹配
int slack[max_n];//松弛函数，以求出顶标的最小改变量d
int match[max_n];//记录匹配关系
int dfs(int girl)
{
vis_girl[girl] = 1;
for(int boy = 0;boy-1)
{
ans += love[match[i]][i];
}
}
return ans;
}
int main()
{
while(~scanf("%d",&n))
{
for(int i = 0; i 参考文章
x_y_q_ ，带权二分图的最佳匹配（KM算法）， https://blog.csdn.net/x_y_q_/article/details/51927054 （偏向原理）
SixDayCoder，二分图的最佳完美匹配——KM算法，https://blog.csdn.net/sixdaycoder/article/details/47720471 （概念向）
wenr，KM算法详解+模板，https://www.cnblogs.com/wenruo/p/5264235.html （算法流程介绍）
kuangbin， Kuhn-Munkres算法（二分图最大权匹配），https://www.cnblogs.com/kuangbin/archive/2012/08/19/2646535.html （总结向）