匈牙利算法求最大匹配（HDU-4185 Oil Skimming）

如下图：要求最多可以凑成多少对对象

大佬博客：

https://blog.csdn.net/cillyb/article/details/55511666

https://blog.csdn.net/denghecsdn/article/details/77619308

https://www.cnblogs.com/wangjunyan/p/5563154.html

模板：

int link[maxn],vis[maxn];
bool dfs(int x)
{
    for(int i = ; i <= num; i++)
    {
        if(!vis[i] && cp[x][i])
        {
            vis[i] = ;
            if(link[i] ==  || dfs(link[i]))
            {
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}

int hunyary()
{
    int sum = ;
    memset(link, , sizeof(link));
    for(int i = ; i <= num; i++)
    {
        memset(vis, , sizeof(vis));
        if(dfs(i))
            sum++;
    }
    return sum;
}

最后的sum值就是最大的匹配值。

附练习题：HDU-4185 Oil Skimming

水平和竖直的两个‘#’，可以连成一条边，求这样的边最多有多少条。

主要是自己如何去建图：

1. 遍历数组mp对‘#’进行标号；
2. 再次遍历遇到‘#’，看他上下左右是否有‘#’，有的话在cp数组中进行连线；
3. 再利用匈牙利算法求最大匹配。

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <queue>
#include <vector>
#include <algorithm>
#define FRE() freopen("in.txt","r",stdin)
#define INF 0x3f3f3f3f

using namespace std;
typedef long long ll;
const int maxn = ;
int n,k,cnt = ,num;
char mp[maxn][maxn];
int temp[maxn][maxn],cp[maxn][maxn];
int link[maxn],vis[maxn];

inline void init()
{
    memset(temp, , sizeof(temp));
    memset(cp, , sizeof(cp));
}

bool dfs(int x)
{
    for(int i = ; i <= num; i++)
    {
        if(!vis[i] && cp[x][i])
        {
            vis[i] = ;
            if(link[i] ==  || dfs(link[i]))
            {
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}

int hunyary()
{
    int sum = ;
    memset(link, , sizeof(link));
    for(int i = ; i <= num; i++)
    {
        memset(vis, , sizeof(vis));
        if(dfs(i))
            sum++;
    }
    return sum;
}

int main()
{
    //FRE();
    scanf("%d",&k);
    while(k--)
    {
        init();
        scanf("%d",&n);
        for(int i = ; i < n; i++)
            scanf("%s",mp[i]);

        num = ;
        for(int i = ; i < n; i++)
        {
            for(int j = ; j < n; j++)
            {
                if(mp[i][j] == '#')
                    temp[i][j] = ++num;
            }
        }

        for(int i = ; i < n; i++)
        {
            for(int j = ; j < n; j++)
            {
                if(mp[i][j] == '#')
                {
                    if(i -  >=  && mp[i-][j] == '#') cp[temp[i][j]][temp[i-][j]] = ;
                    if(i +  < n && mp[i+][j] == '#') cp[temp[i][j]][temp[i+][j]] = ;
                    if(j -  >=  && mp[i][j-] == '#') cp[temp[i][j]][temp[i][j-]] = ;
                    if(j +  < n && mp[i][j+] == '#') cp[temp[i][j]][temp[i][j+]] = ;
                }

            }
        }
        int res = hunyary();
        printf("Case %d: %d\n",++cnt, res / );
    }
    return ;
}