题目背景
第二次世界大战时期..
题目描述
英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的2 名飞行员,其中1 名是英国飞行员,另1名是外籍飞行员。在众多的飞行员中,每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。如何选择配对飞行的飞行员才能使一次派出最多的飞机。对于给定的外籍飞行员与英国飞行员的配合情况,试设计一个算法找出最佳飞行员配对方案,使皇家空军一次能派出最多的飞机。
对于给定的外籍飞行员与英国飞行员的配合情况,编程找出一个最佳飞行员配对方案,使皇家空军一次能派出最多的飞机。
输入输出格式
输入格式:第 1 行有 2 个正整数 m 和 n。n 是皇家空军的飞行员总数(n<100);m 是外籍飞行员数(m<=n)。外籍飞行员编号为 1~m;英国飞行员编号为 m+1~n。
接下来每行有 2 个正整数 i 和 j,表示外籍飞行员 i 可以和英国飞行员 j 配合。最后以 2个-1 结束。
第 1 行是最佳飞行员配对方案一次能派出的最多的飞机数 M。接下来 M 行是最佳飞行员配对方案。每行有 2个正整数 i 和 j,表示在最佳飞行员配对方案中,飞行员 i 和飞行员 j 配对。如果所求的最佳飞行员
对方案不存在,则输出‘No Solution!’
其实这不就是一个网络流的入门题吗233~
二分图匹配一下就好了呀~
唯一有一点问题的就是打出方案。。。
答: 中间两列飞行员之间的边的退流边大于0就代表这条路被操作过了。。。
输出它。。。OK
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int maxn = 208, INF = 0x7fffffff, s = 0, t = 206;
int ans = 0, cnt = -1, m, n, a, b;
int layer[maxn];
struct lpl
{
int to, dis;
}lin;
vector<int> point[maxn];
vector<lpl> edge;
queue<int> q;
inline void connect(int aa, int bb, int w)
{
cnt++; lin.to = bb; lin.dis = w; point[aa].push_back(cnt); edge.push_back(lin);
cnt++; lin.to = aa; lin.dis = 0; point[bb].push_back(cnt); edge.push_back(lin);
}
inline void putit()
{
scanf("%d%d", &m, &n);
for(;;)
{
scanf("%d%d", &a, &b);
if(a == -1) break;
connect(a, b, INF);
}
for(int i = 1; i <= m; ++i)
{
connect(s, i, 1);
}
for(int i = 1 + m; i <= n; ++i)
{
connect(i, t, 1);
}
}
inline bool bfs()
{
memset(layer, 0, sizeof(layer));
q.push(s); layer[s] = 1;
while(!q.empty())
{
int now = q.front(); q.pop();
for(int i = point[now].size() - 1; i >= 0; --i)
{
if(layer[edge[point[now][i]].to] == 0 && edge[point[now][i]].dis != 0)
{
layer[edge[point[now][i]].to] = layer[now] + 1;
q.push(edge[point[now][i]].to);
}
}
}
return layer[t];
}
inline int dfs(int aa, int w, int bb)
{
if(aa == bb) return w;
int ret = 0;
for(int i = point[aa].size() - 1; i >= 0; --i)
{
if(layer[edge[point[aa][i]].to] == layer[aa] + 1 && edge[point[aa][i]].dis != 0)
{
int tmp = dfs(edge[point[aa][i]].to, min(w - ret, edge[point[aa][i]].dis), t);
if(tmp == 0) continue;
ret += tmp; edge[point[aa][i]].dis -= tmp; edge[point[aa][i] ^ 1].dis += tmp;
}
}
return ret;
}
inline int Dinic()
{
int ret = 0;
while(bfs())
{
ret += dfs(s, INF, t);
}
return ret;
}
int main()
{
putit();
printf("%d\n", Dinic());
for(int i = 1; i <= m; ++i)
{
for(int j = point[i].size() - 1; j >= 0; --j)
{
if(edge[point[i][j]].to != s && edge[point[i][j] ^ 1].dis != 0)
{
printf("%d %d\n", i, edge[point[i][j]].to);
}
}
}
return 0;
}
5 10 1 7 1 8 2 6 2 9 2 10 3 7 3 8 4 7 4 8 5 10 -1 -1
4 1 7 2 9 3 8 5 10