[KM算法 O(n^ 3) ] 带权二分图最大匹配

最新推荐文章于 2024-01-15 21:06:32 发布

原创最新推荐文章于 2024-01-15 21:06:32 发布 · 345 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

模板同时被 2 个专栏收录

47 篇文章

订阅专栏

图论

7 篇文章

订阅专栏

本文探讨了如何在奔小康赚大钱的过程中，利用KM算法解决O(n^3)复杂度的带权二分图最大匹配问题。通过算法解析，帮助读者理解并应用该算法来优化资源分配。

奔小康赚大钱

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 18257 Accepted Submission(s): 7702

Problem Description

传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革：重新分配房子。
这可是一件大事,关系到人民的住房问题啊。村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住（如果有老百姓没房子住的话，容易引起不安定因素），每家必须分配到一间房子且只能得到一间房子。
另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.(当然是在他们的经济范围内).现在这个问题就是村领导怎样分配房子才能使收入最大.(村民即使有钱购买一间房子但不一定能买到,要看村领导分配的).

Input

输入数据包含多组测试用例，每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量)，接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。

Output

请对每组数据输出最大的收入值，每组的输出占一行。

Sample Input

2 100 10 15 23

Sample Output

123

Source

HDOJ 2008 Summer Exercise（4）- Buffet Dinner

#include <bits/stdc++.h>
using namespace std;

const int mn = 305, inf = 2e9;

int val[mn][mn], match[mn], aex[mn], bex[mn], slack[mn];
bool avis[mn], bvis[mn];

bool dfs(int a, int n)
{
    avis[a] = 1;
    for (int b = 1; b <= n; b++)
    {
        if (bvis[b]) continue;
        int gap = aex[a] + bex[b] - val[a][b];
        if (gap == 0)
        {
            bvis[b] = 1;
            if (match[b] == -1 || dfs(match[b], n))
            {
                match[b] = a;
                return 1;
            }
        }
        else
            slack[b] = min(slack[b], gap);
    }
    return 0;
}

int km(int n)
{
    memset(match, -1, sizeof match);
    memset(bex, 0, sizeof bex);	/// 第二列的期望
    
    for (int i = 1; i <= n; i++) // 第一列期望为边权最大值
    {
        aex[i] = val[i][1];
        for (int j = 2; j <= n; j++)
            aex[i] = max(aex[i], val[i][j]);
    }
    
    for (int i = 1; i <= n; i++)
    {
        for (int i = 1; i <= n; i++) slack[i] = inf;	// 松弛变量
        while (1)
        {
            memset(avis, 0, sizeof avis);
            memset(bvis, 0, sizeof bvis);
            
            if (dfs(i, n)) break;	// 可匹配
            
            int d = inf;	// 需要降低的期望的最小值
            for (int j = 1; j <= n; j++)
                if (!bvis[j]) d = min(d, slack[j]);
            
            for (int j = 1; j <= n; j++)
            {
                if (avis[j]) aex[j] -= d;
                if (bvis[j]) bex[j] += d;
                else slack[j] -= d;
            }
        }
    }
    int res = 0;	// 最大匹配值
    for (int i = 1; i <= n; i++)
        res += val[match[i]][i];
    return res;
}

int main()
{
    int n;
    while (~scanf("%d", &n))
    {
        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                scanf("%d", &val[i][j]);
        printf("%d\n", km(n));	/// KM 带权最大图匹配
    }
    return 0;
}