HDU 2255奔小康赚大钱

HDU 2255奔小康赚大钱

题目链接：HDU 2255

题意:

村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住
每家必须分配到一间房子且只能得到一间房子
由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,
比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.
每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量)，
接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。


分析:

将n个村民分配到n个房子，同时使得获得的金钱最多，即为完备匹配下的最优匹配。

直接上模板

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

const int maxn=1005;
int g[maxn][maxn];
bool visitx[maxn], visity[maxn];
int lx[maxn], ly[maxn];
int slack[maxn];
int match[maxn];
int n,m;

bool Hungary(int u) //匈牙利算法
{
    visitx[u] = true;
    for(int i = 1; i <= n; ++i)
    {
        if(visity[i])
            continue;
        if(lx[u] + ly[i] == g[u][i])
        {
            visity[i] = true;
            if(match[i] == -1 || Hungary(match[i]))
            {
                match[i] = u;
                return true;
            }
        }
        else //不在相等子图
            slack[i] = min(slack[i], lx[u] + ly[i] - g[u][i]);
    }
    return false;
}

void KM()
{
    int temp;
	memset(match, -1, sizeof(match));
    memset(lx, 0, sizeof(lx)); //初始化顶标
    memset(ly, 0, sizeof(ly)); //ly[i]为0

    for(int i = 1; i <= n; ++i) //lx[i]为权值最大的边
        for(int j = 0; j <= m; ++j)
            lx[i] = max(lx[i], g[i][j]);

    for(int i = 1; i <= n; i++) //对n个点匹配
    {
        for(int j = 1; j <= n; j++)
            slack[j] = INT_MAX;
        while(1)
        {
            memset(visitx, false, sizeof(visitx));
            memset(visity, false, sizeof(visity));
            if(Hungary(i)) //匹配成功
                break;
            else //匹配失败，找最小值
            {
                temp = INT_MAX;
                for(int j = 1; j <= n; j++)
                    if(!visity[j])
                        if(temp > slack[j])
                            temp = slack[j];
                for(int j = 1; j <= n; j++) //更新顶标
                {
                    if(visitx[j])
                        lx[j] -= temp;
                }
                for(int j=1;j<=m;++j)
                {
                   if(visity[j])
                        ly[j] += temp;
                    else
                        slack[j] -= temp;
                }
            }
        }
    }
	int ans = 0;
    for(int i = 1; i <= m; ++i) //权值相加
        ans += g[match[i]][i];
    printf("%d\n", ans);
}

int main()
{
	// freopen("in.txt","r",stdin);
    while(scanf("%d", &n) != EOF)
    {
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
				scanf("%d",&g[i][j]);
		}
		m=n;
		KM();
    }
    return 0;
}