二分图匹配(一)

文章目录

  • 什么是二分图:
  • 例题:
    • NC111768 CF741C
      • 题目描述:
      • 题解:
      • 代码:
  • 二分图最大匹配
    • 匈牙利算法
      • 算法思想:
      • 代码:
    • König定理
  • 二分图最优匹配
    • KM(Kuhn-Munkres)算法
      • 算法思路:
      • 具体操作
      • 代码:
      • 复杂度
    • 代码:

什么是二分图:

二分图充分必要条件条件:
至少有两个顶点且没有奇环

二分图判断:
通过黑白染色

例题:

NC111768 CF741C

题目描述:

有n对情侣(2n个人)围成一圈坐在桌子边上,每个人占据一个位子,要求情侣不能吃同一
种食物,并且桌子上相邻的三个人的食物必须有两个人是不同的,只有两种食物(1或者是2),问一种可行分配方式。

题解:

如果我们能把不能吃同一种食物的人连边,问题就变成二分图黑白染色
• 所以情侣关系等价于两者之间连一条边
• “每连续的三个人不能都一样”怎么办?
• 让第2i个人和第2i+1个人不能吃一样的食物即可(即1连2,3连4,5连6以此类推)
• 这样肯定是个二分图——2i和2i-1分别连了他两的情侣,情侣又分别连他两的一个邻居……
每次都是给这个可能存在的环加两个点,所以有环就一定不是奇环

构造方法如下:首先情侣连边,然后在原图的相邻点对上,隔一对连一条边

代码:

二分图最大匹配

促成更多的点配对

匈牙利算法

有一场宴会,男孩和女孩组成舞伴,并且他们必须互相喜欢才能成为舞伴,一个男孩可能喜欢0个或多个女孩,一个女孩也可能喜欢0个或多个男孩,但一个男孩和他喜欢地女孩成为舞伴之后就不能和其他他喜欢地女孩成为舞伴,女孩亦是如此。请问最多可以有多少对舞伴。

算法思想:

对于一个男孩子x,如果他喜欢女孩y,且y还没有舞伴——让他们配对
• 如果y有了舞伴,x还是会去尝试邀请y,如果y发现她的舞伴z可以换一个舞伴,y就主动抛弃掉z(在确定z可以和别人牵手之后),和x成为舞伴

增广路(交错路)
路径的起点和终点都是还没有匹配过的点,并且路径经过的连线是一条没被匹配、一条已经匹配过,再下一条又没匹配这样交替地出现。找到这样的路径后,显然路径里没被匹配的连线比已经匹配了的连线多一条,于是修改匹配图,把路径里所有匹配过的连线去掉匹配关系,把没有匹配的连线变成匹配的,这样匹配数就比原来多1个。

代码:

时间复杂度:
邻接矩阵O(n^3)
邻接表:O(n*m)

bool dfs(int x) 
{for(int i=1;i<=n;i++)//遍历女生 {if(!a[x][i]||vis[i])continue;//没连线或者已被找过 vis[i]=1;if(link[i]==0||dfs(link[i]))//link[i]表示第i个女生的男伴是谁 {link[i]=x;return 1;}}return 0;
}int main()
{for(int i=1;i<=n;i++)//遍历男生 {memset(vis,0,sizeof(vis));//vis[i]第i个女生在这次找搭档的过程中被邀请 sum+=dfs(i);}
}

König定理

二分图中最大匹配数等于这个图中的最小点覆盖数

最小点覆盖:
每个点覆盖以它为端点的所有边,选择最少的点来覆盖所有的边
(选最少的人,联系到其他的所有的人)
证明:
视频时间:01:03:42

二分图最优匹配

存在边权,问最大匹配下的最大边权

KM(Kuhn-Munkres)算法

最优匹配是建立在完美匹配的基础上的,如果不存在完美匹配,那么本算法失效(但是,我们可以人为连一些权值为0的边,甚至加点,使得没有匹配的节点们最后都有一个“虚假”的匹配)。

算法思路:

最开始将每个左边节点连的权值最大的边视为有效边,在匹配给过程当中如果某个点找不到匹配点,则选择将一些无效边改成有效边继续去匹配,选择的这些无效边其实是换掉的原来的某条有效边,那么我们肯定选换掉的代价最小的。
先让所有人选最佳搭档,当发生冲突时,让降低标准最少的人来换搭档

具体操作

