标签:二分图最大权匹配
2 100 10 15 23
123
#include <iostream> #include <cstring> #include <cstdio> using namespace std; const int maxn = 310; const int inf = 1<<30; int weight[maxn][maxn]; //记录完全二分图的所有边的权值 int match[maxn];//记录完全二分图的完美匹配 bool visx[maxn], visy[maxn];//记录x,y集合中的顶点是否出现在增广路中 int lx[maxn], ly[maxn];//记录x,y各顶点的标记量 int slack[maxn];//记录x在增广路,y不在增广路的每条边lx[x]+ly[y]-weight[x][y]到y的最小值放在slack[y] //用来扩大相等子图(二分图的子图)。 bool dfs(const int& x,const int& n) { visx[x] = 1; //x在增广路中出现 for (int i = 0; i < n; i ++) { if (visy[i]) continue; int t = lx[x]+ly[i]-weight[x][i]; if (t == 0) { //x在增广路,y在增广路 visy[i] = 1; if (match[i] == -1 || dfs(match[i], n)) { match[i] = x; return true; } } else { //x在增广路,y不在增广路 slack[i] = min(slack[i], t); } } return false; } int KM(const int& n) { memset(match, -1, sizeof(match)); //初始化,二分图的最大匹配为0 for (int i = 0; i < n; i ++) { lx[i] = 0, ly[i] = 0; for (int j = 0; j < n; j ++) { lx[i] = max(lx[i], weight[i][j]);//初始化lx为x的邻接边最大权值,ly为0 } } for (int i = 0; i < n; i ++) { for (int j = 0; j < n; j ++) slack[j] = inf; while (1) { memset(visx, 0, sizeof(visx)); //初始化x,y都不在增广路中 memset(visy, 0, sizeof(visy)); if (dfs(i, n)) break; //匹配成功,则继续求下个x点的最大匹配边,失败更新各点的标量,增加可行边的数量 int d = inf; for (int j = 0; j < n; j ++) { //求起点在增广路,而终点不在增广路的边,端点标量和与边权差的最小值。 if (!visy[j]) d = min(d, slack[j]); } for (int j = 0; j < n; j ++) { if (visx[j]) lx[j] -= d; } //更新增广路中x的标量 for (int j = 0; j < n; j ++) { if (visy[j]) ly[j] += d; //更新增广路中y的标量 else slack[j] -= d; //更新slack的最小值。(因为与Y中j相连的增广路x的标量都减小了d) } } } int res = 0; for (int i = 0; i < n; i ++) { //Y集合中可能有没有匹配的点 if (match[i] > -1)res += weight[match[i]][i]; } return res; } int main() { int n; while (scanf("%d", &n)!=EOF) { for (int i = 0; i < n; i ++) { for (int j = 0; j < n; j ++) { scanf("%d", &weight[i][j]); } } int ans = KM(n); printf("%d\n", ans); } return 0; }
标签:二分图最大权匹配
原文地址:http://blog.csdn.net/wuhuajunbao/article/details/41978583