染色法
用途:判断一个图是否为二分图
二分图 :所有的点可以分为两个集合V1、V2,图中任意一条边的两个顶点分别在不同的集合,即任意一个集合的内部不能有边。
图论中一个很重要的性质:一个图是二分图,当且仅当图中不含奇数环(边数为奇数的环)
染色法:
时间复杂度:深度优先遍历O(m+n)
思路:将一个图按上述二分图的定义划分为两个集合的过程,可以理解为一个染色的过程,最终我们会将一个图中的所有点染成两种颜色,代表两个集合,如果有邻接点颜色相同(在一个集合),则染色失败。而染色失败的情况一定是因为图中含有奇数环。
(图中不含奇数环,则染色过程中一定没有矛盾)
举个简单的例子:我们用数字”1“和“2”代表两种颜色。显然,只要有一个点的颜色确定了,那么其它点的颜色都可以确定(相邻点不能染相同颜色,因为集合内部不能用有边)
典例:AcWing 860.染色法判定二分图
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1e5 + 10,M = 2e5 + 10;
int h[N],e[M],ne[M],idx; //用邻接表来存储
int n,m;
int color[N]; //标记所染颜色
void add(int a,int b)
{
e[idx] = b,ne[idx] = h[a],h[a] = idx ++ ;
}
bool dfs(int u,int c)
{
color[u] = c; //当前点u的颜色为c
for(int i = h[u]; i != -1; i = ne[i])
{
int j = e[i]; //用j来存储当前点的编号
if(!color[j]) //当前点没被染色
{
//染成另外一种颜色(3-c可以把1变2,2变1)
if(!dfs(j,3-c)) return false;
}
else if(color[j] == c) return false;
}
return true;
}
int main()
{
cin >> n >> m;
memset(h,-1,sizeof h);
while(m -- )
{
int a,b;
cin >> a >> b;
add(a,b),add(b,a); //无向图
}
bool flag = true; //标记染色过程中是否矛盾
for(int i = 1; i <=n; i ++ )
if(!color[i]) //i点没有被染过
{
if(!dfs(i,1)) //如果dfs返回false,我们就认为有矛盾发生
{
flag = false;
break;
}
}
if(flag) puts("Yes");
else puts("No");
return 0;
}
匈牙利算法
理论时间复杂度O(mn),实际运行时间一般远小于该值。
理论时间复杂度高,但一般实际运行效果很好的算法还有很多,如:最大流算法。
用途:给定一个二分图,求最大匹配。
优势:用比较快的时间复杂度,求得左右集合中匹配成功的数量
成功的匹配:不存在两条边共用同一顶点
最大匹配:在成功匹配的基础上边数最大
举个有趣的“栗子”
我们可以把匈牙利算法的匹配过程类比成“找对象”,如下图,左边集合表示Boys,右边集合表示Girls,我们看按照匈牙利算法的思想,可以凑成几对儿情侣。
男生A看上了女生F,女生F恰巧没有男朋友,于是,他们成功牵手;<A,F>
男生B看上了女生E,女生E恰巧没有男朋友,于是他们成功走到了一起;<B,E>
男生C看上了女生F,但是女生F已经和A在一起了,男生C得知,去找男生A,男生A是个软蛋(同时也是个渣男,他其实早就对女生H有想法),于是男生A抛弃了女生F和H在一起了,男生C也成功得到了女生F;<A,H>,<C,F>
男生D和女生G都没有对象,他们也一见钟情,幸福地走到一起。<D,G>
综上,图中有四对情侣,也就是有四组匹配。
匈牙利算法就是这样,只有尝试了各种各样的方法,都达不到想要的“匹配”,才会放弃。
典例:AcWing 861.二分图的最大匹配
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 510,M = 1e5 + 10;
int n1,n2,m;
int h[N],e[M],ne[M],idx;
int match[N]; //右边点对应的点(存储每个妹子和哪个男生在一起)
bool st[N]; //判重,每次不要重复搜一个点
void add(int a,int b)
{
e[idx] = b,ne[idx] = h[a],h[a] = idx ++ ;
}
//判断是否能够匹配
bool find(int x)
{
for(int i = h[x]; i != -1; i = ne[i]) // 遍历 x 男孩喜欢所有的女孩
{
int j = e[i]; //用j表示“女生”集合中点的编号
if(!st[j])
{
st[j] = true;
//当前的女生还没有男朋友,或女孩的男朋友可以和其他喜欢的女生在一起
if(match[j] == 0 || find(match[j]))
{
match[j] = x;
return true;
}
}
}
return false;
}
int main()
{
cin >> n1 >> n2 >> m;
memset(h,-1,sizeof h);
while(m -- )
{
int a,b;
cin >> a >> b;
add(a,b);
}
int res = 0; //匹配的数量
for(int i = 1; i <= n1; i ++ ) //从前往后分析每个男生该找哪个女生
{
memset(st,false,sizeof st); //把所有右边集合清空(保证每个妹子考虑一遍)
if(find(i)) res ++ ;
}
cout << res << endl;
return 0;
}