二分匹配Hopcroft-Karp算法

关注

原创最新推荐文章于 2025-12-28 09:24:01 发布 734 阅读 AI 写同款· GEO 优化 ›

深入探讨二分匹配算法及其应用实例在图论中，二分图（Bipartite Graph）是一种特殊类型的图，在该图中顶点可以被划分为两个互不相交的集合，使得图中的每条边连接的两个顶点分别属于这两个集合。这种结构常见于许多实际应用中，例如社交网络中的好友关系或企业与员工之间的雇佣关系。首先，我们需要将问题转化为二分图。在图中，每个用户表示为一个顶点，并且用户之间的好友关系表示为连接两个顶点的边。因此，如果两个用户是好友，那么在图中就会有一条边连接他们。阅读详情

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int N = 110;
const int INF=0x7fffffff;
int bmap[N][N];//每次记得清空
//对于两边各50000个点，200000条边的二分图最大匹配可以在1s内出解，效果很好：）
struct HK{//用HK算法求解二分最大匹配问题(时间复杂度为sqrt(v)*e)
    int nx,ny,dis;
    int cx[N];//cx[i]表示左集合i顶点所匹配的右集合的顶点序号
    int cy[N]; //cy[i]表示右集合i顶点所匹配的左集合的顶点序号
    int dx[N];
    int dy[N];
    int Q[N];
    bool bmask[N];
    bool searchpath(){//极大最短增广路集
        int front,rear;
        front=rear=0;//初始化队列
        dis=INF;
        memset(dx, -1, sizeof(dx));
        memset(dy, -1, sizeof(dy));
        for(int i=1;i<=nx;i++){
            if(cx[i]==-1)//左边未匹配
            {
                Q[rear++]=i;//将左边节点放入队列
                dx[i]=0;//分层 所有潜在起点构成0层
            }
        }
        
        while(front<rear){
            int u=Q[front++];//u都是未匹配的
            if(dx[u]>dis)break;
            for(int v=1;v<=ny;v++){
                if(bmap[u][v]&&dy[v]==-1){//相邻未分层构成i+1层。
                    dy[v]=dx[u]+1;//v对应的距离 为u对应距离加1
                    if(cy[v]==-1)dis=dy[v];//右侧未匹配，更新最短距离
                    else{
                        dx[cy[v]]=dy[v]+1;//右侧已经匹配，根据已匹配边回溯。分层
                        Q[rear++]=cy[v];//左侧已匹配点入队。
                    }
                }
            }
        }
        return dis!=INF;
    }
    int findpath(int u){//只搜索极大最短增广路集中的边
        for(int v=1;v<=ny;v++){
            if(!bmask[v]&&bmap[u][v]&&dy[v]==dx[u]+1){//反向搜索的过程
                bmask[v]=1;
                if(cy[v]!=-1&&dy[v]==dis)continue;
                if(cy[v]==-1||findpath(cy[v])){
                    cy[v]=u;
                    cx[u]=v;
                    return 1;
                }
            }
        }
        return 0;
    }
    int Max_match(){
        int res=0;
        memset(cx,-1,sizeof(cx));
        memset(cy,-1,sizeof(cy));
        while(searchpath()){
            memset(bmask, 0, sizeof(bmask));
            for(int i=1;i<=nx;i++){
                if(cx[i]==-1){
                    res+=findpath(i);
                }
            }
        }
        return res;
    }
};
HK hk;
//调用的时候确定hk.nx,hk.ny,然后直接调用Max_match()就行

收录于