acwing-2175. 飞行员配对方案问题(二分图|最大流)

第二次世界大战时期，英国皇家空军从沦陷国征募了大量外籍飞行员。

由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的 2 名飞行员，其中 1 名是英国飞行员，另 1 名是外籍飞行员。

在众多的飞行员中，每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。

如何选择配对飞行的飞行员才能使一次派出最多的飞机。

对于给定的外籍飞行员与英国飞行员的配合情况，试设计一个算法找出最佳飞行员配对方案，使皇家空军一次能派出最多的飞机。

输入格式
第 1 行有 2 个正整数 m 和 n。m 是外籍飞行员数；n 是皇家空军的飞行员总数。

外籍飞行员编号为 1∼m；英国飞行员编号为 m+1∼n。

接下来每行有 2 个正整数 i 和 j，表示外籍飞行员 i 可以和英国飞行员 j 配合。

文件最后以 2 个 −1 结束。

输出格式
第 1 行是最佳飞行员配对方案一次能派出的最多的飞机数 M。

接下来 M 行是最佳飞行员配对方案。

每行有 2 个正整数 i 和 j，表示在最佳飞行员配对方案中，外籍飞行员 i 和英国飞行员 j 配对。

如果有多种配对方案，则输出任意一种即可，方案内部配对输出顺序随意。

数据范围
1<m<n<100

输入样例：
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

题解
二分图，匈牙利算法，时间复杂度O(nm)

#include<bits/stdc++.h>
using namespace std;
int n,m;
const int N = 1e2 + 10;
int match[N],vis[N];
int g[N][N];
bool find(int x){
    for(int i = 1;i <= n;i ++){
        if(vis[i] || !g[x][i])continue;
        vis[i] = true;
        if(!match[i] || find(match[i])){
            match[i] = x;
            return true;
        }
    }
    return false;
}
int main(){
    cin>>m>>n;
    int x,y;
    while(cin>>x>>y,x != -1 && y != -1){
        g[x][y - m] = 1;
    }
    int res = 0;
    for(int i = 1;i <= n;i ++){
        memset(vis,0,sizeof vis);
        if(find(i))res ++;
    }
    cout<<res<<endl;
    for(int i = 1;i <= n - m;i ++){
        if(match[i])cout<<match[i]<<" "<<i + m<<endl;
    }
    return 0;
}

2 题解
最大流

#include<bits/stdc++.h>
using namespace std;
const int N = 210;
const int M = 1e5 + 10;
const int INF = 0x3f3f3f3f;
struct Edge{
    int v,next,w;
}edge[M];
int head[N],cnt;
int n,m,s,e;
void add(int u,int v,int w){
    edge[cnt].v = v;
    edge[cnt].w = w;
    edge[cnt].next = head[u];
    head[u] = cnt ++;
}
int q[N],hh = 0,tt = 0;
int d[N],cur[N];
bool bfs(){
    tt = hh = 0;
    memset(d,-1,sizeof d);
    q[tt ++] = s,d[s] = 0,cur[s] = head[s];
    while(hh < tt){
        int t = q[hh ++];
        for(int i = head[t];~i;i = edge[i].next){
            int v = edge[i].v,w = edge[i].w;
            if(d[v] == -1 && w > 0){
                d[v] = d[t] + 1;
                cur[v] = head[v];
                if(v == e)return true;
                q[tt ++] = v;
            }
        }
    }
    return false;
}
int dfs(int u,int limit){
    if(u == e)return limit;
    int flow = 0;
    for(int i = cur[u];~i && flow < limit;i = edge[i].next){
        int v = edge[i].v,w = edge[i].w;
        cur[u] = i;
        if(d[v] == d[u] + 1 && w > 0){
            int t = dfs(v,min(limit - flow,w));
            if(!t)d[v] = -1;
            flow += t,edge[i].w -= t,edge[i ^ 1].w += t;
        }
    }
    return flow;
}
int dinic(){
    int maxflow = 0,flow = 0;
    while(bfs())
        while(flow = dfs(s,INF))maxflow += flow;
    return maxflow;
}
int main(){
    memset(head,-1,sizeof head);
    cin>>m>>n;
    s = 0,e = n + 1;
    for(int i = 1;i <= m;i ++)add(0,i,1),add(i,0,0);
    for(int i = m + 1;i <= n;i ++)add(i,n + 1,1),add(n + 1,i,0);
    int x,y;
    while(cin>>x>>y,x != -1 && y != -1){
        add(x,y,1);
        add(y,x,0);
    }
    cout<<dinic()<<endl;
    for(int i = 0;i < cnt;i += 2){
        if(edge[i].w == 0 && edge[i].v > m && edge[i].v < n + 1){
            cout<<edge[i ^ 1].v<<" "<<edge[i].v<<endl;
        }
    }
    return 0;
}