链式前向星版本:
int n,cnt,used[55],match[55];
int find(int x){
for(int i=head[x];i!=-1;i=edge[i].next){
int v=edge[i].to;
if(!used[v]){
used[v]=true;
if(!match[v]||find(match[v])){
match[v]=x;
return true;
}
}
}
return false;
}
洛谷P3386模板:
#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define per(i,n,a) for(int i=n;i>=a;i--)
#define endl '\n'
#define eps 0.000000001
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))
#define IO ios::sync_with_stdio(false);cin.tie(0);
using namespace std;
const int INF=0x3f3f3f3f;
const ll inf=0x3f3f3f3f3f3f3f3f;
const int mod=1e9+7;
const int maxn=1e5+5;
int n,m,e,line[505][505],used[505],match[505];
int find(int x){
rep(i,1,m){
if(!used[i]&&line[x][i]){
used[i]=true;
if(match[i]==0||find(match[i])){
match[i]=x;
return true;
}
}
}
return false;
}
int main(){
cin>>n>>m>>e;
rep(i,1,e){
int u,v;cin>>u>>v;
line[u][v]=1;
}
int ans=0;
rep(i,1,n){
mem(used,0);
if(find(i)) ++ans;
}
cout<<ans<<endl;
}
以hdoj2063为例封装模板:
#include<bits/stdc++.h>
#define ll long long
#define endl '\n'
#define mem(a) memset(a,0,sizeof(a))
#define IO ios::sync_with_stdio(false);cin.tie(0);
using namespace std;
const int maxn=505;
ll line[maxn][maxn],boy[maxn],used[maxn];
ll k,n,m;
ll Find(ll x){
for(ll i=1;i<=m;i++){
if(line[x][i]==1&&used[i]==0){
used[i]=1;
if(boy[i]==0||Find(boy[i])){
boy[i]=x;
return true;
}
}
}
return false;
}
int main(){
while(cin>>k>>n>>m&&k){
memset(line,0,sizeof(line));
memset(boy,0,sizeof(boy));
memset(used,0,sizeof(used));
for(ll i=1;i<=k;i++){
ll u,v;cin>>u>>v;
line[u][v]=1;
}
ll sum=0;
for(ll i=1;i<=n;i++){
memset(used,0,sizeof(used));
if(Find(i)) sum++;
}
cout<<sum<<endl;
}
}
以上就是代码的注释(不太喜欢封装完后的代码还有注释乱七八糟的就自己另写了一份注释)。
匈牙利算法是二分图匹配的一种,是基于网络流的一个算法。主要是递归和回溯的过程比较难以理解需要加以时间去理解。