[网络流]NEW！飞行员配对方案问题

题目描述

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

如果所求的最佳飞行员配对方案不存在，则输出‘No Solution!’。

输入输出样例

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

分析

数据居然没有No solution的情况
裸网络流
飞行员之间的最大流量为1，然后建个超级源点连外籍，超级汇点连英国的即可

#include <iostream>
#include <cstdio>
#include <queue>
#define rep(i,a,b) for (i=a;i<=b;i++)
using namespace std;
int n,m;
int s,t;
int list[101],next[20001],u[20001],v[20001],w[20001];
int d[101];
int cur[101];
int ans;
void init()
{
    int x,y,o=3,i;
    scanf("%d%d",&m,&n);
    scanf("%d%d",&x,&y);
    while (!(x==y&&x==-1))
    {
        u[o]=x;v[o]=y;
        w[o]=1;
        next[o]=list[x];
        list[x]=o;
        u[o^1]=y;v[o^1]=x;
        next[o^1]=list[y];
        list[y]=o^1;
        o+=2;
        scanf("%d%d",&x,&y);
    }
    rep(i,1,m)
    {
        v[o]=i;
        w[o]=1;
        next[o]=list[0];
        list[0]=o;
        u[o^1]=i;
        next[o^1]=list[i];
        list[i]=o^1;
        o+=2;
    }
    rep(i,m+1,n)
    {
        u[o]=i;v[o]=n+1;
        w[o]=1;
        next[o]=list[i];
        list[i]=o;
        u[o^1]=n+1;v[o^1]=i;
        next[o^1]=list[n+1];
        list[n+1]=o^1;
        o+=2;
    }
    s=0;t=n+1;
}
bool bfs()
{
    int i,x;
    rep(i,s,t) d[i]=0;
    queue<int> q;
    while (!q.empty()) q.pop();
    d[0]=1;q.push(0);
    while (!q.empty())
    {
        x=q.front();q.pop();
        for (i=list[x];i;i=next[i])
        if (w[i]&&!d[v[i]])
        {
            d[v[i]]=d[x]+1;
            if (v[i]==t) return true;
            q.push(v[i]);
        }
    }
    return false;
}
int dfs(int p,int mf)
{
    int ret=0,flow;
    if (p==t||!mf) return mf;
    for (int &i=cur[p];i;i=next[i])
    if (d[p]+1==d[v[i]]&&w[i])
    {
        flow=dfs(v[i],min(w[i],mf-ret));
        w[i]-=flow;
        w[i^1]+=flow;
        ret+=flow;
        if (mf==ret) break;
    }
    return ret;
}
void dinic()
{
    int i;
    while (bfs())
    {
        rep(i,s,t) cur[i]=list[i];
        ans+=dfs(s,2147483647);
    }
}
int main()
{
    init();
    dinic();
    printf("%d",ans);
}