Hopcroft-Carp 算法模板 自用

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int maxn = , INF = 0xfffffff;
int cx[maxn], cy[maxn], dx[maxn], dy[maxn], used[maxn], g[maxn][maxn];
int nx, ny, dis;
int bfs()
{
    queue<int> Q;
    dis = INF;
    mem(dx,-);
    mem(dy,-);
    for(int i=; i<=nx; i++)
    {
        if(cx[i] == -)
        {
            Q.push(i);
            dx[i] = ;
        }
    }

    while(!Q.empty())
    {
        int u = Q.front(); Q.pop();
        if(dx[u] > dis) break;
        for(int i=; i<=ny; i++)
        {
            if(g[u][i] && dy[i] == -)
            {
                dy[i] = dx[u] + ;
                if(cy[i] == -) dis = dy[i];
                else
                {
                    dx[cy[i]] = dy[i] + ;
                    Q.push(cy[i]);

                }
            }
        }
    }
    return dis != INF;
}

int dfs(int u)
{
    for(int i=; i<=ny; i++)
    {
        if(g[u][i] && !used[i] && dy[i] == dx[u] + )
        {
            used[i] = ;
            if(cy[i] != - && dis == dy[i]) continue;
            if(cy[i] == - || dfs(cy[i]))
            {
                cy[i]= u;
                cx[u] = i;
                return ;
            }
        }
    }
    return ;
}
int hc()
{
    int res = ;
    mem(cx,-);
    mem(cy,-);
    while(bfs())
    {
        mem(used,);
        for(int i=; i<=nx; i++)
        {
            if(cx[i] == - && dfs(i))
                res++;
        }
    }
    return res;
}
int main()
{

    return ;
}