UVa 二分图匹配 Biginners
UVa 1045 - The Great Wall Game 最小权匹配
题意:给你一个n*n的棋盘,上面有n个棋子,要求通过移动各个棋子使得棋子在同一行或者同一列或者对角线上,求最小移动次数。
思路:直接对于所有可能情况构造二分图,X集合为最初棋子,Y集合为移动后的棋子方位,边权为移动的次数。然后KM算法求最小权匹配。
/* **********************************************
Author : JayYe
Created Time: 2013-8-18 15:55:41
File Name : zzz.cpp
*********************************************** */ #include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std; int max(int a, int b) { return a>b?a:b; }
int min(int a, int b) { return a>b?b:a; } const int maxn = 22; struct PP {
int x, y;
}a[maxn]; int n, slack[maxn], lx[maxn], ly[maxn], match[maxn], w[maxn][maxn];
bool S[maxn], T[maxn]; bool dfs(int i) {
S[i] = 1;
for(int j = 1;j <= n; j++) if(!T[j])
slack[j] = min(slack[j], w[i][j] - lx[i] - ly[j]);
for(int j = 1;j <= n; j++) if(w[i][j] == lx[i] + ly[j] && !T[j]) {
T[j] = 1;
if(!match[j] || dfs(match[j])) {
match[j] = i;
return true;
}
}
return false;
} void update() {
int delta = 1<<22;
for(int i = 1;i <= n; i++) if(!T[i])
delta = min(delta, slack[i]);
for(int i = 1;i <= n; i++) {
if(S[i]) lx[i] += delta;
if(T[i]) ly[i] -= delta;
}
} void KM() {
int i, j;
for(i = 1;i <= n; i++) {
ly[i] = match[i] = 0;
lx[i] = 1<<22;
for(j = 1;j <= n; j++)
lx[i] = min(lx[i], w[i][j]);
}
for(i = 1;i <= n; i++) {
while(true) {
for(j = 1;j <= n; j++) S[j] = T[j] = 0, slack[j] = 1<<22;
if(dfs(i)) break;
else update();
}
}
} int solve() {
int i, j, k; for(i = 1;i <= n; i++)
scanf("%d%d" ,&a[i].x, &a[i].y);
int ans = 1<<22;
// 棋子在同一行的情况
for(i = 1;i <= n; i++) {
for(j = 1;j <= n; j++)
for(k = 1;k <= n; k++)
w[j][k] = abs(i - a[j].x) + abs(k - a[j].y);
KM();
int cur = 0;
for(j = 1;j <= n; j++) cur += lx[j] + ly[j];
ans = min(ans, cur);
}
// 棋子在同一列的情况
for(i = 1;i <= n; i++) {
for(j = 1;j <= n; j++)
for(k = 1;k <= n; k++)
w[j][k] = abs(k - a[j].x) + abs(i - a[j].y);
KM();
int cur = 0;
for(j = 1;j <= n; j++) cur += lx[j] + ly[j];
ans = min(ans, cur);
}
// 棋子在对角线的两种情况
for(i = 1;i <= n; i++)
for(j = 1;j <= n; j++)
w[i][j] = abs(j - a[i].x) + abs(j - a[i].y);
KM();
int cur = 0;
for(i = 1;i <= n; i++) cur += lx[i] + ly[i];
ans = min(ans, cur); for(i = 1;i <= n; i++)
for(j =1;j <= n; j++)
w[i][j] = abs(j - a[i].x) + abs(n-j+1 - a[i].y);
KM();
cur = 0;
for(i = 1;i <= n; i++) cur += lx[i] + ly[i];
ans = min(ans, cur);
return ans;
} int main() {
int cas = 1;
while(scanf("%d", &n) != -1 && n) {
printf("Board %d: %d moves required.\n\n", cas++, solve());
}
return 0;
}
UVa 12168 - Cat vs. Dog 最大独立集
根据题意直接构造二分图,X集合表示喜欢狗的人,Y集合表示喜欢猫的人,如果X集合里的人不喜欢某只猫,就与Y集合里喜欢该猫的人连边,反之也一样。那么要使得尽可能多的人满足愿望,也就是求最大独立集。
最大独立集 = N - 最大匹配
/* **********************************************
Author : JayYe
Created Time: 2013-8-18 16:53:58
File Name : zzz.cpp
*********************************************** */ #include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std; const int maxn = 500+10;
struct PP {
int x, y;
PP() {}
PP(int x, int y) : x(x), y(y) {}
}cat[maxn], dog[maxn]; int n, m, match[maxn];
bool vis[maxn], mp[maxn][maxn]; bool dfs(int i) {
for(int j = 1;j <= m; j++) if(mp[i][j] && !vis[j]) {
vis[j] = 1;
if(!match[j] || dfs(match[j])) {
match[j] = i;
return true;
}
}
return false;
} int main() {
int i, j, t, c, d;
scanf("%d", &t);
while(t--) {
scanf("%d%d%d", &c, &d, &n);
int n1 = 0, n2 = 0, x, y;
for(i = 1;i <= n; i++) {
char ch1[11], ch2[11];
scanf("%s%s", ch1, ch2);
if(ch1[0] == 'C') {
sscanf(ch1+1, "%d", &x);
sscanf(ch2+1, "%d", &y);
cat[++n1] = PP(x, y);
}
else {
sscanf(ch1+1, "%d", &x);
sscanf(ch2+1, "%d", &y);
dog[++n2] = PP(x, y);
}
}
n = n1, m = n2;
if(n == 0 || m == 0) {
printf("%d\n", n+m); continue;
}
for(i = 1;i <= n; i++) {
for(j = 1;j <= m; j++) {
if(cat[i].y == dog[j].x || cat[i].x == dog[j].y)
mp[i][j] = 1;
else
mp[i][j] = 0;
}
}
for(i = 1;i <= m; i++) match[i] = 0;
int ans = 0;
for(i = 1;i <= n; i++) {
for(j = 1;j <= m; j++) vis[j] = 0;
if(dfs(i)) ans++;
}
printf("%d\n", n+m - ans);
}
return 0;
}
Uva 1349 - Optimal Bus Route Design
题意:给你一个n个点的有向带环图,要你找出几个圈使得每个结点只属于一个圈,要求输出最小的总的长度,如果没有这样的方案,输出N。
构造二分图,把所有的结点拆成两个,放在X集合的为i, 放在Y集合的为 i ',如果有边i -> j,则在图中引入边i -> j',这样子构造好后实际上就是求最小权完美匹配,如果没有完美匹配则无解。
这里需要注意的是输入的边可能有好多条是重复的但是权值不同,需要取最小权,这个wa了我好几发。。。
/* **********************************************
Author : JayYe
Created Time: 2013-8-18 17:31:11
File Name : zzz.cpp
*********************************************** */ #include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std; int max(int a, int b) { return a>b?a:b; }
int min(int a, int b) { return a>b?b:a; } const int INF = 1<<22;
const int maxn = 100+10;
int n, slack[maxn], match[maxn], lx[maxn], ly[maxn], w[maxn][maxn];
bool S[maxn], T[maxn]; bool dfs(int i) {
S[i] = 1;
for(int j = 1;j <= n; j++) if(!T[j])
slack[j] = min(slack[j], w[i][j] - lx[i] - ly[j]);
for(int j = 1;j <= n; j++) if(w[i][j] == lx[i] + ly[j] && !T[j]) {
T[j] = 1;
if(!match[j] || dfs(match[j])) {
match[j] = i;
return true;
}
}
return false;
} void update() {
int delta = INF;
for(int i = 1;i <= n; i++) if(!T[i])
delta = min(delta, slack[i]);
for(int i = 1;i <= n; i++) {
if(S[i]) lx[i] += delta;
if(T[i]) ly[i] -= delta;
}
} void KM() {
int i, j;
for(i = 1;i <= n; i++) {
ly[i] = match[i] = 0;
lx[i] = INF;
for(j = 1;j <= n; j++)
lx[i] = min(lx[i], w[i][j]);
} for(i = 1;i <= n; i++) {
while(true) {
for(j = 1;j <= n; j++) S[j] = T[j] = 0, slack[j] = INF;
if(dfs(i)) break;
else update();
}
}
} void solve() {
int i, j;
for(i = 1;i <= n; i++)
for(j = 1;j <= n; j++)
w[i][j] = INF;
for(i = 1;i <= n; i++) {
while(scanf("%d", &j) && j) {
int dis;
scanf("%d", &dis);
// 同一条路要取最小值
w[i][j] = min(w[i][j], dis);
}
}
KM();
int ans = 0;
for(i = 1;i <= n; i++) ans += lx[i] + ly[i];
if(ans > INF-10) puts("N");
else printf("%d\n", ans);
} int main() {
while(scanf("%d", &n) != -1 && n) {
solve();
}
return 0;
}
UVa 二分图匹配 Biginners的更多相关文章
- UVa 二分图匹配 Examples
这些都是刘汝佳的算法训练指南上的例题,基本包括了常见的几种二分图匹配的算法. 二分图是这样一个图,顶点分成两个不相交的集合X , Y中,其中同一个集合中没有边,所有的边关联在两个集合中. 给定一个二分 ...
- POJ 2289 Jamie's Contact Groups / UVA 1345 Jamie's Contact Groups / ZOJ 2399 Jamie's Contact Groups / HDU 1699 Jamie's Contact Groups / SCU 1996 Jamie's Contact Groups (二分,二分图匹配)
POJ 2289 Jamie's Contact Groups / UVA 1345 Jamie's Contact Groups / ZOJ 2399 Jamie's Contact Groups ...
- uva 12083 Guardian of Decency (二分图匹配)
uva 12083 Guardian of Decency Description Frank N. Stein is a very conservative high-school teacher. ...
- UVA 12549 - 二分图匹配
题意:给定一个Y行X列的网格,网格种有重要位置和障碍物.要求用最少的机器人看守所有重要的位置,每个机器人放在一个格子里,面朝上下左右四个方向之一发出激光直到射到障碍物为止,沿途都是看守范围.机器人不会 ...
- UVA 1663 Purifying Machine (二分图匹配,最大流)
题意: 给m个长度为n的模板串,模板串由0和1和*三种组成,且每串至多1个*,代表可0可1.模板串至多匹配2个串,即*号改成0和1,如果没有*号则只能匹配自己.问:模板串可以缩减为几个,同样可以匹配原 ...
- UVA 11045-My T-shirt suits me(二分图匹配)
题意:有N件T恤,N是6的倍数,因为有6种型号,每种件数相同,有M个人,每个人有两种型号的T恤适合他,每个人可以挑其中的一种,问能否所有的人都能分配到T恤. 解析:典型的二分图匹配,每N/6为同种T恤 ...
- POJ 1274 裸二分图匹配
题意:每头奶牛都只愿意在她们喜欢的那些牛栏中产奶,告诉每头奶牛愿意产奶的牛棚编号,求出最多能分配到的牛栏的数量. 分析:直接二分图匹配: #include<stdio.h> #includ ...
- BZOJ1433 ZJOI2009 假期的宿舍 二分图匹配
1433: [ZJOI2009]假期的宿舍 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2375 Solved: 1005[Submit][Sta ...
- HDU1281-棋盘游戏-二分图匹配
先跑一个二分图匹配,然后一一删去匹配上的边,看能不能达到最大匹配数,不能这条边就是重要边 /*----------------------------------------------------- ...
随机推荐
- support STL Viewer with WordPress On SAE
由于SAE不支持本地代码目录写入, 我把WordPress的uploads路径改到了Storage中, 使用Domain来存放非代码资源. 这导致STL Viewer插件无法正常使用. 解决方法: 把 ...
- Square(hdu 1511)
题目描述: Problem Description Given a set of sticks of various lengths, is it possible to join them end- ...
- poj 1068 Parencodings(栈)
题目链接:http://poj.org/problem?id=1068 思路分析:对栈的模拟,将栈中元素视为广义表,如 (((()()()))),可以看做 LS =< a1, a2..., a1 ...
- POJ 3518 Prime Gap(素数题)
[题意简述]:输入一个数,假设这个数是素数就输出0,假设不是素数就输出离它近期的两个素数的差值,叫做Prime Gap. [分析]:这题过得非常险.由于我是打的素数表. 由于最大的素数是1299709 ...
- 详细解析BluetoothAdapter的详细api
(1)开关状态值 (2)扫描状态值 (3)蓝牙操作接收的广播 (4)蓝牙操作请求的广播 (5)附加域 (6)错误码 (1)获取蓝牙适配器 (2)获取state状态方法 (3)蓝牙是否可用 (4)打开蓝 ...
- 驱动: oops
linux驱动调试--段错误之oops信息分析 http://blog.chinaunix.net/xmlrpc.php?r=blog/article&uid=29401328&id= ...
- [Swust OJ 566]--开N方数(牛顿切线法解高次方程)
题目链接:http://acm.swust.edu.cn/problem/0566/ Time limit(ms): 1000 Memory limit(kb): 65535 Descriptio ...
- [POJ 2588]--Snakes(并查集)
题目链接:http://poj.org/problem?id=2588 Snakes Time Limit: 1000MS Memory Limit: 65536K Description B ...
- 触发器应用 trigger
首先有一张表: create table T_SALARY ( name VARCHAR2(20), age NUMBER(2), salary NUMBER(5) ); insert into t_ ...
- Hadoop Error:Name node is in safe mode的解决方法
Error:name node is in safe mode 解决方法:hadoop dfsadmin -safemode leave(见图)