前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >Uva – 11383 – Golden Tiger Claw

Uva – 11383 – Golden Tiger Claw

作者头像
全栈程序员站长
发布2021-06-22 19:52:24
2410
发布2021-06-22 19:52:24
举报

题意:一个N*N的矩阵,第i行第j列的元素大小为w[i][j],每行求一个数row[i],每列求一个数col[j],使得row[i] + col[j] >= w[i][j],且所有的row[]与所有的col[]和总和最小( N <= 500, 其它输入数为正整数且 <= 100)。

题目链接:http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2378

——>>row[i] + col[j] >= w[i][j],这个恰恰是二分图最佳完美匹配的一个式子,所以,以行row为X结点,以列col为Y结点,权值即为对应元素w[i][j]的值建图,跑一次KM就好。

另外发现:用scanf(“%d”, &N) == 1比用~scanf(“%d”, &N)快了3ms。。。

代码语言:javascript
复制
#include <cstdio>
#include <algorithm>

using namespace std;

const int maxn = 500 + 10;
const int INF = 0x3f3f3f3f;

int N, w[maxn][maxn], lx[maxn], ly[maxn], fa[maxn];
bool S[maxn], T[maxn];

bool match(int i){
    S[i] = 1;
    for(int j = 1; j <= N; j++) if(lx[i] + ly[j] == w[i][j] && !T[j]){
        T[j] = 1;
        if(!fa[j] || match(fa[j])){
            fa[j] = i;
            return 1;
        }
    }
    return 0;
}

void update(){
    int a = INF;
    for(int i = 1; i <= N; i++) if(S[i])
        for(int j = 1; j <= N; j++) if(!T[j])
            a = min(a, lx[i] + ly[j] - w[i][j]);
    for(int i = 1; i <= N; i++){
        if(S[i]) lx[i] -= a;
        if(T[i]) ly[i] += a;
    }
}

void KM(){
    for(int i = 1; i <= N; i++){
        fa[i] = lx[i] = ly[i] = 0;
        for(int j = 1; j <= N; j++) lx[i] = max(lx[i], w[i][j]);
    }
    for(int i = 1; i <= N; i++)
        while(1){
            for(int j = 1; j <= N; j++) S[j] = T[j] = 0;
            if(match(i)) break;
            else update();
        }
}

void read(){
    for(int i = 1; i <= N; i++)
        for(int j = 1; j <= N; j++) scanf("%d", &w[i][j]);
}

void solve(){
    for(int i = 1; i < N; i++) printf("%d ", lx[i]); printf("%d\n", lx[N]);
    for(int i = 1; i < N; i++) printf("%d ", ly[i]); printf("%d\n", ly[N]);
    int sum = 0;
    for(int i = 1; i <= N; i++) sum += lx[i] + ly[i];
    printf("%d\n", sum);
}

int main()
{
    while(scanf("%d", &N) == 1){
        read();
        KM();
        solve();
    }
    return 0;
}

发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/101184.html原文链接:https://javaforall.cn

本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档