飞行员配对方案问题
Description
第二次世界大战时期,英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的 2 名飞行员,其中 1 名是英国飞行员,另 1 名是外籍飞行员。在众多的飞行员中,每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。如何选择配对飞行的飞行员才能使一次派出最多的飞机。对于给定的外籍飞行员与英国飞行员的配合情况,试设计一个算法找出最佳飞行员配对方案,使皇家空军一次能派出最多的飞机。
对于给定的外籍飞行员与英国飞行员的配合情况,编程找出一个最佳飞行员配对方案,使皇家空军一次能派出最多的飞机。
Input
文件第 1 行有 2 个正整数 m 和 n。n 是皇家空军的飞行员总数(n<100);m 是外籍飞行员数。外籍飞行员编号为1~m;英国飞行员编号为 m+1~n。接下来每行有 2 个正整数 i 和 j,表示外籍飞行员 i 可以和英国飞行员 j 配合。文件最后以2个-1结束。
Output
程序运行结束时,将最佳飞行员配对方案输出到文件 output.txt 中。第 1 行是最佳飞行员配对方案一次能派出的最多的飞机数 M。接下来 M 行是最佳飞行员配对方案。每行有 2个正整数 i 和 j,表示在最佳飞行员配对方案中,飞行员 i 和飞行员 j 配对。如果所求的最佳飞行员配对方案不存在,则输出‘No Solution!’。
Sample Input
5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1
Sample Output
4
1 7
2 9
3 8
5 10
题目大意
求二分图最大匹配。
题解
建立一个超级源一个超级汇,源向每个外籍飞行员连一条容量为1的边,每个英国飞行员向汇连一条容量为1的边,每个可以组队的外籍飞行员和英国飞行员间连一条容量为1的边。然后跑最大流,满流边就是一种匹配。
但是由于不会写spj,只能用手上的数据删去输出方案测试通过。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 200 + 10, inf = 0x3f3f3f3f;
int n, m, tot, ans;
struct Edge{
int fr, to, cap, flow;
}edg[N*N];
int hd[N], nxt[N*N];
int s, t;
int d[N], q[N], vis[N], dfn, cur[N];
void insert(int u, int v, int w){
edg[tot].fr = u, edg[tot].to = v, edg[tot].cap = w;
nxt[tot] = hd[u]; hd[u] = tot;
tot++;
edg[tot].fr = v, edg[tot].to = u;
nxt[tot] = hd[v]; hd[v] = tot;
tot++;
}
void init(){
memset(hd, -1, sizeof(hd));
scanf("%d%d", &m, &n);
s = 0, t = n + 1;
int x, y;
while(scanf("%d%d", &x, &y) && x != -1) insert(x, y, 1);
for(int i = 1; i <= m; i++) insert(0, i, 1);
for(int i = m + 1; i <= n; i++) insert(i, n+1, 1);
}
bool bfs(){
int head = 1, tail = 1;
d[s] = 0; q[1] = s; vis[s] = ++dfn;
while(head <= tail){
int u = q[head++];
for(int i = hd[u]; i >= 0; i = nxt[i]){
Edge &e = edg[i];
if(vis[e.to] == dfn || e.cap <= e.flow) continue;
vis[e.to] = dfn;
d[e.to] = d[u] + 1;
q[++tail] = e.to;
}
}
return vis[t] == dfn;
}
int dfs(int x, int a){
if(x == t || a == 0) return a;
int flow = 0, f;
for(int &i = cur[x]; i >= 0; i = nxt[i]){
Edge &e = edg[i];
if(d[e.to] == d[x] + 1 && (f = dfs(e.to, min(a, e.cap-e.flow))) > 0){
e.flow += f;
edg[i^1].flow -= f;
flow += f;
a -= f;
if(a == 0) break;
}
}
return flow;
}
void work(){
while(bfs()){
for(int i = 0; i <= n+1; i++) cur[i] = hd[i];
ans += dfs(s, inf);
}
printf("%d\n", ans);
}
int main(){
freopen("prog81.in", "r", stdin);
freopen("prog81.out", "w", stdout);
init();
work();
return 0;
}