P2756 飞行员配对方案问题(最大流)

题目背景

第二次世界大战期间,英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的两名飞行员,其中一名是英国飞行员,另一名是外籍飞行员。在众多的飞行员中,每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。

题目描述

一共有 nn 个飞行员,其中有 mm 个外籍飞行员和 (n - m)(nm) 个英国飞行员,外籍飞行员从 11 到 mm 编号,英国飞行员从 m + 1m+1 到 nn 编号。 对于给定的外籍飞行员与英国飞行员的配合情况,试设计一个算法找出最佳飞行员配对方案,使皇家空军一次能派出最多的飞机。

输入格式

输入的第一行是用空格隔开的两个正整数,分别代表外籍飞行员的个数 mm 和飞行员总数 nn。
从第二行起到倒数第二行,每行有两个整数 u, vu,v,代表外籍飞行员 uu 可以和英国飞行员 vv 配合。
输入的最后一行保证为 -1 -1,代表输入结束。

输出格式

本题存在 Special Judge。
请输出能派出最多的飞机数量,并给出一种可行的方案。
输出的第一行是一个整数,代表一次能派出的最多飞机数量,设这个整数是 kk。
第 22 行到第 k + 1k+1 行,每行输出两个整数 u, vu,v,代表在你给出的方案中,外籍飞行员 uu 和英国飞行员 vv 配合。这 kk 行的 uu 与 vv 应该互不相同。

 

题解

源点向每个外籍飞行员建边,容量为1。

每个英国飞行员向汇点建边,容量为1。

有关系的外籍飞行员和英国飞行员建边,容量为inf。

判断一条边是否有贡献的办法是枚举每条边,看反向边的容量,如果不为0则说明有贡献,这样就可以输出方案。

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+100;
const int inf=1e9;
int n,m;
struct node {
    int u,v,w,nxt;
}edge[maxn];
int tot;
int head[maxn];
void addedge (int u,int v,int w) {
    edge[tot].u=u;
    edge[tot].v=v;
    edge[tot].w=w;
    edge[tot].nxt=head[u];
    head[u]=tot++;
    
    edge[tot].u=v;
    edge[tot].v=u;
    edge[tot].w=0;
    edge[tot].nxt=head[v];
    head[v]=tot++;
} 

int dep[maxn];
int inq[maxn];
int cur[maxn];
int wjm;
int maxflow=0;
int s,t;
bool bfs () {
    for (int i=0;i<=n+1;i++) {
        cur[i]=head[i];
        dep[i]=inf;
        inq[i]=0;
    }
    dep[s]=0;
    queue<int> q;
    q.push(s);
    while (!q.empty()) {
        int u=q.front();
        q.pop();
        inq[u]=0;
        for (int i=head[u];i!=-1;i=edge[i].nxt) {
            int v=edge[i].v;
            if (dep[v]>dep[u]+1&&edge[i].w) {
                dep[v]=dep[u]+1;
                if (inq[v]==0) {
                    q.push(v);
                    inq[v]=1;
                } 
            }
        }
    }
    if (dep[t]!=inf) return 1;
    return 0;
}
int dfs (int u,int flow) {
    int increase=0;
    if (u==t) {
        wjm=1;
        maxflow+=flow;
        return flow;
    }
    int used=0;
    for (int i=cur[u];i!=-1;i=edge[i].nxt) {
        cur[u]=i;
        int v=edge[i].v;
        if (edge[i].w&&dep[v]==dep[u]+1) {
            if (increase=dfs(v,min(flow-used,edge[i].w))) {
                used+=increase;
                edge[i].w-=increase;
                edge[i^1].w+=increase;
                if (used==flow) break;
            }
        }
    }
    return used;
}
int Dinic () {
    maxflow=0;
    while (bfs()) {
        wjm=1;
        while (wjm==1) {
            wjm=0;
            dfs(s,inf);
        }
    }
    return maxflow;
}


int main () {
    scanf("%d%d",&m,&n);
    for (int i=0;i<=n+1;i++) head[i]=-1;
    s=0,t=n+1;
    while (1) {
        int x,y;
        scanf("%d%d",&x,&y);
        if (x==-1&&y==-1) break;
        addedge(x,y,inf);
    }
    for (int i=1;i<=m;i++) addedge(s,i,1);
    for (int i=m+1;i<=n;i++) addedge(i,t,1);
    int ans=Dinic();
    
    printf("%d
",ans);
    for (int i=1;i<=tot;i+=2) {
        if (edge[i].v==s||edge[i].v==t) continue;
        if (edge[i].u==s||edge[i].u==t) continue;
        if (edge[i].w!=0) {
            printf("%d %d
",edge[i].v,edge[i].u);
        }
    } 
    
}
原文地址:https://www.cnblogs.com/zhanglichen/p/13324999.html