（匹配）Antenna Placement --POJ --3020

匈牙利算法解析

最新推荐文章于 2019-11-08 22:12:50 发布

转载最新推荐文章于 2019-11-08 22:12:50 发布 · 80 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/YY56/p/4722163.html

本文介绍了一种使用匈牙利算法解决最大匹配问题的方法。通过分析矩阵中的元素关系，找到最优匹配路径，实现对特定问题的有效求解。文章提供了一个具体的实现案例，包括输入输出示例、核心代码片段及解释。

链接：

http://poj.org/problem?id=3020

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=82834#problem/H

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 510
#define INF 0x3f3f3f3f


// un是匹配左边的定点数， vn是匹配右边的定点数
int n, m, un, vn, used[N], p[N], hash[N][N], g[N][N];
char G[N][N];

//匈牙利算法， 从左边开始找增广路
int Find(int u)
{
    for(int j=0; j<vn; j++)
    {
        if(!used[j] && g[u][j])
        {
            used[j] = 1;
            if(p[j]==-1 || Find(p[j]))
            {
                p[j] = u;
                return true;
            }
        }
    }
    return false;
}


//最大匹配数
int hungary()
{
    int ans = 0;

    memset(p, -1, sizeof(p));
    for(int i=0; i<un; i++)
    {
        memset(used, 0, sizeof(used));
        if(Find(i)) ans++;
    }
    return ans;
}


int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        int i, j, tol=0;

        scanf("%d%d", &n, &m);
        memset(G, 0, sizeof(G));
        memset(hash, 0, sizeof(hash));
        memset(g, 0, sizeof(g));
        for(i=0; i<n; i++)
        {
            scanf("%s", G[i]);
            for(j=0; j<m; j++)
                if(G[i][j]=='*')
                hash[i][j]=tol++;
        }

       for(i=0; i<n; i++)
       for(j=0; j<m; j++)
       {
           if(G[i][j]=='*')
           {
               if(i>0 && G[i-1][j]=='*')     g[hash[i][j]][hash[i-1][j]]=1;
               if(i<n-1 && G[i+1][j]=='*')   g[hash[i][j]][hash[i+1][j]]=1;
               if(j>0 && G[i][j-1]=='*')     g[hash[i][j]][hash[i][j-1]]=1;
               if(j<m-1 && G[i][j+1]=='*')   g[hash[i][j]][hash[i][j+1]]=1;
           }
       }
        un = vn = tol;
        printf("%d\n", tol-hungary()/2);

    }
    return 0;
}