【网络流24题】飞行员配对方案问题

Description

第二次世界大战时期，英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的2 名飞行员，其中1 名是英国飞行员，另1 名是外籍飞行员。在众多的飞行员中，每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。如何选择配对飞行的飞行员才能使一次派出最多的飞机。对于给定的外籍飞行员与英国飞行员的配合情况，试设计一个算法找出最佳飞行员配对方案，使皇家空军一次能派出最多的飞机。

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

Input

由文件input.txt提供输入数据。文件第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!’。

Sample Input

5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1

Sample Output

4
1 7
2 9
3 8
5 10

题解

二分图最大匹配裸题，网络流最大流的模板题
二分图的话就是建边之后跑匈牙利算法
网络流就是在二分图的边的基础上再加上原点和汇点，并将所有的外籍飞行员与原点建边，英国飞行员与汇点建边，再随便用最大流的模板搞一搞就可以了

code

最大流dinic

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 10010
#define inf 0x7fffffff
using namespace std;
struct ss
{
    int next,to,s;
};
ss Edge[N<<1];
const int T=10000;
int n,m,tot=1,ans;
int head[N];
bool used[N];
int h[N];
int to[N];

void add(int x,int y,int z)
{
    Edge[++tot].next=head[x];
    Edge[tot].to=y;
    Edge[tot].s=z;
    head[x]=tot;
}

void ins(int x,int y)
{
    add(x,y,1),add(y,x,0);
}

bool bfs()
{
    memset(h,-1,sizeof(h));
    queue<int>Q;
    Q.push(0);
    h[0]=0;
    while(!Q.empty())
    {
        int now=Q.front();
        Q.pop();
        for(int i=head[now];i;i=Edge[i].next)
        {
            int to=Edge[i].to,s=Edge[i].s;
            if(s&&h[to]==-1)
            {
                h[to]=h[now]+1;
                Q.push(to);
            }
        }
    }
    if(h[T]==-1)return false;
    return true;
}

int dfs(int u,int f)
{
    if(u==T)return f;
    int ss=0;
    for(int i=head[u];i;i=Edge[i].next)
    {
        int to=Edge[i].to,s=Edge[i].s;
        if(!s||h[to]!=h[u]+1)continue;
        int w=f-ss;
        w=dfs(to,min(w,s));
        Edge[i].s-=w;
        Edge[i^1].s+=w;
        ss+=w;
        if(ss==f)return f;
    }
    if(!ss)h[u]=-1;
    return ss;
}

void dinic()
{
    while(bfs())
        ans+=dfs(0,inf);
}

void getans()
{
    for(int i=1;i<=n;i++)
    {
        for(int j=head[i];j;j=Edge[j].next)
        {
            if(Edge[j].s||!Edge[j].to)continue;
            printf("%d %d\n",i,Edge[j].to);
            break;
        }
    }
}

int main()
{
    cin>>n>>m;
    int x,y;
    while(scanf("%d%d",&x,&y)&&(x!=-1||y!=-1))
    {
        ins(x,y);
        if(!used[x])ins(0,x),used[x]=true;
        if(!used[y])ins(y,T),used[y]=true;
    }
    dinic();
    cout<<ans<<endl;
    getans();
    return 0;
}

匈牙利算法

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 10010
using namespace std;
struct node
{
    int next,to;
};
node Edge[N<<1];
int n,m,tot;
int head[N];
int f[N];
bool used[N];

void add(int x,int y)
{
    Edge[++tot].next=head[x];
    Edge[tot].to=y;
    head[x]=tot;
}

bool dfs(int u)
{
    for(int i=head[u];i;i=Edge[i].next)
    {
        int to=Edge[i].to;
        if(!used[to])
        {
            used[to]=true;
            if(f[to]==-1||dfs(f[to]))
            {
                f[to]=u;
                return true;
            }
        }
    }
    return false;
}

int hungary()
{
    int cnt=0;
    memset(f,-1,sizeof(f));
    for(int i=1;i<=n;i++)
    {
        memset(used,0,sizeof(used));
        if(dfs(i))cnt++;
    }
    return cnt;
}

int main()
{
    cin>>n>>m;
    int x,y;
    while(scanf("%d%d",&x,&y)&&(x!=-1||y!=-1))
        add(x,y);
    printf("%d\n",hungary());
    for(int i=n+1;i<=m;i++)
        if(f[i]!=-1)printf("%d %d\n",f[i],i);
    return 0;
}