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

问题描述

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

编程任务

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

---

数据输入

由文件input.txt提供输入数据。文件第1 行有2个正整数m和n。n是皇家空军的飞行员总数(n<100)；m是外籍飞行员数。外籍飞行员编号为1~m；英国飞行员编号为m+1~n。
接下来每行有2 个正整数i和j，表示外籍飞行员i可以和英国飞行员j配合。文件最后以2个-1 结束。

---

结果输出

程序运行结束时，将最佳飞行员配对方案输出到文件output.txt 中。第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

---

解题思路

网络流模板题，二分图最大匹配。
另外设一个源点S，向每个左边的点连容量为1的边
以及设一个汇点T，每个右边的点向它连一条容量为1的边
中间的按照题目数据连容量为1的边。

---

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;

const int N=2100;
const int M=2000000;
const int INF=0x7fffffff;

struct note
{
    int x,y,c,next;
};

note side[M];
int last[N],cur[N],dis[N];
int n,m,ans,l=1,s,t;


void add(int x,int y)
{
    side[++l].x=x; side[l].y=y; side[l].c=1; 
    side[l].next=last[x]; last[x]=l;
    side[++l].x=y; side[l].y=x; side[l].c=0;
    side[l].next=last[y]; last[y]=l;
}

void init()
{
    memset(last,0,sizeof(last));
    scanf("%d%d",&n,&m);
    s=0; t=m+n+1;
    for (int i=1;i<=n;++i) add(s,i);
    for (int i=n+1;i<=m+n;++i) add(i,t);
  int x=0,y=0;
  while (x!=-1 && y!=-1)
  {
    scanf("%d%d",&x,&y);
    if (x==-1 && y==-1) break;
    add(x,y);
  }
}

bool bfs()
{
    memset(dis,0,sizeof(dis));
    queue<int> Q;
    dis[s]=1;
    Q.push(s);
    while (!Q.empty())
    {
        int now=Q.front();
        Q.pop();
        for (int i=last[now];i!=0;i=side[i].next)
            if (side[i].c && !dis[side[i].y])
            {
                int j=side[i].y;
                dis[j]=dis[now]+1;
                //if (j==t) return true;
                Q.push(j);
            }
  }
  if (dis[t]==0) return false;
    else return true;
}

int dfs(int x,int maxf)
{
    if (x==t || maxf==0) return maxf;
    int ret=0;
    for (int i=cur[x];i!=0;i=side[i].next)
      if (side[i].c && dis[side[i].y]==dis[x]+1)
      {
        int f=dfs(side[i].y,min(maxf-ret,side[i].c));
        side[i].c-=f;
        side[i^1].c+=f;
        ret+=f;
        if (ret==maxf) break;
      }
    return ret;
}

void dinic()
{
    while (bfs())
    {
        for (int i=s;i<=t;++i)
            cur[i]=last[i];
        ans+=dfs(s,INF);
    }
}

int main()
{
    init();
    ans=0;
    dinic();
    if (ans==0) printf("No Solution!\n");
      else printf("%d\n",ans);
    return 0;
}

---

PS

边集数组一定要定好初值，让i和i^1为一对反向边，不然就会错。
                        ↑↑↑↑↑↑↑↑↑↑ 
因为这个调了n久，发现后非常不爽。。。。