runxinzhi.com
首页
百度搜索
带权的二分匹配
#include <queue> #include <algorithm> /*********************************************************/ //这些都是KM模板 const int N = 305;//二分图中每一个子图的点的最大数 const int INF = 1<<28;//正无穷 bool xckd[N], yckd[N];//在一次DFS中,Xi与Yi是否在交错树上 int n;//点的个案 int edge[N][N];//二维权值信息用矩阵来存储 int xmate[N], ymate[N];//保存匹配结果 int lx[N], ly[N];//Xi与Yi和标号,即解说中的A[]和B[] int slack[N];//松弛量 int prev[N];//? queue<int> Q; bool bfs();//寻找增广路径 void agument(int); int KMMatch();//KM算法 /*********************************************************/ bool bfs() { while(!Q.empty()) { int p = Q.front(), u = p>>1; Q.pop(); if(p&1) { if(ymate[u] == -1) { agument(u); return true; } else { xckd[ymate[u]] = true; Q.push(ymate[u]<<1); } } else { for(int i = 0; i < n; i++) { if(yckd[i])continue; else if(lx[u] + ly[i] != edge[u][i]) { int ex = lx[u] + ly[i] - edge[u][i]; if(slack[i] > ex) { slack[i] = ex; prev[i] = u; } } else { yckd[i] = true; prev[i] = u; Q.push((i<<1)|1); } } } } return false; } void agument(int u) { while(u != -1) { int pv = xmate[prev[u]]; ymate[u] = prev[u]; xmate[prev[u]] = u; u = pv; } } int KMMatch() { int i, j, mn; memset(ly, 0, sizeof(ly)); for(i = 0; i < n; i++) { lx[i] = -INF; for(j = 0; j < n; j++) { lx[i] = lx[i]>edge[i][j]?lx[i]:edge[i][j]; } } memset(xmate, -1, sizeof(xmate)); memset(ymate, -1, sizeof(ymate)); bool agu = true; for(mn = 0; mn < n; mn++) { if(agu) { memset(xckd, 0, sizeof(xckd)); memset(yckd, 0, sizeof(yckd)); for(i = 0; i < n; i++) slack[i] = INF; while(!Q.empty()) Q.pop(); xckd[mn] = true; Q.push(mn<<1); } if(bfs()) { agu = true; continue; } int ex = INF; mn--; agu = false; for(i = 0; i < n; i++) if(!yckd[i]) ex = ex<slack[i]?ex:slack[i]; for(i = 0; i < n; i++) { if(xckd[i]) lx[i] -= ex; if(yckd[i]) ly[i] += ex; slack[i] -= ex; } for(i = 0; i < n; i++) if(!yckd[i] && slack[i] == 0) { yckd[i] = true; Q.push((i<<1)|1); } } int cost = 0; for(i = 0; i < n; i++) cost += edge[i][xmate[i]]; return cost; }
相关阅读:
Codeforces Round #353 (Div. 2)
Codeforces Round #304 (Div. 2)
Codeforces Round #250 (Div. 2)D
Codeforces Round #368 (Div. 2)
hdu4348区间更新的主席树+标记永久化
poj3468线段树标记永久化
Educational Codeforces Round 35 (Rated for Div. 2)
一维数组取一部分
序列化数组。
禅道常识
原文地址:https://www.cnblogs.com/windmissing/p/2559824.html
最新文章
紫书 例题 10-9 UVa 1636 (概率计算)
紫书 例题 10-8 UVa 1262 (暴力枚举)
紫书 例题 10-7 UVa 10820 (欧拉函数)
紫书 例题 10-6 UVa 1635 (二项式定理+唯一分解定理)
紫书 例题 10-5 UVa 12716 (枚举方式)
紫书 例题 10-4 UVa 10791(唯一分解定理)
紫书 例题 10-3 UVa 10375 (唯一分解定理)
Linux的终端与进程
(从终端看linux-1)linux tty pty pts 概念 区别
研究系统IO和glib IO的关系
热门文章
open/write/read
linux文件分类
linux 新建用户、用户组 以及为新用户分配权限
SSH 远程连接
SSH应该使用密钥还是密码?
ssh-keygen的使用方法
加密算法
hdu1695莫比乌斯反演模板题
逆元板子
Codeforces Round #220 (Div. 2)
Copyright © 2020-2023
润新知