匈牙利算法（二分图最大匹配）

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int i,j,k,m,n,e;
int hd[2001],temp;
int mat[2001],b[2001];
int ans;

struct data
{
int y,nxt;
}a[2000001];

int r()
{
int aans=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')
f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
aans*=10;
aans+=ch-'0';
ch=getchar();
}
return aans*f;
}

void add(int xx,int yy)
{
a[++temp].y=yy;
a[temp].nxt=hd[xx];
hd[xx]=temp;
}

bool dfs(int x)
{
for(int p=hd[x];p;p=a[p].nxt)
{
int son=a[p].y;
if(!b[son])
{
b[son]=1;//是否遍历过
if(!mat[son]||dfs(mat[son]))//如果目标搭档没有搭档或者目标搭档可以被抢走
{
mat[son]=x;//那么把让这两个成为一对
mat[x]=son;
return 1;//返回当前可行
}
}
}
return 0;//返回不可行
}

int main()
{
//  freopen("in.txt","r",stdin);
n=r(),m=r(),e=r();
int u,v;
for(i=1;i<=e;i++)
{
u=r(),v=r();
if(v<=m&&u<=n&&u>0&&v>0)//判断数据
{add(u,v+n);//建双向边
add(v+n,u);}
}

for(i=1;i<=n;i++)
{
if(!mat[i])//是否有搭档
{
memset(b,0,sizeof(b));
if(dfs(i))
ans++;
}
}

cout<<ans;

return 0;
}
/*
3 4 6
1 1
1 2
1 3
2 2
2 4
3 1
3 3
*/