首页 > 代码库 > [kuangbin带你飞]专题十 匹配问题 二分图多重匹配
[kuangbin带你飞]专题十 匹配问题 二分图多重匹配
二分图的多重匹配问题不同于普通的最大匹配中的“每个点只能有最多一条边” 而是“每个点连接的边数不超过自己的限定数量”
最大匹配所解决的问题一般是“每个人都有一群想加入的团体 并且一个团体只能收一个人 问有多少人可以加入一个自己喜欢的团体”
而多重匹配是 “每个人都有一群想加入的团体 每个团体可以收给定的人数 问有多少人可以加入一个自己喜欢的团体”
解决这个问题 目前看貌似有三个办法
1 拆点 一个团体可以招x个人 就把它拆成x个只能招一个人的团体 进行最大匹配
2 网络流 s -> 人 cap : 1 人 -> 自己想加的团体 cap : 1 每个团体 -> t cap : 这个团体的容纳人数
3 改一改 xyl的形式 过去我们写linker[]来记录 每个团体的招人情况
现在我们创造一个数组 记录这个团体 能招多少人 现招多少人 招人的情况 这个网址讲的很简单 http://wenku.baidu.com/view/60b301c00c22590102029db2.html
专题里面的这三道题都差不多 都是用二分求一个最大最小值
M poj2289
模板题 需要做一个字符串处理
#include<stdio.h> #include<string.h> #include<algorithm> #include<math.h> #include<map> #include<string> #include<vector> #include<queue> #include<iostream> using namespace std; #define L long long int n , m ; char s[1050]; struct node{ int many ; int pp[1050]; int cnt ; }linker[505]; bool vis[505]; vector<int >q[1050]; bool fin(int u){ for(int i = 0 ; i< q[u].size() ; i ++ ){ int v = q[u][i]; if(vis[v]){ vis[v] = false ; if(linker[v].cnt < linker[v].many) { linker[v].pp[++linker[v].cnt] = u ; return true ; } else { for(int j = 1; j <= linker[v].many ;j ++ ){ if(fin(linker[v].pp[j])) { linker[v].pp[j] = u ; return true ; } } } } } return false ; } bool xyl(int many) { int res = 0 ; for(int i = 0 ; i < m ; i ++) { linker[i].many = many ; linker[i].cnt = 0 ; } for(int i = 1 ; i <= n ; i ++ ){ memset(vis , true , sizeof(vis )) ; if(fin(i)){ res ++ ; } } return (res == n) ; } int main(){ while(scanf("%d%d\n",&n,&m)!=EOF){ if(n == 0 && m == 0)break; for(int i = 1; i<= n ; i ++)q[i].clear() ; for(int i = 1; i<= n ; i ++ ){ gets(s) ; int len = strlen(s); int j ; for(j = 0 ; j < len ;j ++ ){ if(s[j] == ‘ ‘){ j ++ ; break ; } } int res = 0; for( ; j < len ; j ++ ){ if(s[j ] == ‘ ‘){ q[i].push_back(res ); res = 0 ; } else { res *= 10 ; res += ( s[j] - ‘0‘ ); } } q[i].push_back(res); res = 0; } int ans = -1 ; int l = 1 ; int r = 1010; while(l <= r){ int mid = (l + r) / 2 ; if(xyl(mid)){ ans = mid ; r = mid - 1 ; } else { l = mid + 1 ; } } printf("%d\n",ans); } }
N poj2112
画风突然不对了 细节很多
1 给出的邻接矩阵形式 并不是最短距离 而是类似于 给出一群路来 连接每一个地点
2 给出的 0 代表这两个地点没有路
需要二分一下最大距离
#include<stdio.h> #include<string.h> #include<algorithm> #include<math.h> #include<map> #include<string> #include<vector> #include<queue> #include<iostream> using namespace std; #define L long long int n , a , b , m ; struct node{ int many ; int pp[300]; int cnt ; }linker[300]; bool vis[300]; int c[300][300]; bool fin(int u , int jl){ for(int v = 1; v <= a ; v ++ ){ if(c[u][v] > jl || c[u][v] == 0){ continue ; } if(vis[v]){ vis[v] = false ; if(linker[v].cnt < linker[v].many){ linker[v].pp[++linker[v].cnt] = u ; return true ; } else { for(int j = 1 ; j <= linker[v].many ; j ++ ){ if(fin(linker[v].pp[j] , jl)){ linker[v].pp[j] = u ; return true ; } } } } } return false ; } bool xyl(int jl){ int res = 0 ; for(int i = 1; i <= a ;i ++){ linker[i].cnt = 0 ; linker[i].many = m ; } for(int i = a + 1 ; i <= n ; i ++ ){ memset(vis , true , sizeof(vis)) ; if(fin(i , jl)){ res ++ ; } } return (res == b) ; } int main(){ while(scanf("%d%d%d",&a , &b , & m)!= EOF){ n = a + b; for(int i = 1; i <= n ; i ++ ){ for(int j = 1 ; j <= n ; j ++ ){ scanf("%d" , &c[i][j]) ; } } for(int i = 1 ; i <= n ; i ++ ){ for(int j = 1 ; j <= n; j ++ ){ for(int k = 1; k <= n ; k ++ ){ if(i == j || i == k || k == j)continue; if(c[j][i] != 0 && c[i][k] != 0){ if(c[j][k] > 0) c[j][k] = min(c[j][i] + c[i][k] , c[j][k]) ; else { c[j][k] = c[j][i] + c[i][k] ; } } } } } int l = 1 ; int r = 999999999 ; int ans = -1 ; while(l <= r){ int mid = (l + r) / 2 ; if(xyl(mid)){ ans = mid ; r = mid - 1 ; } else { l = mid + 1 ; } } printf("%d\n",ans) ; } }
O poj3189
比较考验读题能力...
题意是 给出每头牛对每个谷仓的loverank
求出一个最小的范围 让每头牛最后分配的谷仓的loverank都在这个大小范围内
所以二分范围 枚举起点
#include<stdio.h> #include<string.h> #include<algorithm> #include<math.h> #include<map> #include<string> #include<vector> #include<queue> #include<iostream> using namespace std; #define L long long int n , b ; struct node{ int many ; int pp[1050]; int cnt ; }linker[25]; bool vis[25]; int many[25]; int c[1050][25] ; bool fin(int u , int jl , int st){ for(int i = st; i <= st + jl - 1; i ++ ){ int v = c[u][i] ; if(vis[v]) { vis[v] = false ; if(linker[v].cnt < linker[v].many) { linker[v].pp[++linker[v].cnt] = u ; return true ; } else { for(int j = 1 ; j <= linker[v].many ; j ++ ){ if(fin(linker[v].pp[j] , jl , st)){ linker[v].pp[j] = u ; return true ; } } } } } return false ; } bool xyl(int jl , int st){ int res = 0 ; for(int i = 1; i <= b ;i ++ ){ linker[i].cnt = 0; linker[i].many = many[i] ; } for(int i = 1 ; i <= n ; i ++ ){ memset(vis , true , sizeof(vis)) ; if(fin(i , jl , st)){ res ++ ; } } return (res == n) ; } int main(){ while(scanf("%d%d",&n , &b)!=EOF){ for(int i = 1; i<= n ; i ++ ){ for(int j = 1 ; j <= b ;j ++ ){ scanf("%d" , &c[i][j]) ; } } for(int i = 1; i<= b ; i ++ ){ scanf("%d" , &many[i]) ; } int l = 1; int r = b ; int ans = -1 ; while(l <= r){ int mid = (l + r) / 2 ; int i ; for(i = 1; i + mid - 1 <= b ; i ++ ){ if(xyl(mid , i)){ ans = mid ; r = mid - 1 ; break ; } } if(i + mid - 1 > b){ l = mid + 1 ; } } printf("%d\n",ans) ; } }
可以看到Dinic跑多重匹配和跑最大匹配 都是一样的写法 只是cap不同
看起来的话 拆点法和xyl的时间复杂度都是很高的 .. Dinic不错的样子
看起来 小虎的担子 又重了一分
[kuangbin带你飞]专题十 匹配问题 二分图多重匹配