UVA - 1045 The Great Wall Game(二分图最佳完美匹配)

题目大意:给出棋盘上的N个点的位置。现在问将这些点排成一行或者一列,或者对角线的最小移动步数(每个点都只能上下左右移动,一次移动一个)

解题思路:暴力+二分图最佳完美匹配

#include <cstdio>
#include <cstring>

#define N 20
#define INF 0x3f3f3f3f
#define abs(x) ((x) > 0 ? (x) : (-(x)))
#define max(a,b)((a)>(b)?(a):(b))
#define min(a,b)((a)<(b)?(a):(b))

struct Node {
    int x, y;
}node[N];

int w[N][N], left[N], Lx[N], Ly[N], slack[N];
int ans, n;
bool S[N], T[N]; 

bool match(int i) {
    S[i] = true;
    for (int j = 1; j <= n; j++) {
        if (Lx[i] + Ly[j] == w[i][j] && !T[j]) {
            T[j] = true;
            if (!left[j] || match(left[j])) {
                left[j] = i;
                return true;
            }
        }
        else slack[j] = min(slack[j], Lx[i] + Ly[j] - w[i][j]);
    }
    return false;
}

void update() {
    int a = 1 << 30;
    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[i] - 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++) {
        left[i] = Ly[i] = 0;
        Lx[i] = -INF;
        for (int j = 1; j <= n; j++)
            Lx[i] = max(Lx[i], w[i][j]);
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= n; j++)
            slack[j] = INF;
        for (;;) {
            for (int j = 1; j <= n; j++) S[j] = T[j] = 0;
            if (match(i)) break;
            int a = INF;
            for (int j = 1; j <= n; j++)
                if (!T[j])
                    a = min(a, slack[j]);
            for (int j =  1; j <= n; j++) {
                if (S[j]) Lx[j] -= a;
                if (T[j]) Ly[j] += a;   
            }
        }
    }
    int t = 0;
    for (int i = 1; i <= n; i++)
        t += Lx[i] + Ly[i];
    ans = max(ans, t);
}

void init() {
    for (int i = 1; i <= n; i++)
        scanf("%d%d", &node[i].x, &node[i].y);
}

int cas = 1;
void solve() {
    ans = -INF;
    for (int j = 1; j <= n; j++) {
        for (int k = 1; k <= n; k++) {
            for (int i = 1; i <= n; i++) {
                w[i][k] =  abs(node[i].x - j) + abs(node[i].y - k);
                w[i][k] = -w[i][k];
            }
        }
        KM();
    }

    for (int j = 1; j <= n; j++) {
        for (int k = 1; k <= n; k++) 
            for (int i = 1; i <= n; i++) {
                w[i][k] = abs(node[i].x - k) + abs(node[i].y - j);
                w[i][k] = -w[i][k];
            }
        KM();
    }

    for (int j = 1; j <= n; j++)
        for (int i = 1; i <= n; i++) {
            w[i][j] = abs(node[i].x - j) + abs(node[i].y - j);
            w[i][j] = -w[i][j];
        }
    KM();

    for (int j = 1; j <= n; j++)
        for (int i = 1; i <= n; i++) {
            w[i][j] = abs(node[i].x - j) + abs(node[i].y - (n - j + 1));
            w[i][j] = -w[i][j];
        }
    KM();
    printf("Board %d: %d moves required.\n\n", cas++, abs(ans));
}

int main() {
    while (scanf("%d", &n) != EOF && n) {
        init();
        solve();
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值