Luogu 2756
題意
給出n個英國飛行員和m個外籍飛行員,給出每個英國飛行員能配合的外籍飛行員編號,求最多可以選出多少對能互相配合的飛行員
題解
最經典的二分圖匹配問題。其實用匈牙利算法更合適,在時間複雜度上會更加優秀。用網絡流來做的話,超級源點向每個英國飛行員連流量爲1的邊,每個英國飛行員向能配合的外籍飛行員連流量爲1的邊,每個外籍飛行員向超級匯點連一條流量爲1的邊,這樣保證每個人僅會被選一次。跑一邊最大流即可。
代碼
最好不用該代碼當板子......
#include<cstdio>
#include<cstring>
#include<vector>
#define MAXN 20000+30
#define INF 0x3fffffff
using namespace std;
int m,n,K,Du[MAXN],dis[MAXN],dl[23333];
vector <int> f[MAXN],FB[MAXN],Flow[MAXN];
int cur[MAXN];
int a[30][30];
int Load[30],Num[30];
void ADD(int x,int y,int z)
{
Du[x]++;Du[y]++;
f[x].push_back(y);Flow[x].push_back(z);FB[x].push_back(Du[y]-1);
f[y].push_back(x);Flow[y].push_back(0);FB[y].push_back(Du[x]-1);
}
bool BFS(int Begin,int End)
{
int t=0,w=1,x,X;
memset(dis,0xff,sizeof(dis));
dis[Begin]=0;dl[1]=Begin;
do
{
x=dl[++t];
for(int i=0;i<Du[x];i++)
{
X=f[x][i];
if(Flow[x][i]<=0||dis[X]>=0) continue;//不連通或已訪問
dis[X]=dis[x]+1;dl[++w]=X;
}
}while(t<w);
if(dis[End]>0) return true;
else return false;
}
int Find(int x,int MFLOW,int y)
{
if(x==y) return MFLOW;
int X,h;
for(int i=cur[x];i<Du[x];i++)
{
cur[x]=i;
X=f[x][i];
if(Flow[x][i]>0&&dis[x]+1==dis[X]&&(h=Find(X,min(MFLOW,Flow[x][i]),y)))
{
Flow[x][i]-=h;
Flow[X][FB[x][i]]+=h;
return h;
}
}
return 0;
}
int Solve(int x,int y)
{
int X,Ans=0;
while(1)
{
if(!BFS(x,y)) break;
memset(cur,0,sizeof(cur));
while((X=Find(x,0x7fffffff,y)))
Ans+=X;
}
return Ans;
}
int S,T;
int main()
{
scanf("%d%d",&n,&m);
S=0;T=m+1;
for(int i=1;i<=n;i++)
ADD(S,i,1);
for(int i=n+1;i<=m;i++)
ADD(i,T,1);
int x,y;
while(1)
{
scanf("%d%d",&x,&y);
if(x==-1&&y==-1) break;
ADD(x,y,1);
}
printf("%d\n",Solve(S,T));
for(int i=1;i<=n;i++)
{
for(int j=0;j<Du[i];j++)
{
if(f[i][j]>n&&Flow[i][j]==0) printf("%d %d\n",i,f[i][j]);
}
}
return 0;
}