uva 1349 - Optimal Bus Route Design(完美匹配)
2015-09-23 22:27
471 查看
题目链接:uva 1349 - Optimal Bus Route Design
权值变为负的即便为求最小值。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn = 105;
const int inf = 0x3f3f3f3f;
int N, L[maxn], Lx[maxn], Ly[maxn], W[maxn][maxn], slack[maxn];
bool S[maxn], T[maxn];
int match (int u) {
S[u] = true;
for (int i = 1; i <= N; i++) if (!T[i]) {
if (Lx[u] + Ly[i] == W[u][i]) {
T[i] = true;
if (!L[i] || match(L[i])) {
L[i] = u;
return true;
}
} else
slack[i] = min(slack[i], Lx[u]+Ly[i]-W[u][i]);
}
return false;
}
void update () {
int a = inf;
for (int i = 1; i <= N; i++) if (!T[i])
a = min(a, slack[i]);
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++) {
L[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++) {
for (int j = 1; j <= N; j++) slack[j] = inf;
while (true) {
for (int j = 1; j <= N; j++) S[j] = T[j] = false;
if (match(i)) break;
else update();
}
}
}
void init () {
for (int i = 1; i <= N; i++)
for (int j = 1; j <= N; j++)
W[i][j] = -inf;
int u, v;
for (int i = 1; i <= N; i++) {
while (scanf("%d", &u) == 1 && u) {
scanf("%d", &v);
W[i][u] = max(W[i][u], -v);
}
}
}
int main () {
while (scanf("%d", &N) == 1 && N) {
init();
KM();
bool flag = false;
for (int i = 1; i <= N; i++) {
if (L[i] && W[L[i]][i] == -inf)
flag = true;
}
if (flag) printf("N\n");
else {
int ans = 0;
for (int i = 1; i <= N; i++)
ans += (Lx[i] + Ly[i]);
printf("%d\n", -ans);
}
}
return 0;
}
权值变为负的即便为求最小值。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn = 105;
const int inf = 0x3f3f3f3f;
int N, L[maxn], Lx[maxn], Ly[maxn], W[maxn][maxn], slack[maxn];
bool S[maxn], T[maxn];
int match (int u) {
S[u] = true;
for (int i = 1; i <= N; i++) if (!T[i]) {
if (Lx[u] + Ly[i] == W[u][i]) {
T[i] = true;
if (!L[i] || match(L[i])) {
L[i] = u;
return true;
}
} else
slack[i] = min(slack[i], Lx[u]+Ly[i]-W[u][i]);
}
return false;
}
void update () {
int a = inf;
for (int i = 1; i <= N; i++) if (!T[i])
a = min(a, slack[i]);
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++) {
L[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++) {
for (int j = 1; j <= N; j++) slack[j] = inf;
while (true) {
for (int j = 1; j <= N; j++) S[j] = T[j] = false;
if (match(i)) break;
else update();
}
}
}
void init () {
for (int i = 1; i <= N; i++)
for (int j = 1; j <= N; j++)
W[i][j] = -inf;
int u, v;
for (int i = 1; i <= N; i++) {
while (scanf("%d", &u) == 1 && u) {
scanf("%d", &v);
W[i][u] = max(W[i][u], -v);
}
}
}
int main () {
while (scanf("%d", &N) == 1 && N) {
init();
KM();
bool flag = false;
for (int i = 1; i <= N; i++) {
if (L[i] && W[L[i]][i] == -inf)
flag = true;
}
if (flag) printf("N\n");
else {
int ans = 0;
for (int i = 1; i <= N; i++)
ans += (Lx[i] + Ly[i]);
printf("%d\n", -ans);
}
}
return 0;
}
相关文章推荐
- HDU 5458(Stability-LCT维护连通图2点间割边个数)
- 全面讲解Tomcat下SSL证书的配置(五)
- ubuntu下source、sh、bash、./执行脚本的区别
- linux(ftp)
- Linux DRM KMS 驱动简介
- Linux复习笔记(一) -- Bash的基本操作
- 云计算大赛准备过程之安装虚拟机及linux系统ubuntu过程
- nginx配置语法
- linux(知识)
- 分析OpenCV图像输入输出基本框架
- shell 解释器
- no ocijdbc11 in java.library.path linux
- Docker awesome commands
- linux中断原理专题
- Linux使用笔记: 定制core dump文件的文件名
- Android访问百度网站
- Install R and RStudio on Ubuntu
- linux面试题
- shell第一个例子清理日志
- Linux下面的终端编程 做一个简单的菜单