- 染色法 O ( n + m ) O(n +m) O(n+m),二分图的判定
- 匈牙利算法 O ( m n ) O(mn) O(mn),实际运行时间一般远小于 O ( m n ) O(mn) O(mn)
- 最大流
什么是二分图
二分图:如果一个图的 N (N >= 2)个节点,可以分成 A, B两个非空集合,其中 A ⋂ B = ∅ A\bigcap B = \varnothing A⋂B=∅,并且在同一集合的点之间都没有边相连,那么称这个图为二分图。AB为图的左部和右部
染色法
首先说明一个性质:一个图是二分图,当且仅当图中不含奇数环(环边数是奇数)
读者可自行画圈证明
所以,如果我们对一个图进行染色,对一个点进行染白色,其相邻的点染成黑色,即对某个点染相反的颜色
如果一个图成功染色,说明它不含奇数环,说明这个图是二分图
由于某个点染色成功不代表整个图就是二分图,因此只有某个点染色失败才能立刻break/return
// Problem: 染色法判定二分图
// Contest: AcWing
// URL: https://www.acwing.com/problem/content/862/
// Memory Limit: 64 MB
// Time Limit: 1000 ms
// Code by: ING__
//
// Powered by CP Editor (https://cpeditor.org)
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 1e5 + 10, M = 2 * N;
int n, m;
int h[N], e[M], ne[M], idx;
int color[N];
void add(int u, int v){
e[idx] = v;
ne[idx] = h[u];
h[u] = idx++;
}
bool dfs(int u, int c){
color[u] = c;
for(int i = h[u]; i != -1; i = ne[i]){
int j = e[i];
if(!color[j]){ // 那个边连过来的这个点如果没有被染过
if(!dfs(j, 3 - c)) return 0; // 我们就染一下这个点看看行不行
}
else if(color[j] == c) return 0;
}
return 1;
}
int main(){
memset(h, -1, sizeof(h));
cin >> n >> m;
for(int i = 1; i <= m; i++){
int x, y;
cin >> x >> y;
add(x, y);
add(y, x);
}
bool t = 1;
for(int i = 1; i <= n; i++){
if(!color[i]){
if(!dfs(i, 1)){
t = 0;
break;
}
}
}
if(t){
cout << "Yes";
}
else{
cout << "No";
}
}
二分图最大匹配
任意两条边都没有公共断点的边的集合被称为图的一组匹配。在二分图中,包含边数最多的一组匹配被称为二分图的最大匹配
二分图的一组匹配 S 是最大匹配,当且仅当图中不存在 S 的增广路
关于匈牙利算法解释我解释不大清楚,看一下这个吧https://zhuanlan.zhihu.com/p/208596378
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 510;
const int M = 1e5 + 10;
int n1, n2;
int m;
int match[N];
bool vis[N];
int h[N], e[M], ne[M], idx;
void add(int u, int v){
e[idx] = v;
ne[idx] = h[u];
h[u] = idx++;
}
bool dfs(int x){
for(int i = h[x]; i != -1; i = ne[i]){
int j = e[i];
if(!vis[j]){ // 之前没有考虑过
vis[j] = 1;
if(!match[j] || dfs(match[j])){ // 可以为j匹配好的匹配另一个
match[j] = x;
return 1;
}
}
}
return 0;
}
int main(){
memset(h, -1, sizeof(h));
cin >> n1 >> n2 >> m;
while(m--){
int x, y;
cin >> x >> y;
add(x, y);
}
int ans = 0;
for(int i = 1; i <= n1; i++){
memset(vis, 0, sizeof(vis));
if(dfs(i)) ans ++;
}
cout << ans;
return 0;
}