给每个点预设一个顶标,只有两个端点的顶标加起来等于边权的边,我们才认为是有效边,
即L[x]+R[y]==w[x][y]的时候才是有效边
• 最开始左集合的每个点的顶标为他连出去权值最大的边的权值,右集合每个点顶标为0,当匹配失败的时候,我们遍历之前的左集合本次尝试匹配遍历过的点,在他们的连向右集合未遍历的点的边中找一个与顶标和差值最小的边,即delta=w[x][y]-L[x]-R[y]最小的边(可理解为找出一条损失最小的找出一条增广路。)找到之后修改顶标:左侧所有遍历过的点减去
delta,右边所有遍历过的点加上delta,这样原来的有效边还是有效边,而我们新加的边也已经加上了。
• 重复找匹配+修改顶标的操作,一直到找到一个完美匹配即可
在这里插入图片描述
详细过程:
降低的标准x,左边-x,右边+x
在这里插入图片描述
每个人都要为了团队牺牲自己,甘愿降低标准

代码:

lx[]//左端点,赋值为0x7f 
ly[]//右端点,赋值为0
link[i]=j//第i个人匹配的是j 
int dfs(int x)
{visx[x]=1; for(int i=1;i<=m;i++){if(!visy[i]&&lx[x]+ly[i]==w[x][i])//有边且为有效边 {visy[i]=1;if(link[i]==-1||dfs(link[i])){link[i]=t;return 1;}}}return 0;
}
memset(ly,0,sizeof(ly));
memset(lx,0xf7,sizeof(lx));
memset(link,-1,sizeof(link));
for(int i=1;i<=n;i++)
{for(int j=1;j<=n;j++){lx[i]=max(w[i][j],lx[i]);//最心动的边 }
} 
for(int i=1;i<=n;i++)
{while(1)//全部找到完匹配才能结束 {memset(visx,0,sizeof(visx)); memset(visy,0,sizeof(visy)); if(dfs(i))break;//已经给第i人找到匹配int d=0x7f7f7f7f; for(int j=1;j<=n;j++){if(visx[j])//要降低标准的男生 {for(int k=1;k<=m;k++) {if(!visy[k])//引入新的女生 d=min(d,lx[j]+ly[k]-w[j][k]);//求最小的代价 }}}for(d==0x7f7f7f7f)return -1;//匹配失败 for(int j=1;j<=n;j++)if(visx[j])lx[j]-=d;for(int j=1;j<=m;j++)if(visy[j])ly[j]+=d;}
}
//n<=m 
//保证所有男生都匹配(左点),女生不一定 

复杂度

在这里插入图片描述
·如果卡n^2*m可以用网络流来做

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#define N 505
#define M 250005
#define INF 9990365505
#define ll long long
using namespace std;
int n,m,x,y,z,tot,tim,l,r,q[N],fr[N],nxt[M],d1[M],d2[M];
int pre[N],visx[N],visy[N],mchx[N],mchy[N];
ll ex[N],ey[N],slack[N];
void add(int x,int y,int z)
{tot++;d1[tot]=y;d2[tot]=z;nxt[tot]=fr[x];fr[x]=tot;
}
void modify(int cur)
{for (int x=cur,lst;x;x=lst)lst=mchx[pre[x]],mchx[pre[x]]=x,mchy[x]=pre[x];
}
void bfs(int cur)
{for (int i=1;i<=n;i++)slack[i]=INF,pre[i]=0;l=1,r=0;q[++r]=cur;++tim;for (;;){while (l<=r){int u=q[l];l++;visx[u]=tim;for (int i=fr[u];i;i=nxt[i]){int v=d1[i],cost=d2[i];if (visy[v]==tim)continue;ll del=ex[u]+ey[v]-cost;if (del<slack[v]){slack[v]=del;pre[v]=u;if (!del){visy[v]=tim;if (!mchy[v]){modify(v);return;}q[++r]=mchy[v];}}}}ll del=INF;for (int i=1;i<=n;i++)if (visy[i]!=tim)del=min(del,slack[i]);l=1,r=0;for (int i=1;i<=n;i++){if (visx[i]==tim)ex[i]-=del;if (visy[i]==tim)ey[i]+=del; elseslack[i]-=del;}for (int i=1;i<=n;i++)if (visy[i]!=tim && !slack[i]){visy[i]=tim;if (!mchy[i]){modify(i);return;}q[++r]=mchy[i];}}
}
void KM()
{for (int i=1;i<=n;i++)bfs(i);ll ans=0;for (int i=1;i<=n;i++)ans+=ex[i]+ey[i];printf("%lld\n",ans);for (int i=1;i<=n;i++)printf("%d ",mchy[i]);putchar('\n');
}
int main()
{scanf("%d%d",&n,&m);for (int i=1;i<=n;i++)ex[i]=-INF,ey[i]=0;for (int i=1;i<=m;i++){scanf("%d%d%d",&x,&y,&z);add(x,y,z);ex[x]=max(ex[x],(ll)z);}KM();return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318645.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

不止代码:恐狼后卫(ybtoj-区间dp)

文章目录题目描述解析代码thanks for reading!题目描述 一代炉石的眼泪啊 解析 用dp[i][j]表示i与j之间的全部消掉&#xff08;不含两端&#xff09;的最小花费 然后枚举中间最后杀死的狼就行了 本题没有一次AC&#xff0c;因为一开始dp定义成了包含两端&#xff0c;然后因为…

牛客题霸 [二叉树的之字形层序遍历] C++题解/答案

牛客题霸 [二叉树的之字形层序遍历] C题解/答案 题目描述 给定一个二叉树&#xff0c;返回该二叉树的之字形层序遍历&#xff0c;&#xff08;第一层从左向右&#xff0c;下一层从右向左&#xff0c;一直这样交替&#xff09; 例如&#xff1a; 给定的二叉树是{3,9,20,#,#,15…

【.NET Core项目实战-统一认证平台】第九章 授权篇-使用Dapper持久化IdentityServer4...

上篇文章介绍了IdentityServer4的源码分析的内容&#xff0c;让我们知道了IdentityServer4的一些运行原理&#xff0c;这篇将介绍如何使用dapper来持久化Identityserver4&#xff0c;让我们对IdentityServer4理解更透彻&#xff0c;并优化下数据请求&#xff0c;减少不必要的开…

不止代码:生日欢唱(ybtoj-区间dp)

文章目录题目描述解析代码thanks for reading&#xff01;题目描述 解析 这题挺好的 思路&#xff1a;dp[i][j]表示必须把i和j配对&#xff0c;可达到的最大值 首先&#xff1a; dp[i][j]dp[i-1][j-1]a[i]*b[j];然后可以分别尝试把男生或女生往前放弃一段&#xff1a; for(i…

打造自己的.NET Core项目模板

前言每个人都有自己习惯的项目结构&#xff0c;有人的喜欢在项目里面建解决方案文件夹&#xff1b;有的人喜欢传统的三层命名&#xff1b;有的人喜欢单一&#xff0c;简单的项目一个csproj就搞定。。反正就是萝卜青菜&#xff0c;各有所爱。可能不同的公司对这些会有特定的要求…

超好用的C#控制台应用模板

默认模板之缺在工作学习中&#xff0c;我们经常需要创建一些简单的控制台应用(Console App)去验证某个想法&#xff0c;或者作为小工具交付给其他同事。通常我们的选择是 Visual Studio 自带的 Console App 模板&#xff0c;这个经典模板只有预设好的 csproj 文件和空荡荡的 Ma…

不止代码:区间圆数(ybtoj-数位DP)

文章目录题目描述解析理解万岁&#xff01;代码题目描述 解析 一寸山河一寸血 理解万岁&#xff01; 首先&#xff0c;这题统计[l,r]的个数&#xff0c;可以用[1,r]-[1,l-1]来实现 接下来就是如何统计出[1,n]的个数了 首先&#xff0c;用dp[pos][s0]来表示一个二进制pos位有s…

.NET Core 2.2 新增部分功能使用尝鲜

前言美国当地时间12月4日&#xff0c;微软2019开发者大会中发布了一系列的重磅消息&#xff0c;包含了软硬件和开源社区的各种好消息是铺天盖地&#xff0c;作为一名普通的开发者&#xff0c;我第一时间下载了 .NET Core 2.2 的源码&#xff0c;针对发布说明逐条浏览&#xff0…

不止代码:机器分配(动态规划)

题目描述 解析 头疼 什么破题 就是一个dp寻找最优性决策的常规题 但是要输出过程&#xff0c;可以使用递归输出 但是&#xff01; 这题数据的意思是&#xff1a;存在a[i]a[i1]的情况&#xff0c;且在不影响利润的情况下&#xff0c;机器要尽可能的用完 这河里吗 可能是机器多了…

2021中国大学生程序设计竞赛(CCPC)- 网络选拔赛 GCD on Sequence(线段树)

https://blog.csdn.net/qq_45863710/article/details/120121607 对于每个左端点 l 维护尚未确定 v 值的最大的右端点 mx[l] #include<bits/stdc.h>using namespace std; using lllong long; const int N100010; vector<int> vec[N]; int n,a[N],p[N]; ll ans[N]; …

Docker最全教程——从理论到实战(六)

本篇教程持续编写了3个星期左右并且一直在完善、补充具体的细节和实践&#xff0c;预计全部完成需要1到2个月的时间。由于编写的过程中极其费时&#xff0c;并且还需要配合做一些实践&#xff08;有些实践存在一些坑&#xff0c;而且极其费时费事&#xff09;。因此目前产出的速…

.NET Core实战项目之CMS 第八章 设计篇-内容管理极简设计全过程

写在前面上一篇文章.NET Core实战项目之CMS 第七章 设计篇-用户权限极简设计全过程中我带着大家进行了权限部分的极简设计&#xff0c;也仅仅是一个基本的权限设计。不过你完全可以基于这套权限系统设计你的更复杂的权限系统&#xff0c;当然更复杂的权限系统要根据你的业务来进…

2021“MINIEYE杯”中国大学生算法设计超级联赛(2)I love max and multiply(转化)

I love max and multiply Code 代码抄的std #include<bits/stdc.h> using namespace std; using lllong long; template <class Tint> T rd() {T res0;T fg1;char chgetchar();while(!isdigit(ch)) {if(ch-) fg-1;chgetchar();}while( isdigit(ch)) res(res<&…

人工智能第二课:认知服务和机器人框架探秘

这是《人工智能系列笔记》的第二篇&#xff0c;我利用周六下午完成课程学习。这一方面是因为内容属于入门级&#xff0c;并且之前我已经对认知服务和机器人框架比较熟悉。如有兴趣&#xff0c;请关注该系列 https://aka.ms/learningAI 但是学习这门课程还是很有收获&#xff0c…

粉刷木板(ybtoj-单调队列)

题目描述 解析 头疼 定义dp[i]:只用前i块板的最大价值 对于新加入的一个木匠&#xff1a; 不难写出dp转移式&#xff1a; dp[i]max(dp[k](i-k)*p)k表示开始刷的前一个 其中i>s i-k<l; 要是这么转移会是n^2m 所以我就不废了。。。 其实离正解很接近了 把上面移一下项&am…

使用Dapper持久化IdentityServer4

最近研究dotnet core,微软将IdentityServer4作为推荐的服务授权和验证的组件,其独立性特别适合微服务或者分布式的服务扩展验证,所以非常受广大dotnet开发人员的青睐.默认的IdentityServer4默认使用内存对象的验证和授权,而在IdentityServer的官方推荐只有Entity Framework cor…

再不学习我们就out了

前不久我们组来了个Graduate Developer&#xff0c;刚毕业&#xff0c;经验不多&#xff0c;有一次闹了个乌龙&#xff0c;把Stage数据库直接删掉了……好在Azure有备份&#xff0c;不然就算Stage没有重要数据&#xff0c;也得花点时间重建&#xff0c;其他的开发、测试工作都得…

矩阵快速幂一篇通

文章目录概述快速幂解析代码矩阵运算定义加法乘法单位矩阵一、斐波拉契&#xff08;基础模板&#xff09;题目描述解析代码二、行为方案&#xff08;实际应用&#xff09;题目描述解析代码三、矩阵求和&#xff08;子矩阵作为矩阵元素&#xff09;题目描述解析代码四、最短路径…

玉米田(加加强版)【插头dp】

前言 水解警告&#xff0c;数据水勉强卡过的 正题 题目大意 n∗mn*mn∗m的网格里面有些格子被禁止&#xff0c;现在求选取若干个不相邻的格子的方案数。 1≤n≤120,1≤m≤211\leq n\leq 120,1\leq m\leq 211≤n≤120,1≤m≤21 解题思路 听说是插头dpdpdp然后想了一下觉得比插…

牛客题霸 [将字符串转化为整数] C++题解/答案

牛客题霸 [将字符串转化为整数] C题解/答案 题目描述 实现函数 atoi 。函数的功能为将字符串转化为整数 提示&#xff1a;仔细思考所有可能的输入情况。这个问题没有给出输入的限制&#xff0c;你需要自己考虑所有可能的情况。 题解&#xff1a; 题目很简单&#xff0c;但是…