听说这个算法在NOIP也比较重要?那就学习一下吧,不过由于时间有限,就直接转发别人的吧。以后有时间了(如果还没退役的话),我再重新写一份吧。
二分图的最大匹配
二分图最大匹配总结
下面是用匈牙利算法做的二分图的最大匹配(省略掉了一个数组)
#include <bits/stdc++.h>
using namespace std;
const int maxn=12306;
int lnk[maxn],used[maxn];
vector<int> divide[maxn];//领接矩阵好写啊
int n,m,e,ans;
inline int read() {//忽略掉这个快读
int x=0; char ch=getchar();
while(ch<'0'||ch>'9')ch=getchar();
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x;
}
bool Find(int x,int f){
for(int i=0;i<divide[x].size();i++){
int v=divide[x][i];
if(used[v]!=f){//当前匹配的不是f,本轮还没用过
used[v]=f;
if(lnk[v]==-1||Find(lnk[v],f)){//可以让步或名花无主
lnk[v]=x;//匹配
return 1;
}
}
}
return 0;
}
int main(){
n=read();m=read();e=read();
memset(lnk,-1,sizeof(lnk));
for(int i=1,a,b;i<=e;i++){
a=read();b=read();
if(b>m||a>n) continue;
divide[a].push_back(b);
}
for(int i=1;i<=n;i++){
if(Find(i,i))//省掉了memset
ans++;
}
printf("%d\n",ans);
return 0;
}
下面是没有省略那个vis数组的代码
#include <bits/stdc++.h>
using namespace std;
const int maxn=1001;
int n,m,v,tot,last[maxn*2],link[maxn],vis[maxn];
struct near{
int to,pre;
}tree[maxn*maxn];
bool dfs(int x){
int i=last[x];
while (i){
int p=tree[i].to;
if (!vis[p]){
vis[p]=1;//这个点走过了
if (!link[p]||dfs(link[p])){//如果右边的某p点还没有人匹配 或者 他已经有人匹配 但是与它匹配的那个点能找到其它的匹配(也就是不和p匹配)
link[p]=x;//这里才脱单!
return 1;
}
}
i=tree[i].pre;
}
return 0;
}
int main(){
scanf("%d %d %d",&n,&m,&v);
int x,y;
for (int i=1;i<=v;i++){
scanf("%d %d",&x,&y);
if (x>n||y>m) continue;
tree[++tot]=(near){y,last[x]};//加一边就行,以左边为基准进行匹配
last[x]=tot;
}
int ans=0;
for (int i=1;i<=n;i++){
memset(vis,0,sizeof(vis));//这个是防止点重复走的,保证每次最多遍历一遍,并不代表匹配上了(什么脱单了之类的)
if (dfs(i)) ans++;//看到目前为止这个点能不配对上
}
printf("%d\n",ans);
return 0;
}