简单的二分图匹配问题
从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;
}*/ //输出方案
}
}