(ищу алгоритм) Хопкрофта-Карпа

#include <algorithm>
#include <cstdio>
#include <vector>
using namespace std;
 
const int N = 15000;
int distx[N], disty[N], matex[N], matey[N], q[N];
vector<int> g[N];
 
bool BFS(int n)
{
    bool found = false;
    int fore = 0, rear = 0;
    for (int i = 0; i < n; ++i)
    {
        if (matex[i] == -1)
            q[rear++] = i;
        distx[i] = 0;
        disty[i] = 0;
    }
    for (; fore < rear; ++fore)
    {
        int x = q[fore];
        for (vector<int>::const_iterator it = g[x].begin(); it != g[x].end(); ++it)
            if (!disty[*it])
            {
                disty[*it] = distx[x]+1;
                if (matey[*it] == -1)
                    found = true;
                else
                    distx[matey[*it]] = disty[*it]+1, q[rear++] = matey
 
[*it];
            }
    }
    return found;
}
 
bool DFS(int x)
{
    for (vector<int>::const_iterator it = g[x].begin(); it != g[x].end(); ++it)
        if (disty[*it] == distx[x]+1)
        {
            disty[*it] = 0;
            if (matey[*it] == -1 || DFS(matey[*it]))
                return matex[x] = *it, matey[*it] = x, true;
        }
    return false;
}
 
int Hopcroft_Karp(int n)
{
    int res = 0;
    fill_n(matex, n, -1);
    fill_n(matey, n, -1);
    while (BFS())
        for (int i = 0; i < n; ++i)
            if (matex[i] == -1 && DFS(i))
                ++res;
    return res;
}
 
int main()
{
    int m;
    scanf("%d%d", &n, &m);
    for (; m; --m)
    {
        int u, v;
        scanf("%d%d",&u, &v);
        g[u].push_back(v);
    }
    printf("%d\n", Hopcroft_Karp(n));
}