【网络流/二分图匹配】飞行员配对

最新推荐文章于 2021-04-29 02:48:02 发布

Kemlkyo

最新推荐文章于 2021-04-29 02:48:02 发布

阅读量771

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/kemlkyo/article/details/24178925

版权

网络流专栏收录该内容

3 篇文章 0 订阅

订阅专栏

简单的二分图匹配问题

从S向每个正驾驶员连一条边，从可以配对的正驾驶员向副驾驶员连边，从每个副驾驶员向T连边

以上都是有向边，容量为1

跑一边最大流就是答案

需要输出方案的话就是找每一条飞行员直接的满流边

/*Author : Soap
  Problem : 网络流24题 01
  Statu : COGS AC 
  Solution : 二分图匹配
  */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 500
#define INF  0x3f3f3f3f
struct Edge{
       int to,next,cap,op;
       }edge[maxn]; 
 
int tot,last[maxn],ans,m,n,S,T,dep[maxn],gap[maxn];
 
void add(int x,int y,int w){
      edge[++tot].to=y;
      edge[tot].next=last[x];
      edge[tot].cap=w;
      edge[tot].op=tot+1;
      last[x]=tot;
      edge[++tot].to=x;
      edge[tot].next=last[y];
      edge[tot].cap=0;
      edge[tot].op=tot-1;
      last[y]=tot;
}
 
void build(){
    scanf("%d%d",&n,&m);
    S=0; T=n+1;
    n+=2;
    int a,b;
    for (int i=1;i<=m;i++) add(S,i,1);
    for (int i=m+1;i<=n;i++) add(i,T,1);
    while (scanf("%d%d",&a,&b)!=EOF){
          if (a==-1 && b==-1) break;
          add(a,b,1);
          }
    }
    
int SAP(int now,int delta){
    if (now==T) return delta;
    int sum=0,mindis=n; 
    for (int i=last[now];i;i=edge[i].next){
        if (edge[i].cap>0 && dep[edge[i].to]+1==dep[now]){
           int save=SAP(edge[i].to,min(edge[i].cap,delta-sum));
           edge[i].cap-=save;
           edge[edge[i].op].cap+=save;
           sum+=save;
           if (sum==delta || dep[S]>=n) return sum;
           }
        if (edge[i].cap>0) mindis=min(mindis,dep[edge[i].to]);
    }
    if (sum==0){
       if (!--gap[dep[now]]) dep[S]=n;
         else ++gap[dep[now]=mindis+1];
         }
    return sum;
}   
           
 
int main(){
    freopen("flyer.in","r",stdin);
    freopen("flyer.out","w",stdout);
    build();
    memset(dep,0,sizeof(dep));
    gap[0]=n;
    while (dep[S]<n) ans+=SAP(S,INF);
    if (ans==0) printf("No Solution!\n");
      else  {
            printf("%d\n",ans);
           /* bool v[maxn];
            memset(v,false,sizeof(v));
            for (int i=1;i<=m;i++)
               for (int j=last[i];j;j=edge[j].next)
                   if (edge[j].cap==0 && !v[edge[j].to] && edge[j].to!=S)
                     {
                        printf("%d %d\n",i,edge[j].to);
                        v[edge[j].to]=true;
                        break;
                        }*/ //输出方案
            }      
}