二分图最佳完美匹配——KM算法总结
2017-06-03 21:14
225 查看
KM 算法
求解二分图最佳完美匹配的算法。先来看一道例题Hdu 2255。
显然是KM的裸题。假设我们要匹配集合X和Y的点,先给每个点一个顶标Lx和Ly。
为什么要给顶标?
首先顶标是我们限制边的一个条件,我们走过的边一道要满足dst[i][j]==Lx[i]+Ly[j],其次所有顶标的加和就是我们答案,我们肯定希望顶标的加和越大越好,所以Lx的初值给的就是连向Y集合所有边的最大值,那么Ly先给0。
然后我们希望每个人都有一套房子,所以每个人枚举过去,用DFS验证能不能给这个人一套房子,如果能,那太好了,如果不能,我们就要让顶标减小一点点,放宽限制,使一些没房子的人有房子。那么减小多少呢,我们不能减少太多,因为答案要大,但是太小又有可能没有什么效果,所以我们将顶标减小到刚好能多容纳至少一条边,这个值显然就是匹配不成功的边所相差的期望值的最小值(顶标可以看成一种对答案的期望,这个值显然在寻找答案的时候就可以刷出,我们假设它为mi)。还要注意寻找答案时用S数组标记X的集合的节点有没有走过,T数组标记Y。修改顶标时我们这样做:
for (int j=1;j<=n;j++){ if (S[j]) Lx[j]-=mi; if (T[j]) Ly[j]+=mi; }
不难发现这样做不会影响到已经匹配好的边,并且至少会放宽一条新边的限制,而且这样使答案减小的幅度也是最小的。
看到现在,这个算法的复杂度也不难知道,最多放宽限制M次,在乘以寻找匹配点find的复杂度。我们可以先认为这是O(n*M)的,M=n^2,所以复杂度为O(n^3),但是略比这个大。KM还有一种就是用BFS实现,非常强大,但是……我不会QAQ!!!!!!!
在考虑一些特殊的情况,就是题目中给出的图并不能完美匹配,这样我们的程序会被卡死,所以可以加入一些边权为0的新边,如果Y集合的点都不够,我们就要加入新的虚拟节点。这样题目的图在不影响答案的条件下变完美了 >_<。
#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int maxn=405,MAXINT=2147483647; int dst[maxn][maxn],n1,n2,m,n,Lx[maxn],Ly[maxn],mi,who[maxn],w[maxn]; bool S[maxn],T[maxn]; bool find(int x){ S[x]=1; for (int i=1;i<=n;i++) if (!T[i]){ int s=Lx[x]+Ly[i]-dst[x][i]; if (!s){T[i]=1;if (!who[i]||find(who[i])){who[i]=x;return 1;}} else mi=min(mi,s); } return 0; } void KM(){ for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) Lx[i]=max(Lx[i],dst[i][j]); for (int i=1;i<=n;i++){ while (1){ memset(S,0,sizeof(S));memset(T,0,sizeof(T)); mi=MAXINT; if (find(i)) break; for (int j=1;j<=n;j++){ if (S[j]) Lx[j]-=mi; if (T[j]) Ly[j]+=mi; } } } long long ans=0; for (int i=1;i<=n;i++) ans+=Lx[i]+Ly[i],w[who[i]]=i; printf("%lld\n",ans); for (int i=1;i<=n1;i++) if (dst[i][w[i]]) printf("%d ",w[i]);else printf("0 "); } inline int _read(){ int num=0;char ch=getchar(); while (ch<'0'||ch>'9') ch=getchar(); while (ch>='0'&&ch<='9') num=num*10+ch-48,ch=getchar(); return num; } int main(){ freopen("exam.in","r",stdin); freopen("exam.out","w",stdout); n1=_read(),n2=_read(),m=_read();n=max(n1,n2); for (int i=1;i<=m;i++){int x=_read(),y=_read();dst[x][y]=_read();} KM(); return 0; }
相关文章推荐
- Uva 1411 ants(KM算法--二分图最佳完美匹配)
- 二分图的最佳完美匹配——KM算法
- KM算法(二分图的最佳完美匹配)
- 完美匹配【KM算法-hdu总结】
- 二分图带权匹配、最佳匹配与KM算法
- UVALive 4043 Ants 蚂蚁(二分图最佳完美匹配,KM算法)
- UVa1349 Optimal Bus Route Design(二分图最佳完美匹配)
- 【转】二分图的完美匹配(最佳匹配)之KM算法
- KM算法 二分图的最佳匹配
- KM算法二分图完美匹配 模板
- 1076. Trash(KM算法 二分最佳完美匹配)
- [BZOJ1070] [SCOI2007]修车 && 二分图最佳完美匹配 或 费用流
- 二分图最佳完美匹配
- UVA 11383 Golden Tiger Claw(最佳二分图完美匹配)
- HDU-#2255 奔小康赚大钱(二分图最佳完美匹配+KM)
- 二分图的最佳完美匹配(和一些变形…
- 二分图带权匹配、最佳匹配与KM算法
- hdu - 3488 - Tour(二分图最佳完美匹配)
- Ural 1076 Trash (二分图最佳完美匹配)
- UVA - 1045 The Great Wall Game(二分图最佳完美匹配)