二分图算法总结(C/C++)

一. 二分图判定——染色法

1. 算法介绍

  如何判定一个图是否是二分图 ?
  解答 ：只需要我们从一个图的第一个点认为是颜色 1，然后利用深度优先搜索的方法不断去染色其他方块为不同颜色 2，依次循环如果检查到矛盾则不是二分图，否则二分图成立。

2. 关键分析

bool dfs(int u, int c){
	// 首先将此点染色为颜色 c 
    color[u] = c;
    //再逐个对相邻的点搜索
    //若没染色-检查染色为相反颜色是否成立，若已经染色-判断颜色是否矛盾
    for(int i = h[u]; ~i; i = ne[i] ){
        int j =e[i];
        if(!color[j]){ // 没染色-检查是否可以染色为相反颜色
            if(!dfs(j, 3 - c) ) return false;
        }
        else if(color[j] != 3 - c) return false; // 染色-检查是否矛盾
    }
    return true;
}

3. 完整代码

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
using namespace std;
const int N = 1e5 + 10;
int h[N], e[N << 1], ne[N << 1], idx = 0;
int n, m, color[N];
void add(int x, int y){
    e[idx] = y, ne[idx] = h[x], h[x] = idx ++;
}

bool dfs(int u, int c){
    color[u] = c;
    for(int i = h[u]; ~i; i = ne[i] ){
        int j =e[i];
        if(!color[j]){
            if(!dfs(j, 3 - c) ) return false;
        }
        else if(color[j] != 3 - c) return false;
    }
    return true;
}

int main(){
    ios::sync_with_stdio(false); cin.tie(0);
    memset(h, -1,sizeof(h));
    cin >> n >> m;
    for(int i = 0; i < m; i ++ ){
        int x, y;
        cin >> x >> y;
        add(x, y), add(y, x);
    }

    bool flag = true;
    for(int i = 1; i <= n; i ++ ){
        if(!color[i]){
            if(!dfs(i, 1)){
                flag = false;
                break;
            }
        }
    }
    if(flag) puts("Yes");
    else puts("No");
    return 0;

}

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

1. 算法介绍

  如何找到最大匹配 (二分图匹配：同一个点最多只能匹配一条边) ？
  解答 ：通俗的说就是两个集合A 和 B，我们每次从集合A开始选择一个元素出来，再将其对应的集合B当中的元素进行逐个匹配，如果刚好这个元素没有被匹配，或者这个B集合中元素虽然已经被匹配但可以被让出来我们都认为其成立。

2. 关键分析

  假如选择元素已经被选择，如何判断B集合当中的这个元素？
  解答：找到这个B集合的元素 y 对应的 A集合当中的元素 x1, x2, x3 依次进行下面的方案操作找到第一个满足的情况即可
 方案 ：请求 x 能否重新选择一个新元素—— x 选择下一个B集合元素 yy(如果成立直接返回， 如果不成立那么申请 yy对应的 xx 能否重新选择) …不断重复这个过程即可
 核心：没有碰撞直接选择，若有碰撞不断循环走到最底层的匹配开始进行让位的操作直到从最初匹配的集合A元素到上一个成功匹配集合A元素都无法避让出集合B当中的位置为止。

/* 
1. 这里的st[]每次会被初始赋值为false ，然后依次假设集合A元素 x 去选择看能否成立
其实st[] 还有一个优点,就是假如你进行前面元素让位发现有些B集合元素b2为true 
假设下一个元素为b2时，为什么可以跳过这个过程呢？
因为b2是在之前遍历中已经被确定是无法避让位置的元素,所以当你再次去遍历遇到b2的时候是可以
直接跳过这个过程的
2. match[j] 表示对应于集合B当中j元素的对应元素(对应于集合A)
*/
bool find(int x){
    for(int i = h[x]; ~ i; i = ne[i] ){
    	// 按顺序选择对应集合B当中的元素，找到第一个能让位或者空位的元素
        int j = w[i];
        if(!st[j] ){ // st[] 表示此元素还未被选择，假设这个元素选择了集合B当中 J 判断是否成立
            st[j] = true;
            if( match[j] == 0 || find(match[j]) ){ // 为空或者可以避让出位置就成立
                match[j] = x;
                return true;
            }
        }
    }
    return false;
}

3. 完整代码

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
using namespace std;
const int N = 1e5 + 10;
int h[N], e[N << 1], ne[N << 1], idx = 0;
int n1, n2, m, match[N] ,res = 0;
bool st[N];

void add(int x, int y){
    e[idx] = y, ne[idx] = h[x], h[x] = idx ++;
}

bool find(int u){
    for(int i = h[u]; ~ i; i = ne[i] ){
        int j = e[i];
        if(!st[j]){
            st[j] = true;
            if(match[j] == 0 || find(match[j]) ){
                match[j] = u;
                return true;
            }
        }
    }
    return false;
}

int main(){
    memset(h, -1,sizeof(h));
    
    cin >> n1 >> n2 >> m;
    for(int i = 0; i < m; i ++ ){
        int x, y;
        cin >> x >> y;
        add(x, y);
    }
    
    for(int i = 1;i <= n1; i ++){
        memset(st, false, sizeof(st) );
        if(find(i)) res ++ ;
    }
    
    cout << res << endl;
    
}