简单二分图判定与最大匹配算法

前言

二分图,又称二部图,英文名叫 Bipartite graph。
通俗一点就是一个无向图如果能划成两个非空点集,且两部分内部没有边,则这是一张二分图。

如果从颜色的角度来说,就是把节点染成黑色/白色,并且使得没有相邻的节点颜色相同。

一张二分图:
在这里插入图片描述

二分图判定

奇环:长度为奇数的环

一张无向图 G G G为二分图的充要条件是:图中不存在奇环
证明:
必要性显然,因为奇环显然不是二分图。
充分性:只需证明没有奇环的图就是二分图
考虑无向图 G G G没有奇环且连通,若不连通可以对每个连通块单独考虑。
从节点 x x x开始染色,把 x x x染成白色,找到节点 x x x到图中所有点的所有路径,对于一条从 x x x y y y的路径 ( u 0 = x , u 1 , u 2 , . . . , u k = y ) (u_0=x,u_1,u_2,...,u_k=y) (u0=x,u1,u2,...,uk=y),我们把 u i ∈ 奇数 u_{i\in 奇数} ui奇数染黑,把 u i ∈ 偶数 u_{i\in 偶数} ui偶数染白。

  • 如果这个操作中一个点只被染成了一种颜色,那我们就构造出了一张二分图的划分
  • 如果存在一个点被染成了两种颜色,那么说明它在某两条路径到 x x x的路径中分别位于奇数位置和偶数位置,这时候我们把这两条路径拼起来,就得到了一个奇环,矛盾。

QED.

但是二分图判定的代码并不是找奇环的,而是dfs/bfs染色,因为图中环的个数是指数级的,找奇环是很亏的。

dfs/bfs判奇环复杂度为 O ( n + m ) O(n+m) O(n+m)

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

由于作者实力不济,本文没有匈牙利算法的证明。

匈牙利算法大概的过程是:

  • 对每个左部点 i i i执行一遍dfs(i),保证每一次期间,执行一个右部点只会被访问一次
  • dfs(u)的过程是:
    • 遍历左部点 u u u的未被访问的右部点 v v v
    • 标记 v v v
    • 检查 v v v是否有match,如果没有抢,如果有,就dfs(match[v])看看能不能抢

换句话说由攻来追受,由于攻比受少,钦定左部点为攻。对于攻 u u u,如果其的受 v v v还没有攻,那就直接匹配。否则找到 v v v对应的攻 m a t c h [ v ] match[v] match[v],检查 m a t c h [ v ] match[v] match[v]的备胎受。如果可以换备胎,那就更换,然后 u , v u,v u,v匹配。如果不行,那就检查 u u u的下一个受 v ′ v' v

时间复杂度 O ( l ⋅ m + r ) O(l\cdot m+r) O(lm+r),其中 l l l是左部点数量, m m m是边数, r r r是右部点数量。

实现

#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
const int N=500;
vector<int> a[N+5];
int match[N+5];
bool vis[N+5];
bool dfs(int u) {
	for(auto&v:a[u])
		if(!vis[v]) {
			vis[v]=1;
			if(!match[v]||dfs(match[v]))
				return (match[v]=u);
		}
	return 0;
}
int main() {
	int l,r,m;
	cin>>l>>r>>m;
	for(int i=1,u,v;i<=m;i++){
		cin>>u>>v;
		a[u].push_back(v);
	}
	int ans=0;
	for(int i=1;i<=l;i++,memset(vis,0,sizeof vis))
		ans+=dfs(i);
	cout<<ans;
	
}

输出方案

#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
const int N=100;
vector<int>a[N+5];
int match[N+5];
bool vis[N+5];
bool dfs(int u){
	for(auto&v:a[u])
		if(!vis[v]){
			vis[v]=1;
			if(!match[v]||dfs(match[v]))
				return (match[v]=u);
		}
	return 0;
}
int main(){
	int l,r;
	cin>>l>>r;
	r-=l;
	int u,v;
	cin>>u>>v;
	while(~v){
		v-=l;
		a[u].push_back(v);
		cin>>u>>v;
	}
	int ans=0;
	for(int i=1;i<=l;i++,memset(vis,0,sizeof vis))
		ans+=dfs(i);
	cout<<ans<<endl;
	for(int i=1;i<=r;i++) 
		if(match[i])
			cout<<match[i]<<' '<<i+l<<endl;
}

后记

于是皆大欢喜。

  • 22
    点赞
  • 23
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值