匈牙利算法

理论依据：hall定理

#include<bits/stdc++.h>
#define maxn 2005
using namespace std;

int n,m,e;
int info[maxn],Prev[maxn*maxn],to[maxn*maxn],cnt_e;
inline void Node(int u,int v){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v; }
int cy[maxn];
bool vis[maxn];

bool dfs(int now)
{
    for(int i=info[now];i;i=Prev[i])
        if(!vis[to[i]])
        {
            vis[to[i]] = 1;
            if(!cy[to[i]] || dfs(cy[to[i]]))
            {
                cy[to[i]] = now;
                return 1;
            }
        }
    return 0;
}

int main()
{
    scanf("%d%d%d",&n,&m,&e);
    for(int i=1,u,v,w;i<=e;i++)
    {
        scanf("%d%d",&u,&v);
        if(v>m) continue;
        Node(u,v+n),Node(v+n,u);
    }
    int ans = 0;
    for(int i=1;i<=n;i++)
    {
        memset(vis,0,sizeof vis);
        ans += dfs(i);
    }
    printf("%d\n",ans);
}

HK算法，BFS优化，时间复杂度O(sqrt(|V|) * |E|)
可以证明在寻找增广路集的每一个阶段所寻找到的最短增广路都具有相等的长度，并且随着算法的进行最短增广路的长度是越来越长的，更进一步的分析可以证明最多只需要增广ceil(sqrt(n))次就可以得到最大匹配（证明在这里略去）。
实际上感觉对于随（水）机（题）数（尿）据（性）优化并不明显，
匈牙利可以n^3过1000

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define maxn 2005
using namespace std;

int n,m,e;
int d[maxn],c[maxn];
int info[maxn],Prev[maxn*maxn],to[maxn*maxn],cnt_e;
inline void Node(int u,int v){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v; }

int Q[maxn],L,R;
int BFS()
{
	int Min = 0x3f3f3f3f;
	memset(d,0,sizeof d);
	L=R=0;
	for(int i=1;i<=n;i++) if(!c[i]) Q[R++] = i,d[i] = 1;
	for(int now;L<R;)
	{
		if(d[now = Q[L++]] > Min) continue;
		for(int i=info[now];i;i=Prev[i])
			if(!d[to[i]]) 
			{
				d[to[i]] = d[now] + 1;
				if(!c[to[i]]) Min = d[to[i]];
				else d[c[to[i]]]=d[to[i]]+1,Q[R++] = c[to[i]];
			}
	}
	return Min;
}

bool dfs(int now)
{
	for(int i=info[now];i;i=Prev[i])
		if(d[to[i]] == d[now] + 1)
		{
			d[to[i]] = 0;
			if(!c[to[i]] || dfs(c[to[i]]))
			{
				c[c[to[i]] = now] = to[i];
				return 1;
			}
		}
	return 0;
}

int main()
{
	scanf("%d%d%d",&n,&m,&e);
	for(int u,v;e--;)
	{
		scanf("%d%d",&u,&v);
		if(u<=m && v<=m) Node(u,v+n),Node(v+n,u);
	}
	
	int ans = 0;
	while(BFS() < 0x3f3f3f3f)
	{
		for(int i=1;i<=n;i++)
			if(!c[i])ans += dfs(i);
	}
	
	printf("%d\n",ans);
}