简述匈牙利算法

二分图的最大匹配——匈牙利算法

【基本概念】：

二分图：

二分图又称作二部图，是图论中的一种特殊模型。 设G=(V,E)是一个无向图，如果顶点V可分割为两个互不相交的子集(A,B)，并且图中的每条边（i，j）所关联的两个顶点i和j分别属于这两个不同的顶点集(i in A,j in B)，则称图G为一个二分图。 

无向图G为二分图的充分必要条件是，G至少有两个顶点，且其所有回路的长度均为偶数。 

最大匹配：

给定一个二分图G，在G的一个子图M中，M的边集中的任意两条边都不依附于同一个顶点，则称M是一个匹配.　　选择这样的边数最大的子集称为图的最大匹配问题，如果一个匹配中，图中的每个顶点都和图中某条边相关联，则称此匹配为完全匹配，也称作完备匹配.  

借某位大牛的话，匈牙利算法就一句话：有机会上，没机会创造机会也要上。

//伪代码
bool 寻找从k出发的增广路
{
    for (k的所有邻接点j)
    {
        if (j不在增广路上)
        {
            把j加入增广路;
            if (j是非-M饱和点 或者 从j的对应项出发有可增广路)
            {
                修改j的对应项为k;
                返回true;
            }
        }
    }
    返回false;
}

void 匈牙利hungary()
{
    for i->1 to n
    {
        if (则从i出发有增广路)
            匹配数++;
    }
    输出 匹配数;
}

#include<cstdio>
#include<cstring>
using namespace std;
const int N = 105;
int matrix[N][N],link[N],vis[N];
int X,Y,m;
int ans;
int dfs(int x){
    for(int i=1;i<=Y;i++){
        if(matrix[x][i]!=0&&!vis[i]){
            vis[i]=1;
            if(!link[i]||dfs(link[i])){
                link[i]=x;
                return 1;
            }
        }
    }
    return 0;
}
void Edmonds(){
    memset(link,0,sizeof(link));
    int i,j;
    for(i=1;i<=X;i++){
        memset(vis,0,sizeof(vis));
        if(dfs(i)) ans++;
    }
    printf("the number of maximal matching is: %d\n",ans);
    printf("this is one of them:\n");
    for(i=1;i<=X;i++){
        if(link[i]!=0){
            printf("x%d->y%d\n",link[i],i);
        }
    }
}
int main(){
    int i,j;
    while(scanf("%d%d%d",&X,&Y,&m)!=EOF){
        int x,y;
        ans=0;
        for(i=1;i<=m;i++){
            scanf("%d%d",&x,&y);
            matrix[x][y]=1;
        }
        //memset(link,0,sizeof(link));
        Edmonds();
    }
    return 0;
}