[模板]二分图算法-Hopcroft-Carp

#include<bits/stdc++.h>
using namespace std;
const int maxn = 5e4+10;
const int INF = 0x3f3f3f3f;
vector<int> G[maxn];
/// mx  是 左边点集所匹配的右边的点 ，同理 my 是 右边点集所匹配的左边的点集，dx是源点到x点集的距离，dy是源点到y点集的距离。
int mx[maxn],my[maxn],dx[maxn],dy[maxn];
bool vis[maxn];
int cntx,cnty;
int dis;
bool SearchPath() {
    dis = INF;
    //bool flag = false;
    queue<int> qu;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=cntx;i++)if(mx[i] == -1) { /// 以左点集(未匹配的点)作为源点集，
        qu.push(i);dx[i] = 0;
    }
    while(!qu.empty())
    {
        int u = qu.front();qu.pop();
        if(dx[u] > dis) break;
        for(int i=0;i<G[u].size();i++) {
            int to = G[u][i];
            if(dy[to] == -1) { /// 如果to还没有访问到
                dy[to] = dx[u] + 1;
                if(my[to] != -1) {  ///如果to已经有了匹配的对象
                    dx[my[to]] = dy[to] + 1; /// 距离加1
                    qu.push(my[to]); /// 将to的对象加入队列.
                }
                else dis = dy[to]; ///to没有匹配的对象，则找到了增广路.
                //else flag = true;
            }
        }
    }
    return flag;
}
bool dfs(int u) { /// 类似与匈牙利的find函数
    for(int i=0;i<G[u].size();i++) {
        int to = G[u][i];
        if(!vis[to] && dy[to] == dx[u] + 1) ///这里不同，只需要距离为1的匹配。
        {
            vis[to] = true;
            if(my[to]!=-1 &&dy[to] == dis)continue;
            if(my[to]==-1||dfs(my[to])) { /// 如果to没有匹配，或者to的匹配节点还能找到另一个匹配.
                mx[u] = to;               /// 改变匹配关系。
                my[to] = u;
                return true;
            }
        }
    }
    return false;
}
int MaxMatch() {
    int ans = 0;
    memset(mx,-1,sizeof(mx));
    memset(my,-1,sizeof(my));
    while(SearchPath()) /// search the path
    {
        memset(vis,0,sizeof(vis));
        for(int i=1;i<=cntx;i++) {   /// dfs is match
            if(mx[i] == -1 && dfs(i)) ans++;
        }
    }
    return ans;
}