【学习】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 具有脚本语言中最丰富和强大的类库,足以支持绝大多数日常应用 ...
随机推荐
- 在caffe-ssd安装编译环境运行make all时候报错:Makefile:572: recipe for target '.build_release/src/caffe/util/hdf5.o' failed make: *** [.build_release/src/caffe/util/hdf5.o] Error 1
解决办法: 修改:Makefile.config INCLUDE_DIRS /usr/include/hdf5/serial/ 修改:Makefile LIBRARIES hdf5_hl and hd ...
- python中的lxml模块
Python中自带了XML的模块,但是性能不太好,相比之下,LXML增加了很多实用的功能. lxml中主要有两部分, 1) etree,主要可以用来解析XML字符串, 内部有两个对象,etree._E ...
- yii2 rules验证规则,ajax验证手机号码是否唯一
<?php namespace frontend\models; use Yii; use yii\base\Model; /** * Signup form */ class SignupFo ...
- pyspider 示例
数据存放目录: C:\Users\Administrator\data 升级版(可加载文章内所有多层嵌套的图片标签) #!/usr/bin/env python # -*- encoding: utf ...
- 转:判断Caps Lock键是否打开,如果打开则关闭
using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; usin ...
- 蓝牙Profile的概念和常见种类(转)
蓝牙Profile Bluetooth的一个很重要特性,就是所有的Bluetooth产品都无须实现全部 的Bluetooth规范.为了更容易的保持Bluetooth设备之间的兼容,Bluetooth规 ...
- AtCoder Beginner Contest 044 A - 高橋君とホテルイージー / Tak and Hotels (ABC Edit)
Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement There is a hotel with ...
- Web开发笔记 #07# Swagger Editor
Swagger Editor是一款可以用yaml格式进行RESTful API设计.可视化.测试的工具,并且能够实时看到自动生成的文档.效果大概是这样的↓ 根据官方网站介绍,如果是团队的话,建议用在线 ...
- hdu some problems in Multi-University Training Contest
hdu 6103 Kirinriki #include<bits/stdc++.h> using namespace std; int n,m,ans; ]; void doit(int ...
- selenium webdriver 实现Canvas画布自动化测试
https://blog.csdn.net/xiaoguanyusb/article/details/80324210 由借鉴意义, 转过来 canvas 是一个画布,定位元素时只能定位到画布上,如下 ...