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

1. [网络流24题] 搭配飞行员
   ★★☆ 输入文件：flyer.in 输出文件：flyer.out 简单对比
   时间限制：1 s 内存限制：128 MB
   【问题描述】
   飞行大队有若干个来自各地的驾驶员，专门驾驶一种型号的飞机，这种飞机每架有两个驾驶员,需一个正驾驶员和一个副驾驶员。由于种种原因，例如相互配合的问题，有些驾驶员不能在同一架飞机上飞行，问如何搭配驾驶员才能使出航的飞机最多。
   
   如图，假设有10个驾驶员，如图中的V1，V2，…，V10就代表达10个驾驶员,其中V1，V2，V3，V4，V5是正驾驶员,V6，V7，V8，V9，V10是副驾驶员。如果一个正驾驶员和一个副驾驶员可以同机飞行，就在代表他们两个之间连一条线,两个人不能同机飞行，就不连。例如V1和V7可以同机飞行，而V1和V8就不行。请搭配飞行员，使出航的飞机最多。注意:因为驾驶工作分工严格,两个正驾驶员或两个副驾驶员都不能同机飞行.
   【输入格式】
   输入文件有若干行
   第一行，两个整数n与n1，表示共有n个飞行员(2<=n<=100),其中有n1名飞行员是正驾驶员.
   下面有若干行,每行有2个数字a,b。表示正驾驶员a和副驾驶员b可以同机飞行。
   注:正驾驶员的编号在前,即正驾驶员的编号小于副驾驶员的编号.
   【输出格式】
   输出文件有一行
   第一行，1个整数，表示最大起飞的飞机数。
   【输入输出样例】
   输入文件名： flyer.in
   10 5
   1 7
   2 6
   2 10
   3 7
   4 8
   5 9
   输出文件名：flyer.out
   4
   配对的飞行员之间连边，如果要用最大流，把源点0与正驾驶员连边，把汇点n+1与副驾驶员连边，流量都为1，如果要用二分图匹配，直接连边完跑一遍匈牙利算法就行了。。
   注意点：while ((x=read(),y=read())!=EOF)不知道为什么会死循环。。
   网络流dinic：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
using namespace std;
int m,n,i,x,y,ans,sum,cnt=1,h[100001]={0},head,tail,q[100001],dis[100001];
struct node
{
    int to,next,v;
};
node edge[100001];

int read()
{
    int w=0,c=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-')
          c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      {
        w=w*10+ch-'0';
        ch=getchar();
      }
    return w*c;
}

void add(int u,int v,int w)
{
    cnt++;
    edge[cnt].next=h[u];
    h[u]=cnt;
    edge[cnt].to=v;
    edge[cnt].v=w;
}

bool bfs()
{
    int j,p;
    memset(dis,-1,sizeof(dis));
    q[0]=0;
    dis[0]=0;
    head=0;
    tail=1;
    while (head<tail)
      {
        head++;
        j=q[head];
        p=h[j];
        while (p)
          {
            if (dis[edge[p].to]<0 && edge[p].v>0)
              {
                dis[edge[p].to]=dis[j]+1;
                tail++;
                q[tail]=edge[p].to;
              }
            p=edge[p].next;
          }
      }
    if (dis[n+1]>0)
      return true;
    else
      return false;
}

int dfs(int x,int f)
{
    int w,used=0,i=h[x];
    if (x==n+1)
      return f;
    while (i)
      {
        if (edge[i].v && dis[edge[i].to]==dis[x]+1)
          {
            w=f-used;
            w=dfs(edge[i].to,min(w,edge[i].v));
            edge[i].v-=w;
            edge[i^1].v+=w;
            used+=w;
            if (used==f)
              return f;
          }
        i=edge[i].next;
      }
    if (!used)
      dis[x]=-1;
    return used;
}

int main()
{
    freopen("flyer.in","r",stdin);
    freopen("flyer.out","w",stdout);
    n=read();
    m=read();
    for (i=1;i<=m;i++)
      {
        add(0,i,1);
        add(i,0,1);
      }
    for (i=m+1;i<=n;i++)
      {
        add(i,n+1,1);
        add(n+1,i,1);
      }
    while ((scanf("%d %d",&x,&y))!=EOF)
      {
        add(x,y,1);
        add(y,x,1);
      }
    ans=0;
    while (bfs())
      while (sum=dfs(0,0x7fffffff))
        ans+=sum;
    printf("%d",ans);  
    fclose(stdin);
    fclose(stdout);
    return 0;
}

二分图匹配：

#include<cstdio>
#include<cstring>
using namespace std;
int cnt=1,head[100001]={0},part[100001]={0},n,m,i,x,y,ans;
struct node
{
    int to,next;
};
node edge[100001]={0};
bool used[100001];

int read()
{
    int w=0,c=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-')
          c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      {
        w=w*10+ch-'0';
        ch=getchar();
      }
    return w*c;
}

void add(int u,int v)
{
    cnt++;
    edge[cnt].next=head[u];
    head[u]=cnt;
    edge[cnt].to=v;
}

bool find(int s)
{
    int k;
    for (k=head[s];k;k=edge[k].next)
      if (!used[edge[k].to])
        {
            used[edge[k].to]=true;
            if (!part[edge[k].to] || find(part[edge[k].to]))
              {
                part[edge[k].to]=s;
                return true;
              }
        }
    return false;
}

int main()
{
    freopen("flyer.in","r",stdin);
    freopen("flyer.out","w",stdout);
    n=read();
    m=read();
    while ((scanf("%d %d",&x,&y))!=EOF)
      {
        add(x,y);
        add(y,x);
      }
    ans=0;
    for (i=1;i<=m;i++)
      {
        memset(used,false,sizeof(used));
        if (find(i))
          ans++;
      }
    printf("%d",ans);
    fclose(stdin);
    fclose(stdout);
    return 0;
}