【学习】Hall’s Marriage Theorem
其实是在做题时遇到这个定理的。
这个定理的图论意义是:
对于一个二分图\(G=\{X+Y,E\}\),它满足:
\(\forall W \subseteq X, \, |W| \leq |N_G(W)|\) \(\iff\)\(X\)中的每个结点都有匹配.
其中\(N_G(W)\)为图\(G\)中所有与\(W\)相连的结点的集合。
必要性显然。
对于充分性,不会……以后再补充。
由这个定理,我们能得到一个推论:
二分图\(G\)的最大匹配\(M\)等于\(|X| - \max (|W| - |N_G(W)|)\).
顺便一提,我们允许\(W = \emptyset\),故\(\max (|W| - |N_G(W)|) \geq 0\)。
我们分两步来证明这个推论。在此,我们假设\(W_0 \subseteq X\)满足\(|W_0| - |N_G(W_0)| = \max (|W| - |N_G(W)|)\)
- 证明:\(M \leq |X| - \max (|W| - |N_G(W)|)\).
考虑\(W_0\),即使\(N_G(W_0)\)中的每个结点都与\(W_0\)中的结点匹配,也会有\(|W_0| - |N_G(W_0)|\)个结点得不到匹配。故\(X\)中至少有\(|W_0| - |N_G(W_0)|\)个结点得不到匹配。因此,\(M \leq |X| - \max (|W|)\)。 - 证明:\(M \geq |X| - \max (|W| - |N_G(W)|)\).
我们在\(Y\)中添加\(\max (|W| - |N_G(W)|\)个结点与\(X\)中所有结点相连得到新图\(G`\),那么,\(|W_0| - |N_{G`}(W_0)| = 0\)。由Hall‘s Marriage Theorem可知,图\(G`\)中\(X\)的每个元素都有匹配。而对于\(W_0\),因为\(|W_0| = |N_{G`}(W_0)|\),所以\(N_{G`}(W_0)\)中的每个元素都与\(W_0\)中的元素匹配,可得我们新加入的结点都有匹配。因此,除\(W_0\)中有\(|W_0|-|N_G(W_0)|\)个结点与新加入的结点匹配之外,\(X\)中其余结点都与原来\(Y\)中的结点匹配。故在删去新加的结点后,\(X\)中至多有\(|W_0| - |N_G(W_0)|\)个结点没有匹配。则\(M \geq |X| - \max (|W| - |N_G(W)|)\)。
那么,这个定理及其推论有什么用呢?因为要枚举\(X\)的所有子集,所以一般是没什么用的。然而,某些有特殊性质的二分图的最大匹配问题,会使用这个定理及其推论来转化问题。
例题
arc076F Exhausted?
题意:给出二分图\(G=\{X+Y,E\}\),\(X\)中的所有结点满足:若其编号为\(i\),则只与\(Y\)中编号小于等于\(L_i\)或编号大于等于\(R_i\)的结点连边。给出\(|X|,|Y|\)和所有的\(L_i,R_i\),求\(G\)的最大匹配\(M\)。(输出\(|X|-M\))
\(|X|,|Y| \leq 2 \times 10^5\).
有了这个定理,问题就简单了。我们只要求\(\max (|W| - |N_G(W)|)\)。对于\(W \subseteq X\),\(N_G(W)\)就是形如\(Y - [l,r]\)的形式,其中\([l,r]\)是\(Y\)中一个编号的区间。我们枚举\(r\),用线段树实现区间加和区间查最大值即可。具体做法不再阐述。
时间复杂度\(O(n \log n)\)。
#include <bits/stdc++.h>
using namespace std;
const int N = 200010;
int n,m,ans,cnt;
struct data {
int l,r;
} ordl[N],ordr[N];
bool cmpl(data a,data b) {
return a.l < b.l;
}
bool cmpr(data a,data b) {
return a.r < b.r;
}
struct node {
int mx,tag;
} t[N << 2];
void puttag(int x,int val) {
t[x].mx += val;
t[x].tag += val;
}
void push_down(int x) {
puttag(x<<1,t[x].tag);
puttag(x<<1|1,t[x].tag);
t[x].tag = 0;
}
void push_up(int x) {
t[x].mx = max(t[x<<1].mx,t[x<<1|1].mx);
}
void modify(int l,int r,int val,int x=1,int lp=1,int rp=m) {
if (l > rp || r < lp) return;
if (lp >= l && rp <= r) {
puttag(x,val);
return;
}
if (t[x].tag) push_down(x);
int mid = (lp + rp) >> 1;
modify(l,r,val,x<<1,lp,mid);
modify(l,r,val,x<<1|1,mid+1,rp);
push_up(x);
}
int query(int l,int r,int x=1,int lp=1,int rp=m) {
if (l > rp || r < lp) return 0;
if (lp >= l && rp <= r) return t[x].mx;
if (t[x].tag) push_down(x);
int mid = (lp + rp) >> 1;
return max(query(l,r,x<<1,lp,mid),query(l,r,x<<1|1,mid+1,rp));
}
int main() {
int a,b;
scanf("%d%d",&n,&m);
ans = max(0,n - m);
cnt = n;
for (int i = 1 ; i <= cnt ; ++ i) {
scanf("%d%d",&a,&b);
a ++;
b --;
if (a > b) {
i --, cnt --;
continue;
}
ordl[i] = (data) {a,b};
ordr[i] = (data) {a,b};
}
sort(ordl+1,ordl+cnt+1,cmpl);
sort(ordr+1,ordr+cnt+1,cmpr);
int pl = 1, pr = 1;
for (int i = 1 ; i <= m ; ++ i) {
modify(1,i,1);
for (; ordl[pl].l <= i && pl <= cnt ; ++ pl)
modify(ordl[pl].l,m,1);
ans = max(ans,query(1,i) - m);
for (; ordr[pr].r <= i && pr <= cnt ; ++ pr)
modify(ordr[pr].l,m,-1);
}
printf("%d\n",ans);
return 0;
}
小结:学了一个没卵用的定理。也挺好奇还有哪些二分图能用到这个定理。
upd 2018.12.26
假如\(G\)满足,\(X\)中的每一个结点与\(Y\)中一个区间内的结点连边,即\(\forall a \in X, N_G(a) = [l_a,r_a]\)。
定义对于\(Y\)中的一个编号区间\([l,r]\),\(f([l,r])\)为\(X\)中只与\([l,r]\)中点连边的点构成的集合,即\(\{ a \in X \ | \ N_G(a) \subseteq [l,r] \}\)。
求证:
\]
证明:
- 必要性。令\(W = f([l,r])\),那显然有\(N_G(W) \subseteq [l,r]\)。因此,\(|f([l,r])| = |W| \leq |N_G(W)| \leq |[l,r]|\)。
- 充分性。考虑任意一个\(W \subseteq X\),那么\(N_G(W)= \bigcup_{a \in W} [l_a,r_a]\)。那么,\(N_G(W)\)就一定可以被表示为若干个互不相交的区间的并。因为每个\(X\)中的结点都只和一个区间内的点连边,所以我们可以把\(W\)划分为若干个集合,其中每个集合对应\(N_G(W)\)中的一段区间。对于我们划分出来的任意一个集合,假如它对应的区间为\([l,r]\),那这个集合就一定是\(f([l,r])\)的子集。因此,对于\(W\)划分出来的每一集合\(S\)都满足\(|S| \leq |N_G(S)|\),故\(W\)也满足\(|W| \leq |N_G(W)|\)。
【学习】Hall’s Marriage Theorem的更多相关文章
- 【CF981F】Round Marriage(二分答案,二分图匹配,Hall定理)
[CF981F]Round Marriage(二分答案,二分图匹配,Hall定理) 题面 CF 洛谷 题解 很明显需要二分. 二分之后考虑如果判定是否存在完备匹配,考虑\(Hall\)定理. 那么如果 ...
- 加州大学伯克利分校Stat2.2x Probability 概率初步学习笔记: Section 4 The Central Limit Theorem
Stat2.2x Probability(概率)课程由加州大学伯克利分校(University of California, Berkeley)于2014年在edX平台讲授. PDF笔记下载(Acad ...
- 关于Hall定理的学习
基本定义 \(Hall\) 定理是二分图匹配的相关定理 用于判断二分图是否存在完美匹配 存在完美匹配的二分图即满足最大匹配数为 \(min(|X|,|Y|)\) 的二分图,也就是至少有一边的点全部被匹 ...
- Matrix_tree Theorem 矩阵树定理学习笔记
Matrix_tree Theorem: 给定一个无向图, 定义矩阵A A[i][j] = - (<i, j>之间的边数) A[i][i] = 点i的度数 其生成树的个数等于 A的任意n ...
- 【CF981F】Round Marriage(二分答案,hall定理)
传送门 题意: 给出一个长度为\(L\)的环,标号从\(0\)到\(L-1\). 之后给出\(n\)个新郎,\(n\)个新娘离起点的距离. 现在新郎.新娘要一一配对,但显然每一对新人的产生都会走一定的 ...
- 【Matrix-tree Theorem学习笔记】
定义度数矩阵\(D(G)\): 定义邻接矩阵\(C(G)\): 定义\(Laplace\)矩阵\(A\) \( A(G) = D(G) - C(G) \) 记图\(G\)的所有生成树权值和为\(t(G ...
- <老友记>学习笔记
这是六个人的故事,从不服输而又有强烈控制欲的monica,未经世事的千金大小姐rachel,正直又专情的ross,幽默风趣的chandle,古怪迷人的phoebe,花心天真的joey——六个好友之间的 ...
- 增强学习(三)----- MDP的动态规划解法
上一篇我们已经说到了,增强学习的目的就是求解马尔可夫决策过程(MDP)的最优策略,使其在任意初始状态下,都能获得最大的Vπ值.(本文不考虑非马尔可夫环境和不完全可观测马尔可夫决策过程(POMDP)中的 ...
- [Python] 学习资料汇总
Python是一种面向对象的解释性的计算机程序设计语言,也是一种功能强大且完善的通用型语言,已经有十多年的发展历史,成熟且稳定.Python 具有脚本语言中最丰富和强大的类库,足以支持绝大多数日常应用 ...
随机推荐
- 关于python中selector问题
在做大型的爬虫时,re表达式往往效率不高,scrapy框架为爬虫提供了很好的爬虫方法 scrapy提取数据时有一套自己的机制,即selectors,一般通过特定的XPath,或者特定 的CSS表达式来 ...
- centos下mysql 5源码安装全过程记录
参考:http://blog.csdn.net/mycwq/article/details/24488691 安装cmake,mysql 5.5以后的版本要通过cmake进行编译 在新装的CentOS ...
- Gambler Bo (高斯消元求特解)
对于图中的每一个点假设点击Xi * m + j 然后每个点都有那么对于每一个点可以列举出一个方程式,n*m个点解n*m个未知数.利用高斯消元就可以解决. 问题就在这个题目可能不止有一个特,所以我们需要 ...
- WEB应用程序:AJAX全套
概述 对于WEB应用程序:用户浏览器发送请求,服务器接收并处理请求,然后返回结果,往往返回就是字符串(HTML),浏览器将字符串(HTML)渲染并显示浏览器上. 1.传统的Web应用 一个简单操作需要 ...
- Spark学习之路 (七)Spark 运行流程
一.Spark中的基本概念 (1)Application:表示你的应用程序 (2)Driver:表示main()函数,创建SparkContext.由SparkContext负责与ClusterMan ...
- mysql安装使用
linux系统 mysql-5.7.14-linux.zip部署包支持在CentOS 6.x/7.x 服务器硬盘大小要求 a) /data/mysql_data 如果存在该独立分区,要求该分区 &g ...
- android安全检测工具,梆梆安全 - 防止反编译|APP安全加固|应用加固|盗版监测
android安全检测工具,梆梆安全 - 防止反编译|APP安全加固|应用加固|盗版监测https://dev.bangcle.com/ 业内专业的应用加固服务供应商 帮助数十万APP抵御破解风险,早 ...
- activiti 报 next dbid
记录一下吧. 今天将生产环境的几个服务节点改成集群模式,其中包含activiti审批服务节点,其中各个服务几点间数据通信采用MQ(与本文无关). 然后报出如题错误. 究其原因就是,在启动activit ...
- Django 创建项目流程
django 项目创建流程 1 创建项目 cmd django-admin startproject 项目名称 pycharm file -- new project -- Django -- 项目名 ...
- sparkStrming 实时插入 mysql 今天使用echart 实现了简单数据展示 很low 但学习必须加深