过山车
Problem Description
RPG girls今天和大家一起去游乐场玩,终于可以坐上梦寐以求的过山车了。可是,过山车的每一排只有两个座位,而且还有条不成文的规矩,就是每个女生必须找个个男生做partner和她同坐。但是,每个女孩都有各自的想法,举个例子把,Rabbit只愿意和XHD或PQK做partner,Grass只愿意和linle或LL做partner,PrincessSnow愿意和水域浪子或伪酷儿做partner。考虑到经费问题,boss刘决定只让找到partner的人去坐过山车,其他的人,嘿嘿,就站在下面看着吧。聪明的Acmer,你可以帮忙算算最多有多少对组合可以坐上过山车吗?
Input
输入数据的第一行是三个整数K , M , N,分别表示可能的组合数目,女生的人数,男生的人数。0<K<=1000
1<=N 和M<=500.接下来的K行,每行有两个数,分别表示女生Ai愿意和男生Bj做partner。最后一个0结束输入。
Output
对于每组数据,输出一个整数,表示可以坐上过山车的最多组合数。
Sample Input
6 3 3
1 1
1 2
1 3
2 1
2 3
3 1
0
Sample Output
3
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <memory.h>
int n,m;
int map[505][505];
int vis[1100];//这轮里是否访问过
int match[10001];//女生匹配的男生
int find(int x)
{
int i;
for(i=1;i<=n;i++)
{
if(map[i][x]&&!vis[i])//女生愿意并且没访问过
{
vis[i]=1;//标记
if(!match[i]||find(match[i]))//没匹配过或匹配的男生还能找到其他女生
{
match[i]=x;
return 1;
}
}
}
return 0;
}
int main()
{
int T,u,v,i;
while(scanf("%d",&T)!=EOF)
{
if(T==0) break;
scanf("%d%d",&n,&m);
memset(map,0,sizeof(map));
memset(match,0,sizeof(match));
while(T--)
{
scanf("%d%d",&u,&v);
map[u][v]=1;
}
int count=0;
for(i=1;i<=m;i++)//循环每一个男生
{
memset(vis,0,sizeof(vis));
if(find(i))//找到配对
count++;
}
printf("%d\n",count);
}
return 0;
}
匈牙利算法几乎是二分图匹配的核心算法,除了二分图多重匹配外均可使用
匈牙利算法实际上就是一种网络流的思想,其核心就是寻找增广路。
“匈牙利算法”指出 通过一条增广路径,通过取反操作,我们就能匹配更多的点。
取反:增广路上的边性质改变,连上的变为可以连的,可以连的变为连上的
什么是增广路径?增广路径是指,由一个未匹配的顶点开始,经过若干个匹配顶点,最后到达对面集合的一个未匹配顶点的路径,即这条路径将两个不同集合的两个未匹配顶点通过一系列匹配顶点相连。
匈牙利算法核心思想:
-
1.初始匹配M为空
-
2.找出一条增广路径p,取反操作得到更大的匹配M’代替M
-
3.重复步骤2,直到找不出增广路为止
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <memory.h>
int n,m;
int map[505][505];
int vis[1100];//这轮里是否访问过
int match[10001];//女生匹配的男生
int find(int u)//就是一个dfs//寻找从u出发的增广路径
{
int v;
for(v=1;v<=n;v++)
{
if(map[v][u]&&!vis[v])
{
vis[v]=1;
if(!match[v]||find(match[v]))//所谓增广路径取反 需要满足w=match[v]的w有新的路径后 将(w,v)改为未匹配,(u,v)设为匹配
{
match[v]=u;
return 1;
}
}
}
return 0;
}
int main()
{
int T,u,v,i;
while(scanf("%d",&T)!=EOF)
{
if(T==0) break;
scanf("%d%d",&n,&m);
memset(map,0,sizeof(map));
memset(match,0,sizeof(match));
while(T--)
{
scanf("%d%d",&u,&v);
map[u][v]=1;
}
int count=0;
for(i=1;i<=m;i++)//循环每一个男生
{
memset(vis,0,sizeof(vis));
if(find(i))//找到配对
count++;
}
for(i=1;i<=n;i++)//循环每一个男生
{
printf("match: %d",match[i]);
}
printf("%d\n",count);
}
return 0;
}
更详尽的解释见二分图匹配——匈牙利算法和KM算法