首页 > 代码库 > 二分图最大权匹配
二分图最大权匹配
#include<string.h> #define MAXN 310 #define inf 1000000000 #define _clr(x) memset(x,-1,sizeof(int)*MAXN) int KM(int m,int n,int mat[][MAXN],int *match1,int *match2) { int s[MAXN],t[MAXN],l1[MAXN],l2[MAXN]; int p,q,i,j,k,ret=0; for(i=0;i<m;i++) { l1[i]=-inf; for(j=0;j<n;j++) l1[i]=mat[i][j]>l1[i]?mat[i][j]:l1[i]; if(l1[i]==-inf) return -1; } for(i=0;i<n;i++) l2[i]=0; _clr(match1); _clr(match2); for(i=0;i<m;i++) { _clr(t); p=0;q=0; for(s[0]=i;p<=q&&match1[i]<0;p++) { for(k=s[p],j=0;j<n&&match1[i]<0;j++) { if(l1[k]+l2[j]==mat[k][j]&&t[j]<0) { s[++q]=match2[j]; t[j]=k; if(s[q]<0) { for(p=j;p>=0;j=p) { match2[j]=k=t[j]; p=match1[k]; match1[k]=j; } } } } } if(match1[i]<0) { i--; p=inf; for(k=0;k<=q;k++) { for(j=0;j<n;j++) { if(t[j]<0&&l1[s[k]]+l2[j]-mat[s[k]][j]<p) p=l1[s[k]]+l2[j]-mat[s[k]][j]; } } for(j=0;j<n;j++) l2[j]+=t[j]<0?0:p; for(k=0;k<=q;k++) l1[s[k]]-=p; } } for(i=0;i<m;i++) ret+=mat[i][match1[i]]; return ret; }
二分图最大权匹配
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。