二分图——染色法与匈牙利算法

二分图：n阶无向图G是二部图当且仅当G中无奇圈。

染色法

思路：由二分图定义， 对于图中每条边，其端点位于不同集合中，因此，对于图中每个连通快，可用dfs为其中每个点打上标记， 如果存在矛盾的点，则该图不为二分图。时间复杂度O(n+M).

#include <iostream>
#include <cstring>

using namespace std;

const int N = 100010, M = 200010;
int e[M], ne[M], h[N], idx;
int color[N];
int n, m;

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

bool dfs(int now, int c)
{
    color[now] = c;
    for(int i = h[now]; i != -1; i = ne[i])
    {
        int j = e[i];
        if(color[j] == c) return false; 
        if(!color[j] && !dfs(j, -c)) return false;   
    }  
    return true;
}

int main()
{
    memset(h, -1, sizeof h);
    cin >> n >> m;
    
    int a, b;
    for(int i = 0; i < m; i++)
    {
        cin >> a >> b;
        add(a, b);
        add(b, a);
    }
    
    bool flag = true;
    for(int i = 1; i <= n; i++)
    {
        if(!color[i] && !dfs(i, 1)) 
        {
            flag = false;
            break;
        }
    }
    
    if(flag) cout << "Yes";
    else cout << "No";
    
    return 0;
}

匈牙利算法

思路：
对于二分图的 v1， v2点集， 欲求其最多匹配边数， 可用一个邻接表存v1各点所能访问的v2的点， 另用一个数组存v2各点当前匹配到的v1的点。遍历每个v1的点，如果寻找到的v2集合中的点未被占用，则使用，如已被占用，则看占用该点的v1集合中的点能否更改。时间复杂度O(nm).
注意：必须使用一个数组记录v2点是否被访问，否则可能出现多个v1中的点“争抢”v2中某点的情况，造成死循环。

#include <iostream>
#include <cstring>
using namespace std;
const int N = 510, M = 100010;
int h[N], e[M], ne[M], idx; // 用链表记录左半部分各点所能到达的点
int match[N]; // 记录右半部分各点所匹配到左半部分的点
bool vis[N]; // 用于判断右边点是否被访问，反正某点被重复改变出现死循环
int n1, n2, m;

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])
    {
        int j = e[i];
        if(!vis[j])
        {
            vis[j] = true;
            if(!match[j] || find(match[j]))
            {
                match[j] = x;
                return true;
            }
        }
    }
    return false;
}

int main()
{
    memset(h, -1, sizeof h);
    cin >> n1 >> n2 >> m;

    int l, r;
    for (int i = 0; i < m; i++)
    {
        cin >> l >> r;
        add(l, r);
    }
    
    int ans = 0;

    for(int i = 1; i <= n1; i++)
    {
        memset(vis, false, sizeof vis);
        if(find(i)) ans++;
    }

    cout << ans;
    return 0;
}