二分图最大匹配算法-Hopcroft-Karp模板

时间复杂度：O((√V)*E)

#include<stdio.h>
#include<string.h>
const int N=,M=,INF=0x3f3f3f3f;
int dx[N],dy[M],sx[N],sy[M],p[N],q[N],a[N][M],l,r,n,m,d;
int bfs()
{
    l=r=;
    memset(dx,-,sizeof(dx));
    memset(dy,-,sizeof(dy));
    int i,u;d=INF;
    for(i=;i<=n;i++)
    {
        if(sx[i]==-)
        {
            q[++r]=i;
            dx[i]=;
        }
    }
    while(l<r)
    {
        u=q[++l];
        if(dx[u]>d)    break;
        for(i=;i<=m;i++)
        {
            if(a[u][i]&&dy[i]==-)
            {
                dy[i]=dx[u]+;
                if(sy[i]==-)    d=dy[i];
                else
                {
                    dx[sy[i]]=dy[i]+;
                    q[++r]=sy[i];
                }
            }
        }
    }
    return d!=INF;
}
int dfs(int u)
{
    for(int i=;i<=m;i++)
    {
        if(a[u][i]&&!p[i]&&dy[i]==dx[u]+)
        {
            p[i]=;
            if(sy[i]!=-&&dy[i]==d)    continue;
            if(sy[i]==-||dfs(sy[i]))
            {
                sy[i]=u,sx[u]=i;
                return ;
            }
        }
    }
    return ;
}
int HK_maxMatch()
{
    int ans=,i;
    memset(sx,-,sizeof(sx));
    memset(sy,-,sizeof(sy));
    while(bfs())
    {
        memset(p,,sizeof(p));
        for(i=;i<=n;i++)
        {
            if(sx[i]==-)
            {
                ans+=dfs(i);
            }
        }
    }
    return ans;
}

H-K