DFS版
#define INF 2147483644
const int N=200;
bool T[N],S[N];
int slack[N];
int w[N][N];
int Link[N];
int Lx[N],Ly[N];
int n;
bool match(int i)
{
S[i]=1;
for(int j=1; j<=n; j++)if(!T[j])
{
int t=Lx[i]+Ly[j]-w[i][j];
if(!t)
{
T[j]=true;
if(!Link[j]||match(Link[j]))
{
Link[j]=i;
return true;
}
}
else slack[j]=min(slack[j],t);
}
return false;
}
void update()
{
int a=INF;
for(int i=1; i<=n; i++)
if(!T[i]) a=min(a,slack[i]);
for(int i=0; i<=n; i++)
{
if(S[i]) Lx[i]-=a;
if(T[i]) Ly[i]+=a;
else slack[i]-=a;
}
}
void KM()
{
for(int i=1; i<=n; i++)
{
Link[i]=Lx[i]=Ly[i]=0;
slack[i]=INF;
for(int j=1; j<=n; j++)
{
Lx[i]=max(Lx[i],w[i][j]);
}
}
for(int i=1; i<=n; i++)
{
for(int j=1; j<=n; j++) slack[j]=INF;
while(true)
{
for(int j=1; j<=n; j++) S[j]=T[j]=0;
if(match(i)) break;
else update();
}
}
}
BFS版
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
#define INF 2147483644
const int N=20;
bool used[N];
bool T[N],S[N];
int fa[N];
int slack[N];
int w[N][N];
int Linky[N],Linkx[N];
int Lx[N],Ly[N];
int n;
void getmatch(int x)
{
while(fa[x])
{
int u=fa[x];
Linky[Linkx[x]]=u;
x=u;
}
}
bool match(int i)
{
queue<int> Q;
Q.push(i);
memset(fa,0,sizeof(fa));
while(!Q.empty())
{
int x=Q.front();
Q.pop();
S[x]=1;
for(int y=1; y<=n; y++)
{
if(T[y])continue;
int t=Lx[x]+Ly[y]-w[x][y];
if(t==0)
{
T[y]=true;
if(!Linky[y])
{
Linky[y]=x;
getmatch(x);
return true;
}
else
{
Q.push(Linky[y]);
fa[Linky[y]]=x;
}
}
else slack[y]=min(slack[y],t);
}
}
return false;
}
void update()
{
int a=INF;
for(int i=1; i<=n; i++)
if(!T[i]) a=min(a,slack[i]);
for(int i=0; i<=n; i++)
{
if(S[i]) Lx[i]-=a;
if(T[i]) Ly[i]+=a;
else slack[i]-=a;
}
}
void KM()
{
for(int i=1; i<=n; i++)
{
Linky[i]=Lx[i]=Ly[i]=0;
slack[i]=INF;
for(int j=1; j<=n; j++)
{
Lx[i]=max(Lx[i],w[i][j]);
}
}
for(int i=1; i<=n; i++)
{
for(int j=1; j<=n; j++) slack[j]=INF;
while(true)
{
for(int j=1; j<=n; j++) S[j]=T[j]=0;
memset(Linkx,0,sizeof(Linkx));
for(int j=1; j<=n; j++) Linkx[Linky[j]]=j;
if(match(i)) break;
else update();
}
}
}