匈牙利算法

匈牙利算法

用于二分图匹配,不可用于一般图,一般图要用带花树.来做.
一个一个点来操作,为它先找一个点(任意可以匹配的一点),若其他点已经使用过了,则递归判断是否可以让使用该点的点再找一个点,并用同样的方法判断能否找到
举例说明:A可以匹配a点,但此时B已经与a匹配,则判断B是否可以换一个点.
若B可以匹配b点,而C已经与b匹配,则判断C是否可以换一个点.
如此递归下去,记得a,b再递归过程中不能再用了.
用邻接矩阵的话时间复杂度最劣为O(n^3),用邻接表的时间复杂度为O(nm)

代码(以洛谷 P2756 飞行员配对方案问题为例)

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int n,m,mm[110][110],pp[110],ans,an[110];
bool use[110];

bool find(int u)
{
    int i,j;
    for(i=1;i<=n;i++)
    {
        if(!mm[u][i]||use[i]) continue;
        use[i]=1;
        if(!pp[i]||find(pp[i]))
        {
            pp[i]=u;
            an[u]=i;
            return 1;
        }
    }
    return 0;
}

int main()
{
    int i,j,p,q;
    cin>>m>>n;
    for(;;)
    {
        scanf("%d%d",&p,&q);
        if(p==-1&&q==-1) break;
        q-=m;
        mm[p][q]=1;
    }
    for(i=1;i<=m;i++)
    {
        memset(use,0,sizeof(use));
        if(find(i)) ans++;
    }
    cout<<ans<<endl;
    for(i=1;i<=m;i++)
    {
        if(an[i])
        {
            printf("%d %d\n",i,an[i]+m);
        }
    }
}

此题也可以用最大流来解.