飞行员配对方案问题

Problem
第二次世界大战时期，英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出
的每一架飞机都需要配备在航行技能和语言上能互相配合的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!’。
输入文件示例 输出文件示例
input.txt output.txt
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

题解：一眼题~~
如果两个飞行员之间有可以搭配，就连一条容量为1的边，再建原点汇点，连边之后跑一遍最大流，即是数量，关于输出方案貌似题目数据不全（存在多种匹配方案满足最优解），下面放代码>_<：

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
#define maxn 100 + 10
#define maxm 10000 + 10
#define INF 0x7f7f7f7f
struct edge{
    int to, next, cap;
}e[maxm];
int n, m, s, t, cur = 1, ans = 0, h[maxn], que[maxn], front[maxn], pa[maxn];
inline void addedge(int u,int v,int c){
     cur++;
     e[cur].to = v;
     e[cur].cap = c;
     e[cur].next = front[u];
     front[u] = cur;
     cur ++;
     e[cur].to = u;
     e[cur].next = front[v];
     front[v] = cur;
}
inline bool bfs(){
    int x, tp = 0, d = 1, p;
    memset(h, -1, sizeof(h));
    que[tp] = h[0] = 0;
    while(tp < d){
        x = que[tp]; tp ++;
        int tmp = front[x];
        while(tmp){
            if(e[tmp].cap && h[e[tmp].to] < 0){
                que[d ++] = e[tmp].to;
                h[e[tmp].to] = h[x] + 1;
            }
            tmp = e[tmp].next;
        }
    }
    if(h[n + 1] == -1) return 0;
    return 1;
}
inline int dfs(int x, int min_adv){
    if(x == n + 1) return min_adv;
    int tmp = front[x];
    int flow = 0, f;
    while(tmp){
        if(e[tmp].cap && h[e[tmp].to] == h[x] + 1) {
            f = min_adv - flow;
            f = dfs(e[tmp].to, min(f, e[tmp].cap));
            if(f) pa[x] = e[tmp].to;
            e[tmp].cap -= f;
            e[tmp ^ 1].cap += f;
            flow += f;
            if(flow == min_adv) return min_adv;
        }
        tmp = e[tmp].next;
    }
    if(! flow) h[x] = -1;
    return flow;
}
inline void dinic(){
     while(bfs()) ans += dfs(0,INF);
}
int main(){
    scanf("%d%d",&m,&n);
    int u, v; while(scanf("%d%d",&u, &v) == 2 && (u > 0)) addedge(u, v, 1);
    for(int i = 1; i <= m; i ++) addedge(0, i, 1);
    for(int i = m + 1;i <= n;i ++) addedge(i, n + 1, 1);
    dinic();
    printf("%d\n",ans);
    for(int i = 1;i <= m;i ++) if(pa[i]) printf("%d %d\n", i, pa[i]);
    return 0;
}