Luogu 1559 运动员最佳匹配问题(带权二分图最大匹配)
Luogu 1559 运动员最佳匹配问题(带权二分图最大匹配)
Description
羽毛球队有男女运动员各n人。给定2 个n×n矩阵P和Q。P[i][j]是男运动员i和女运动员j配对组成混合双打的男运动员竞赛优势;Q[i][j]是女运动员i和男运动员j配合的女运动员竞赛优势。由于技术配合和心理状态等各种因素影响,P[i][j]不一定等于Q[j][i]。男运动员i和女运动员j配对组成混合双打的男女双方竞赛优势为P[i][j]*Q[j][i]。设计一个算法,计算男女运动员最佳配对法,使各组男女双方竞赛优势的总和达到最大。
Input
第一行有1 个正整数n (1≤n≤20)。接下来的2n行,每行n个数。前n行是p,后n行是q。
Output
将计算出的男女双方竞赛优势的总和的最大值输出。
Sample Input
3
10 2 3
2 3 4
3 4 5
2 2 2
3 5 3
4 5 1
Sample Output
52
Http
Luogu:https://www.luogu.org/problem/show?pid=1559#sub
Source
带权二分图的最大匹配
解决思路
这个题目是一道带权二分图最大匹配的题目,我们要用到的就是KM算法。
那么什么是KM算法呢?
通俗点来讲,KM算法就是给匈牙利算法提供一张图,让后者去匹配。
KM算法是基于贪心的思想,选取权值最大的边建图,然后让匈牙利算法去判断其可行性。
我们知道,匈牙利算法的关键就是“让”出边来匹配,那么KM算法的关键就是如何建出这张图。
接下来我们结合代码来看一看。
首先我们对代码中的变量进行一下说明,我们定义二分图的左边这一组为
X,右边一组为Y。定义use_left,use_right分别标记左边或右边的点是否标记,同时标记其是否在增广路中。定义Wx,Wy分别为左边和右边的顶标这是KM算法的关键之处。这个顶标是用来干嘛的呢?这就是用来判断某一条边是否在我们定义的图中的。我们定义只有满足Wx[i]+Wy[j]==W[i][j]的这样一条边i->j才在图中存在。所以初始化的时候,Wx[i]置为从i出发的最大边的权值,而Wy[j]置为0。这就相当于把i->j这条边加入图中。
对于每一个左边的点i,我们首先对其跑一边匈牙利算法,如果发现无法匹配,则进行下面的操作。从剩余的边中选一条出发点i已标记而到达点j未标记(标记操作在匈牙利算法中实现),并满足Wx[i]+Wy[j]-W[i][j]最小。选择这条边,然后把所有已标记的左边的点的Wx加上W[i][j],而右边的已标记的点的Wy分别减去W[i][j]。这相当于把i->j这条边加入到我们要判断的图中,并且原来图中有的边不变。(仔细想一想,结合我们判断某条边是否在图中的方法:Wx[i]+Wy[j]==W[i][j])
推荐一篇关于二分图匹配的KM算法的文章:http://www.cnblogs.com/Lanly/p/6291214.html
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define mem(Arr,x) memset(Arr,x,sizeof(Arr))
class EDGE
{
public:
int v,w;
};
const int maxN=30;
const int inf=2147483647;
int n;
vector<EDGE> E[maxN];//存二分图的边
int Match[maxN];//存右边匹配的左点编号
bool use_left[maxN];//判断是否访问,一是避免重复,二十判断是否在增广路上,以便后面加边
bool use_right[maxN];
int Wx[maxN];//两边的顶标
int Wy[maxN];
int read();
bool Hungary(int u);
int main()
{
int P[maxN][maxN];
int Q[maxN][maxN];
mem(Match,-1);
mem(Wy,0);
cin>>n;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
P[i][j]=read();
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
Q[i][j]=read();
for (int i=1;i<=n;i++)
{
int D=-inf;//初始化左边的顶标
for (int j=1;j<=n;j++)
{
E[i].push_back((EDGE){j,P[i][j]*Q[j][i]});
D=max(D,P[i][j]*Q[j][i]);//预处理出Wx
}
Wx[i]=D;
Wy[i]=0;
}
for (int i=1;i<=n;i++)
{
do
{
mem(use_left,0);
mem(use_right,0);
if (Hungary(i))
break;
int D=inf;
for (int j=1;j<=n;j++)
if (use_left[j]==1)
{
for (int k=0;k<E[j].size();k++)
{
int v=E[j][k].v;
int w=E[j][k].w;
if (use_right[v]==0)
{
D=min(D,Wx[j]+Wy[v]-w);
}
}
}
for (int j=1;j<=n;j++)
if (use_left[j])
Wx[j]=Wx[j]-D;
for (int j=1;j<=n;j++)
if (use_right[j])
Wy[j]=Wy[j]+D;
}
while (1);
}
int Ans=0;
for (int i=1;i<=n;i++)
for (int j=0;j<E[Match[i]].size();j++)
{
if (E[Match[i]][j].v==i)
{
Ans+=E[Match[i]][j].w;
break;
}
}
cout<<Ans<<endl;
return 0;
}
int read()//读入优化
{
int x=0;
int k=1;
char ch=getchar();
while (((ch>'9')||(ch<'0'))&&(ch!='-'))
ch=getchar();
if (ch=='-')
{
k=-1;
ch=getchar();
}
while ((ch<='9')&&(ch>='0'))
{
x=x*10+ch-48;
ch=getchar();
}
return x*k;
}
bool Hungary(int u)//匈牙利算法判断
{
use_left[u]=1;
for (int i=0;i<E[u].size();i++)
{
int v=E[u][i].v;
if ((use_right[v]==0)&&(Wx[u]+Wy[v]==E[u][i].w))
{
use_right[v]=1;
if ((Match[v]==-1) || (Hungary(Match[v])) )
{
Match[v]=u;
return 1;
}
}
}
return 0;
}
Luogu 1559 运动员最佳匹配问题(带权二分图最大匹配)的更多相关文章
- [Luogu 1559]运动员最佳匹配问题
Description 题库链接 求 \(2\times N\) 个点的带权二分图最佳匹配. \(1\leq N\leq 20\) Solution 我还是太菜了啊...到现在才学 \(KM\) . ...
- HDU 2255 奔小康赚大钱(带权二分图最大匹配)
HDU 2255 奔小康赚大钱(带权二分图最大匹配) Description 传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子. 这可是一件大事,关系到人民的住房问题啊 ...
- UVA1349(带权二分图最大匹配 --> KM算法模板)
UVA1349 题意:给定一些有向带权边,求出把这些边构造成一个个环,总权值最小 解法: 对于带权的二分图的匹配问题可以用通过KM算法求解. 要求最大权匹配就是初始化g[i][j]为0,直接跑就可以: ...
- "《算法导论》之‘图’":不带权二分图最大匹配(匈牙利算法)
博文“二分图的最大匹配.完美匹配和匈牙利算法”对二分图相关的几个概念讲的特别形象,特别容易理解.本文介绍部分主要摘自此博文. 还有其他可参考博文: 趣写算法系列之--匈牙利算法 用于二分图匹配的匈牙利 ...
- KM模板 最大权匹配(广搜版) Luogu P1559 运动员最佳匹配问题
KM板题: #include <bits/stdc++.h> using namespace std; inline void read(int &num) { char ch; ...
- 带权二分图最大匹配KM算法
二分图的判定 如果一个图是连通的,可以用如下的染色法判定是否二分图: 我们把X部的结点颜色设为0,Y部的颜色设为1. 从某个未染色的结点u开始,做BFS或者DFS .把u染为0,枚举u的儿子v.如果v ...
- POJ 2195 Going Home | 带权二分图匹配
给个地图有人和房子 保证人==房子,每个人移动到房子处需要花费曼哈顿距离的代价 问让人都住在房子里最小代价 显然是个带权二分图最大匹配 转化成以一个网络,规定w是容量,c是代价 1.S向人连边,w=1 ...
- 运动员最佳匹配问题 KM算法:带权二分图匹配
题面: 羽毛球队有男女运动员各n人.给定2 个n×n矩阵P和Q.P[i][j]是男运动员i和女运动员j配对组成混合双打的男运动员竞赛优势:Q[i][j]是女运动员i和男运动员j配合的女运动员竞赛优势. ...
- KM(Kuhn-Munkres)算法求带权二分图的最佳匹配
KM(Kuhn-Munkres)算法求带权二分图的最佳匹配 相关概念 这个算法个人觉得一开始时有点难以理解它的一些概念,特别是新定义出来的,因为不知道是干嘛用的.但是,在了解了算法的执行过程和原理后, ...
随机推荐
- MPP 一、Greenplum 集群安装
Installating and Initializing a Greenplum Database System... 1 安装说明 1.1 环境说明 名称 版本 下载地址 虚拟机 Oracle V ...
- docker - 关于network的一些理解
docker 提供给我们多种(4种)网络模式,我们可以根据自己的需求来使用.例如我们在一台主机(host)或者同一个docker engine上面运行continer的时候,我们就可以选择bridge ...
- ReactNative学习之Html基础
前言: React Native开发作为一种新型的移动开发方式,个人觉得App的一部分需求会逐步替换成这种方式,也是公司移动开发人员所必须掌握的一种开发技术,所以鉴于这种情况我觉得很有必要学习一下,特 ...
- WebSocket+MSE——HTML5 直播技术解析
作者 | 刘博(又拍云多媒体开发工程师) 当前为了满足比较火热的移动 Web 端直播需求,一系列的 HTML5 直播技术迅速的发展起来. 常见的可用于 HTML5 的直播技术有 HLS.WebSock ...
- java中File类应用:遍历文件夹下所有文件
练习: 要求指定文件夹下的所有文件,包括子文件夹下的文件 代码: package 遍历文件夹所有文件; import java.io.File; public class Test { public ...
- Python成长之路 — 字典
一.字典的定义与创建 字典是Python中唯一内建的映射类型.你可以将其想象成书本的目录,章节名称代表"key",页码则代表"value".书本的目录本质上是也 ...
- 【转】Delphi多线程编程
文章来源: http://liukun966123.my.gsdn.net/2004/10/22/4797/ Delphi中有一个线程类TThread是用来实现多线程编程的,这个绝大多数Delphi书 ...
- 如何用VS进行程序调试
VS是一个强大的IDE,如果你现在只会简单地用它查看一下执行效果,那就太大材小用了. 1. CRT函数报错 首先来说说最常见的一个编译错误.微信里常常收到这个错误的截图提问. CRT(C Runtim ...
- Chapter 1:Introduction
作者:桂. 时间:2017-05-24 08:06:45 主要是<Speech enhancement: theory and practice>的读书笔记,全部内容可以点击这里. 1. ...
- 通过数据,修改金蝶ERP的收料通知单不合格和合格数量,修改生产投料单,委外发出数量
update POInStockEntry set FAuxNotPassQty=不合格数量 where FInterID=(select FInterID from POInStock where ...