ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

匈牙利算法求最大匹配(HDU-4185 Oil Skimming)

2018-10-02 19:34:39  阅读:201  来源: 互联网

标签:


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

大佬博客:https://blog.csdn.net/cillyb/article/details/55511666

模板:

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

int hunyary()
{
    int sum = 0;
    memset(link, 0, sizeof(link));
    for(int i = 1; i <= num; i++)
    {
        memset(vis, 0, 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 = 605;
int n,k,cnt = 0,num;
char mp[maxn][maxn];
int temp[maxn][maxn],cp[maxn][maxn];
int link[maxn],vis[maxn];

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

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

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


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

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

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

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

 

标签:
来源: https://www.cnblogs.com/sykline/p/9737800.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有