2018"百度之星"程序设计大赛 - 资格赛 - 整数规划（最小权匹配-KM算法）

最新推荐文章于 2020-07-21 17:10:27 发布

原创最新推荐文章于 2020-07-21 17:10:27 发布 · 801 阅读

0 ·

CC 4.0 BY-SA版权

HDU 同时被 2 个专栏收录

170 篇文章

订阅专栏

二分图

7 篇文章

订阅专栏

本文详细介绍了使用匈牙利算法解决二分图最大权匹配问题的过程，并提供了完整的C++实现代码。通过调整权值，该算法能够求得匹配的最大权值总和。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://acm.hdu.edu.cn/showproblem.php?pid=6346

#include <stdio.h>
#include <iostream>
#include <algorithm>
using namespace std;
#define  LL long long
const LL N = 222;
const LL inf = 0x3f3f3f3f3f3f3f3f;
LL n;
LL val[N][N];
LL lx[N],ly[N];
LL linky[N];
LL pre[N];
bool vis[N];
bool visx[N],visy[N];
LL slack[N];

void bfs(LL k){
    LL px, py = 0,yy = 0, d;
    memset(pre, 0, sizeof(LL) * (n+2));
    memset(slack, inf, sizeof(LL) * (n+2));
    linky[py]=k;
    do{
        px = linky[py],d = inf, vis[py] = 1;
        for(LL i = 1; i <= n; i++)
            if(!vis[i]){
                if(slack[i] > lx[px] + ly[i] - val[px][i])
                    slack[i] = lx[px] + ly[i] -val[px][i], pre[i]=py;
                if(slack[i]<d) d=slack[i],yy=i;
            }
        for(LL i = 0; i <= n; i++)
            if(vis[i]) lx[linky[i]] -= d, ly[i] += d;
            else slack[i] -= d;
        py = yy;
    }while(linky[py]);
    while(py) linky[py] = linky[pre[py]] , py=pre[py];
}
LL KM(){
    memset(lx, 0, sizeof(LL)*(n+2));
    memset(ly, 0, sizeof(LL)*(n+2));
    memset(linky, 0, sizeof(LL)*(n+2));
    for(LL i = 1; i <= n; i++)
        memset(vis, 0, sizeof(bool)*(n+2)), bfs(i);
    LL ans = 0;
    for(LL i = 1; i <= n; ++i)
        ans += lx[i] + ly[i];
    return ans;
}
int main()
{
    LL T;
    scanf("%lld",&T);
    LL cas=0;
    while(T--){
        scanf("%lld",&n);
        for(LL i=1;i<=n;i++)
            for(LL j=1;j<=n;j++)
                scanf("%lld",&val[i][j]),val[i][j]=-val[i][j];
        printf("Case #%lld: %lld\n",++cas,-KM());
    }
    return 0;
}