网络流24题 P2756 飞行员配对方案问题 最大流或最大匹配 两种方法过

法1：二分图模板题

ac代码：30ms

```

#include <stdio.h>
#include <string.h>

using namespace std;
const int maxn=210;
int h[maxn],e[maxn*4],ne[maxn*4],idx;
int match[maxn];
bool vis[maxn];
int n,m,a,b;

void add(int a,int b){
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
int find(int x){
    for(int i=h[x];~i;i=ne[i]){
        int j=e[i];
        if(!vis[j]){
            vis[j]=1;
            if(!match[j]||find(match[j])){
                match[j]=x;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(h,-1,sizeof h);
    while(scanf("%d%d",&a,&b)&&(~a||~b)){
        add(a,b);
    }
    int ans=0;
    for(int i=1;i<=n;i++){
        memset(vis,0,sizeof vis);
        if(find(i))ans++;
    }
    if(!ans)printf("No Solution!\n");
    else {
        printf("%d\n",ans);
        for(int i=n+1;i<=m;i++){
            if(match[i])printf("%d %d\n",match[i],i);
        }
    }
    return 0;
}

```

 

法2： 最大流

最大流真的神奇，感觉很多最优解，极值都可以用网络流来解，关键是建图，而建图恰好是网络流的难点，路漫漫其修远兮吾将上下而求索，往后要多加思考多加积累才行啊，而这题是经典的二分图转网络流， 引入超级源点s超级汇点 t，s到A组飞行员路径为1，A组飞行员到B组飞行员为inf，B组飞行员到t为1；反向边初始化为0；就可以求最大流就完事了，最后路径遍历一遍所有边，输出A到B反向边为1的就完事了；

这种做法看似比最大匹配代码要多几行，但是匹配的飞行员多几组，最大匹配就很难去写代码了，而网络流依然轻松解决；

以下给出ac代码：30ms

 ```

#include <stdio.h>
#include <string.h>
#include <queue>
using namespace std;
const int maxn=210,maxm=1e5+7,inf=0x3f3f3f3f;
int h[maxn],e[maxn*10],ne[maxn*10],w[maxn*10],idx;
int dep[maxn];
int n,m,s,t;

void add(int a,int b,int c){
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}
bool bfs(int s,int t){
    memset(dep,-1,sizeof dep);
    queue<int>q;q.push(s);dep[s]=1;
    while(q.size()){
        int u=q.front();q.pop();
        for(int i=h[u];~i;i=ne[i]){
            int j=e[i];
            if(dep[j]==-1&&w[i])dep[j]=dep[u]+1,q.push(j);
        }
    }
    return dep[t]!=-1;
}
int dfs(int x,int v){
    if(x==t||v==0)return v;
    int f,res=0;
    for(int i=h[x];~i;i=ne[i]){
        int j=e[i];
        if(dep[j]==dep[x]+1){
            f=dfs(j,min(v,w[i]));
            w[i]-=f,w[i^1]+=f,v-=f,res+=f;
            if(v==0)break;
        }
    }
    return res;
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(h,-1,sizeof h);
    s=0,t=m+1;
    int a,b;
    while(~scanf("%d%d",&a,&b)&&(~a||~b)){
        add(a,b,inf);add(b,a,0);
    }
    for(int i=1;i<=n;i++){
        add(s,i,1);add(i,s,0);
    }
    for(int i=n+1;i<=m;i++){
        add(i,t,1);add(t,i,0);
    }
    int ans=0;
    while(bfs(s,t))ans+=dfs(s,inf);
    printf("%d\n",ans);
    for(int i=0;i<idx;i+=2){
        if(e[i]!=t&&e[i]!=s&&e[i^1]!=s&&e[i^1]!=t&&w[i^1])printf("%d %d\n",e[i^1],e[i]);
    }
    return 0;
}
 

 

```