洛谷P2756 飞行员配对方案问题（网络流）

飞行员配对方案问题

题目背景

第二次世界大战时期..

题目描述

英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的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!’。

输入输出样例

输入样例#1：
5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1
输出样例#1：
4
1 7
2 9
3 8
5 10

分析：将英国飞行员连原点，外籍飞行员连汇点，容量都为1，可以配对的飞行员两两连边，容量为max，跑一遍网络流即可。

代码

#include <cstdio>
#include <queue>
#include <cstring>
#define inf 1e9
#define N 10005
using namespace std;

queue <int>q; 
struct arr
{
    int from,to,nxt,w;  
} a[N];
int ls[N],dis[N],cur[N],l,n,m,ans;

int min(int x,int y) {return x<y?x:y;}

void add(int x,int y,int z)
{
    a[++l].to=y;a[l].nxt=ls[x];a[l].w=z;ls[x]=cur[x]=l;a[l].from=x;
    a[++l].to=x;a[l].nxt=ls[y];a[l].w=0;ls[y]=cur[x]=l;a[l].from=y;
}

bool bfs()
{
    while (!q.empty()) q.pop();
    for (int i=0;i<=n+2;i++) dis[i]=0;
    dis[n+2]=1;
    q.push(n+2);
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=ls[x];i;i=a[i].nxt)
            if (a[i].w&&!dis[a[i].to])
            {
                dis[a[i].to]=dis[x]+1;
                if (a[i].to==n+1) return true;
                q.push(a[i].to);
            }
    }
    return false;
}

int dfs(int x,int maxf)
{
    if (x==n+1||!maxf) return maxf;
    int ret=0;
    for (int &i=cur[x];i;i=a[i].nxt)
        if (a[i].w&&dis[a[i].to]==dis[x]+1)
        {
            int f=dfs(a[i].to,min(a[i].w,maxf-ret));
            a[i].w-=f;
            a[i^1].w+=f;
            ret+=f;
            if (maxf==ret) break;
        }
    return ret;
}

void dinic()
{
    ans=0;
    while (bfs())
    {
        for (int i=1;i<=l;i++)
            cur[i]=ls[i];
        ans+=dfs(n+2,inf);
    }
}

int main()
{
    scanf("%d%d",&m,&n);
    int x,y;
    l=-1;
    while (scanf("%d%d",&x,&y))
    {
        if (x==-1&&y==-1) break;
        add(x,y,inf);
    }
    for (int i=1;i<=m;i++) add(n+2,i,1);
    for (int i=m+1;i<=n;i++) add(i,n+1,1);
    dinic();
    if (!ans) puts("No Solution!");
    else
    {
        printf("%d\n",ans);
        for (int i=0;i<=l;i+=2)
            if (a[i].to!=n+1&&a[i].to!=n+2&&a[i].from!=n+1&&a[i].from!=n+2&&a[i^1].w) 
                printf("%d %d\n",a[i].from,a[i].to);
    }
}