二分图匹配(一)

文章目录

  • 什么是二分图:
  • 例题:
    • 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,一经查实,立即删除!

相关文章

P3980-[NOI2008]志愿者招募【费用流】

正题 题目链接:https://www.luogu.com.cn/problem/P3980 题目大意 nnn天&#xff0c;第iii天需要AiA_iAi​个志愿者。有mmm种志愿者&#xff0c;第iii种从sis_isi​天服务到tit_iti​天&#xff0c;需要cic_ici​元的费用。 求满足条件的最小费用 1≤n≤1000,1≤m≤100001\le…

不止代码:恐狼后卫(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…

P4180-[BJWC2010]严格次小生成树【Kruskal,倍增】

正题 题目链接:https://www.luogu.com.cn/problem/P4180 题目大意 nnn个点mmm条边的一张无向图&#xff0c;求它的严格次小生成树。 1≤n≤105,1≤m≤31051\leq n\leq 10^5,1\leq m\leq 3\times 10^51≤n≤105,1≤m≤3105 解题思路 一定存在一种严格次小生成树和最小生成树只…

【.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…

[BZOJ5312]冒险(势能线段树)

[BZOJ5312]冒险 维护一个长度为 n 的序列&#xff0c;支持 m 次操作&#xff0c;操作包括区间按位或一个数&#xff0c;区间按位与一个数&#xff0c;以及查询区间最大值。 线段树每个节点上维护区间与、区间或和区间最大值。 如果一次操作对区间与的影响和对区间或的影响相同&…

牛客题霸 [表达式求值] C++题解/答案

牛客题霸 [表达式求值] C题解/答案 题目描述 请写一个整数计算器&#xff0c;支持加减乘三种运算和括号。 题解&#xff1a; 没有除法emmm 我们从头开始依次判断每个字符 如果是左括号&#xff0c;我们就找右括号&#xff0c;并截取括号内的数字 记录上一次的符号&#xff…

P6085-[JSOI2013]吃货JYY【状压dp,欧拉回路】

正题 题目链接:https://www.luogu.com.cn/problem/P6085 题目大意 nnn个点的一张无向图&#xff0c;有kkk条必走边&#xff0c;mmm条其他边&#xff0c;求从111出发经过必走边后回到起点的最短路径。 2≤n≤13,0≤k≤78,2≤m≤2002\leq n\leq 13,0\leq k\leq 78,2\leq m\leq 2…

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

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

不止代码:合唱队列(动态规划)

文章目录题目描述解析1.n^2^朴素算法2.队列nlogn算法代码3.二维DP&#xff08;n^2^)代码thanks for reading!题目描述 五一到了&#xff0c;PKU-ACM队组织大家去登山观光&#xff0c;队员们发现山上一个有N个景点&#xff0c;并且决定按照顺序来浏览这些景点&#xff0c;即每次…

#6029. 「雅礼集训 2017 Day1」市场(势能,区间除)

#6029. 「雅礼集训 2017 Day1」市场 用线段树维护数列&#xff0c;区间上维护最大最小值&#xff0c;区间和还有标记&#xff0c;修改时&#xff0c;区间加直接做&#xff0c;而区间除时&#xff0c;递归到线段树上某一区间&#xff0c;如果这一操作等价于区间加&#xff08;也…

牛客题霸 [容器盛水问题] C++题解/答案

牛客题霸 [容器盛水问题] C题解/答案 题目描述 给定一个整形数组arr&#xff0c;已知其中所有的值都是非负的&#xff0c;将这个数组看作一个容器&#xff0c;请返回容器能装多少水。 具体请参考样例解释 题解&#xff1a; 我们找出容器的左右边界&#xff0c;选择边界更低…

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

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

Loj#3026-「ROIR 2018 Day1」管道监控【Trie,费用流】

正题 题目链接:https://loj.ac/p/3026 题目大意 给出nnn个点的一棵外向树&#xff0c;然后mmm个字符串和费用表示你每次可以花费这个费用覆盖路径字符串和给出字符串相等的路径&#xff0c;求覆盖所有边的最小花费&#xff08;可以重复覆盖&#xff09; 输出方案 1≤n≤500,…

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

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

P3224 [HNOI2012]永无乡(并查集+权值线段树合并/平衡树)

[HNOI2012]永无乡 Code1 权值线段树天然支持merge&#xff0c;线段树上二分求第k小 #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;chgetcha…

牛客题霸 [没有重复项数字的所有排列] C++题解/答案

牛客题霸 [没有重复项数字的所有排列] C题解/答案 题目描述 给出一组数字&#xff0c;返回该组数字的所有排列 例如&#xff1a; [1,2,3]的所有排列如下 [1,2,3],[1,3,2],[2,1,3],[2,3,1],[3,1,2], [3,2,1]. &#xff08;以数字在数组中的位置靠前为优先级&#xff0c;按字典…

.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;机器要尽可能的用完 这河里吗 可能是机器多了…