目录
(4)例题:
看了往上诸多博客,感觉讲的稍有欠缺,于是自己写了一篇(耐心浏览)
最小点覆盖 :
(1) 定义:
假如选了一个点就相当于覆盖了以它为端点的所有边。最小顶点覆盖就是选择最少的点来覆盖所有的边。
(2)结论:
最小点覆盖 = 最大匹配数
(3) 证明:
先了解一些概念
匹配点:匹配边两端的端点
增广路径,即:一边的非匹配点到另一边的非匹配点的一条非匹配边和匹配边交替经过的路径. (开头和结尾一定是非匹配点且不在一个集合中)
如图绿色的线就是增广路径:
然后简单说一下匈牙利算法:
(1)找出一条增广路,通过将增广路取反,得到新的M’来代替原M,新的匹配数相对于原匹配数多1
(2)重复(1)操作,直至匹配过程中找不到增广路为止。
进入正题:
证 最小点覆盖 == 最大匹配数
设 最小点覆盖数 为 s ,最大匹配数 为 m
(1)首先 证 s >= m
我们知道 当我们得到最大匹配数 时,所有边是没有交点的,即所有匹配边没有公共点,这个时候我们就可以在这些边中选一端作为覆盖点 ,我们要找到最小个数的点集使其关联所有边,必然包含我们从匹配边中的到的覆盖点,所以 证得 s>=m
(2) 接下来我们证明 m为最小覆盖点数 (通过构造方法)
我们先通过构造得到一个点集,再证明这些点就是最小覆盖点
构造方法如下:
1.我们在一个二分图中先得到最大匹配
如图红色的是匹配边
2.从左边的每一个非匹配点开始做增广路径(即使不能构成增广路径,因为成功的话就不是最大匹配了,最大匹配后没有增广路径),并且对走过的点进行标记
3.选择方式:选出左边所有未被标记的点+右边所有被标记的点(用红色圈出),我们设为 K
我们通过上面的描述可以得到三个性质:
- 右部所有非匹配点一定没有被标记(因为右边非匹配点被标记的话,就会形成增广路径,这会和我们在最开始已经进行最大匹配矛盾)
- 左部非匹配点一定被标记,因为起点增广标记起点就是非匹配点
- 一条匹配边要么同时被标记要么同时不被标记
接下来 我们证明 k == m
所有边可以分为匹配边和非匹配边,所以我们可以分开说明
先证明 这些点覆盖了所有匹配边(所有点都是匹配点)
- 在选择的时候,我们选择左边所有未被标记的点,即匹配点(性质二可知),形如图中的最下面那条匹配边的左边那个点
- 我们还要选择 右边被标记的点,即由于性质一可以得到
- 对于每个匹配边,左右两个点要么同时被标记,要么同时不被标记,
- 同时被标记的点所在的匹配边,由于我们选择了右边所有被标记的点,(因为增广路径,我们从非匹配点出发,左边被标记一定是因为右边的点先被标记了,所以他是根据右边的标记点来的) ,所以这些匹配边我们全选了
- 同时未标记的点所在的匹配边,由于我们选择了左边所有不被标记的点,所以这些匹配边我们全选了
可知,我们的选择是所有的m条匹配边,且每个匹配边我们只会选择左,或者右边一个点
一共m个点
由此我们证明了k个点都是匹配点 , 覆盖了所有匹配边
再者证明 这些点覆盖了所有非匹配边
而非匹配边又可以分为两种情况:
- (1)左边是非匹配点 ------ 右边是匹配点 (形如图中 s1 的边) :
- 每次我们选左边非匹配点作为起点做增广路径,所以右边的点一定被标记,且在选择的时候,我们会选择右边的标记的匹配点, 所以选出来的点在此情况下覆盖了非匹配边
- (2)左边是匹配点------右边是匹配点
- 形如 s2 ,如果存在这样的边,就说明存在增广路径,和我们最开始求了最大匹配矛盾了,所以在最大匹配后不存在这种情况
所以由此证明选出来的点包含所有边,且 个数 等于从匹配边里选出的覆盖点,即最大匹配数(每条边选一个点)
所以,
最小点覆盖 = 最大匹配数 !!!!!!!!!!!
证毕!
心累,点赞吧!!!!
例题:
输入:
5 5 10 0 1 1 1 1 2 2 1 3 3 1 4 4 2 1 5 2 2 6 2 3 7 2 4 8 3 3 9 4 3 0
输出:
3
上代码:
#include<iostream>
#include<cstring>
using namespace std;
const int maxn = 1e3;
int g[maxn][maxn];
int match[maxn];
int vis[maxn];
int n,m;
int find(int x){
for(int i = 1;i<=m;i++){
if(vis[i]||!g[x][i])continue;
vis[i] = 1;
if(match[i] == 0 || find(match[i]) ){
match[i] = x;
return 1;
}
}
return 0;
}
int main(){
int k=0;
while(cin>>n,n){
memset(match,0,sizeof match);
memset(g,0,sizeof g);
cin>>m>>k;
while(k--){
int a,b,c;
cin>>a>>b>>c;
if(!b||!c)continue;
g[b][c] = 1;
}
int res = 0 ;
for(int i = 1;i<=n;i++){
memset(vis,0,sizeof vis);
if(find(i))res++;
}
cout<<res<<endl;
}
}