标签:
Description
Input
Output
Sample Input
3 2 John 0 1 Rose 1 Mary 1 5 4 ACM 1 2 3 ICPC 0 1 Asian 0 2 3 Regional 1 2 ShangHai 0 2 0 0
Sample Output
2 2
我的第一道二分图多重匹配题:
转一个不错的Blog:http://www.cnblogs.com/zhengguiping--9876/p/4728358.html
Jamie有很多联系人,但是很不方便管理,他想把这些联系人分成组,已知这些联系人可以被分到哪个组中去,而且要求每个组的联系人上限最小,即有一整数k,使每个组的联系人数都不大于k,问这个k最小是多
少?
一对多的二分图的多重匹配。二分图的多重匹配算法的实现类似于匈牙利算法,对于集合x中的元素xi,找到一个与其相连的元素yi后,检查匈牙利算法的两个条件是否成立,若yi未被匹配,则将
xi,yi匹配。否则,如果与yi匹配的元素已经达到上限,那么在所有与yi匹配的元素中选择一个元素,检查是否能找到一条增广路径,如果能,则让出位置,让xi与yi匹配。
二分求出limit,知道找到可以构成多重匹配的最小限制limit,在main函数中二分搜索。
#include<stdio.h> #include<string.h> #include<queue> #include<algorithm> using namespace std; #define N 1010 int vis[N], maps[N][N], ans, n, m; struct node { int cnt;///和yi相匹配的个数; int k[N];///和yi相匹配的x的集合; }Linky[N]; bool Find(int u, int limit) { for(int i=1; i<=m; i++) { if(!vis[i] && maps[u][i]) { vis[i] = 1; if(Linky[i].cnt < limit) { Linky[i].k[ Linky[i].cnt++ ] = u; return true; } for(int j=0; j<Linky[i].cnt; j++) { if(Find( Linky[i].k[j], limit )) { Linky[i].k[j] = u; return true; } } } } return false; } bool hungary(int limit)///匈牙利算法; { memset(Linky, 0, sizeof(Linky)); for(int i=1; i<=n; i++) { memset(vis, 0, sizeof(vis)); if(!Find(i, limit))///当前的limit让i没有匹配,所以不能用limit; return false; } return true; } int main() { int x; char s[20], ch; while(scanf("%d %d", &n, &m), m+n) { memset(maps, 0, sizeof(maps)); for(int i=1; i<=n; i++) { scanf("%s", s); while(1) { scanf("%d%c", &x, &ch); maps[i][x+1] = 1; if(ch == '\n') break; } } int L = 1, R = n; ans = n; while(L <= R) { int mid = (L+R)/2; if(hungary(mid))///如果当前mid满足题意; { R = mid-1; ans = mid; } else L = mid+1; } printf("%d\n", ans); } return 0; }
Jamie's Contact Groups(二分图多重匹配+二分)(网络流)
标签:
原文地址:http://blog.csdn.net/u014665013/article/details/51344918