题目链接:https://vjudge.net/contest/366917#problem/K
题目描述:
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
知识点
- 匹配:两两不含公共端点的边集合
- 二分图: 给定2组顶点,每一组内的任意两个顶点间没有边相连,只有2个集合之间存在边,这样构建出来的
如何求最大匹配?
交错路径:给定图G的一个匹配M,一条路径的边交替出现在M中和不出现在M中。
增广路径:一个交错路径,它的2个端点都不与M中的边关联。
已知匹配:2和4。可以得到一条增广路径1-2-3-4-5。同时可以得到2 3 5为未匹配的
可以推出:每找到一条增广路径,匹配数加1
求最大匹配可以转化为:寻找增广路径
最大匹配的匈牙利算法(双向图)
int V;//定点数
vector<int>G[MAX_V];//图的邻接表表示
int match[MAX_V];//所匹配的顶点
bool used[MAX_V];//dfs中用到的访问标记
void add_edge(int u,int v)//向图中的增加一条连接u和v的边
{
G[u].push_back(v);
G[v].push_back(u);
}
//通过dfs寻找增广图
bool dfs(int v)
{
used[v]=true;
for(int i=0;i<G[v].size();i++)
{
int u=G[v][i],w=match[u];
if(w<0||!used[w]&&dfs(w))//没被别人匹配上,或者原来匹配上的给腾了地方,那当前的就可以匹配上了
{
match[v]=u;
match[u]=v;
return true;
}
}
return false;
}
int bipartite_matching()
{
int res=0;
memset(match,-1,sizeof(match));
for(int v=0;v<V;v++)
{
if(match[v]<0)
{
memset(used,0,sizeof(used));
if(dfs(v))
res++;
}
}
}
对if(w<0||!used[w]&&dfs(w))
解析
我们的目标是尽可能给x中的点找到配对
- 刚开始,一个匹配都没有,我们随意选取一条边,(x1, y1)这条边,构建最初的匹配出来,结果如下,已经配对的边用粗线标出
- . 我们给x2添加一个匹配,如下图的(x2, y2)边
- 我们现在想给x3匹配一条边,发现它的另一端y1已经被x1占用了。让y1离开x1。即将被迫分手的x1很委屈,好在它还有其他的选择,于是 x1 妥协了,x1准备去找自己看中的y2。让 y2主动离开 x2 。x2去找了y5。
牵涉到的节点:(x3, y1, x1, y2, x2, y5),很明显,这是一条路径P
完整代码:
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int N=5*1e2+10;
vector<int>G[N];
int match[N];
bool used[N];
int k,m,n;
bool dfs(int v)
{
for(int i=0; i<G[v].size(); i++)
{
int to=G[v][i];
int w=match[to];
if(!used[to])
{
used[to]=true;
if(w<0||dfs(w))
{
match[to]=v;
return true;
}
}
}
return false;
}
int main()
{
while(~scanf("%d",&k))
{
if(k==0)
break;
for(int i=1; i<N; i++)
G[i].clear();
scanf("%d%d",&m,&n);
while(k--)
{
int u,v;
scanf("%d%d",&u,&v);
G[u].push_back(v);
}
int res=0;
memset(match,-1,sizeof(match));
for(int i=1; i<=m; i++)
{
memset(used,false,sizeof(used));
if(dfs(i))
res++;
}
printf("%d\n",res);
}
return 0;
}