求双连通分量的详解。(根据刘汝佳的训练指南p314)
先看如下的两个定义:
- 点-双连通图:一个连通的无向图内部没有割点,那么该图是点-双连通图。
注意:孤立点,以及两点一边这两种图都是点-双连通的,因为它们都是内部无割点。
- 边-双连通图:一个连通的无向图内部没有桥,那么该图就是边-双连通图。
注意:孤立点是边-双连通的,但是两点一边不是边-双连通的。
由上面定义可以知道:点-双连通图不一定是边-双连通的。
对于一张无向图,点-双连通的极大子图称为双连通分量。不难发现,每条边恰好属于一个双连通分量(所以两点一边是一个点-双连通分量)。但不同双连通分量可能会有公共点,可以证明不同双连通分量最多只有一个公共点,且它一定是割顶。另一方面任意割顶都是至少两个不同的点-双连通分量的公共点。
边-双连通的极大子图称为边-双连通分量。除了桥不属于任何边-双连通分量外,其他每条边恰好属于一个边-双连通分量,而且把所有桥删除之后,每个连通分量对应原图中的一个边-双连通分量。
总之:
判断一个图是不是点-双连通的只要看图中是否有割点。(比寻找割点多了个栈而已)
判断一个图是不是边-双连通的只要看图中是否有桥。
如果看到这里还有很多困惑,先画图思考矛盾的地方,然后看代码可以更清楚,前提是能看懂求割点的算法。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
const int maxn=+; int n,m;
int bcc_cnt;
int dfs_clock; //bcc_cnt计数一共有多少个点-双连通分量
int pre[maxn]; //vis标记,同时也标记了是树中的第几个点
bool iscut[maxn];
int bccno[maxn]; //bccno[i]=x表示第i个顶点属于x号点双连通分量
vector<int> G[maxn],bcc[maxn]; //bcc[i]中包含了i号点-双连通分量的所有节点 struct Edge //边的结构体
{
int u,v;
Edge(int u,int v):u(u),v(v){}
};
stack<Edge> S; int dfs(int u,int fa)
{
int lowu=pre[u]=++dfs_clock;
int child=;
for(int i=;i<G[u].size();i++)
{
int v=G[u][i]; //取出点
Edge e = Edge(u,v); //创建这条边 if(!pre[v]) //v没有被访问过
{
S.push(e); //将边入栈
child++;
int lowv=dfs(v,u); //求low先
lowu=min(lowu,lowv);
if(lowv >= pre[u]) //本节点是割点
{
iscut[u]=true;
bcc_cnt++; //注意bcc_cnt从1开始编号
bcc[bcc_cnt].clear(); //清除之前留下的
while(true) //产生一个双连通分量,
{
Edge x=S.top(); //逐次取出边
S.pop();
//1个点可能属于多个连通分量,且它一定是割点。
if(bccno[x.u]!=bcc_cnt) //这个点还没有统计到这个连通分量。
{
bcc[bcc_cnt].push_back(x.u);
bccno[x.u]=bcc_cnt; //预防重复统计
}
if(bccno[x.v]!=bcc_cnt)
{
bcc[bcc_cnt].push_back(x.v);
bccno[x.v]=bcc_cnt;
}
if(x.u==u && x.v==v) //扫到u-v,栈中又没有与u相连的边了。继续试试其他孩子
break;
}
}
}
else if(pre[v]<pre[u]&&v!=fa) //点v在u上面就被访问过,才可以更新,在下面访问过的,不可以!
{
S.push(e); //这个是和u在一起的双连通分量
lowu=min(lowu,pre[v]);
}
} /*
根的孩子必须大于1才会是割点,有割点才会有双连通分量。
(1)那么如果根不是割点呢?
假设根不是割点,那么根最多只有1个孩子,也就是说根的度为1,那么根不可能处于任何1个双连通分量中。
假设根是割点,那么每个孩子各自是一个连通分量。那么就会在上面的代码中被处理为一个双联通分量。
(2)如果有桥呢?比如u-v是桥,那么会怎样?
假设u-v是桥,且u在数中的时间戳比较小。可知v也就是一个割点啦,u-v断开后,与v相连的都成为一个双连通分量了。
回溯到u时,栈中(或顶)没有包含u的边,直到另一个连通分量的产生。
如果u的孩子中没有连通分量了,那么与u相连的孩子肯定有边连到u的上边,他们又形成了一个环了,双连通分量又产生了,由其他割点去解决。
*/
if(fa< && child==) iscut[u]=false;
return lowu;
} void find_bcc(int n)
{
memset(pre,,sizeof(pre));
memset(iscut,,sizeof(iscut));
memset(bccno,,sizeof(bccno));
dfs_clock = bcc_cnt = ;
for(int i=;i<n;i++) //为了防止有多个连通图,全部都得搜
if(!pre[i]) dfs(i,-);
}
int main()
{
while(scanf("%d%d",&n,&m)==&&n)
{
for(int i=;i<n;i++) G[i].clear(); //点集
for(int i=;i<m;i++) //输入边
{
int u,v;
scanf("%d%d",&u,&v);
G[u].push_back(v);
G[v].push_back(u);
}
find_bcc(n); //计算双连通分量的个数
printf("点-双连通分量一共%d个\n",bcc_cnt); for(int i=;i<=bcc_cnt;i++) //输出每个双连通分量。可能点A在第一个双连通分量中输出,又出现在第2个双连通分量中,因为它是割点。
{
printf("第%d个点-双连通分量包含以下点:\n",i);
sort(&bcc[i][],&bcc[i][]+bcc[i].size()); //对vector排序,使输出的点从小到大
for(int j=;j<bcc[i].size();j++)
{
printf("%d ",bcc[i][j]);
}
printf("\n");
}
}
return ;
} 带注释的源码
带注释的源码
需要注意的是:
- 单个点不算是一个双连通分量,比如仅有1个点的图。
- 两个点一条边的子图算是一个双连通分量。比如:a-b-c 这个图中有两个双连通分量:a-b和b-c。
- 算法中已经考虑到根和叶子和非叶子节点,所以不用例外去添加代码。
- 根据需要可删减。比如可以不用vector来存储所有连通分量等等。
求双连通分量的详解。(根据刘汝佳的训练指南p314)的更多相关文章
- [置顶] 刘汝佳《训练指南》动态规划::Beginner (25题)解题报告汇总
本文出自 http://blog.csdn.net/shuangde800 刘汝佳<算法竞赛入门经典-训练指南>的动态规划部分的习题Beginner 打开 这个专题一共有25题,刷完 ...
- Windows XP硬盘安装Ubuntu 12.04双系统图文详解
Windows XP硬盘安装Ubuntu 12.04双系统图文详解 Ubuntu 12.04 LTS版本于2012年4月26日发布,趁着五一放假,赶紧在自己的Windows XP的电脑上安装下Ubun ...
- Windows XP硬盘安装Ubuntu 16.04双系统图文详解
需要下载的东西有两个,一个是grub4dos,另一个是Ubuntu 16.04 LTS的镜像文件,具体下载地址如下: 1 2 3 1.grub4dos 点击下载 grub4dos 2 ...
- c++20701除法(刘汝佳1、2册第七章,暴搜解决)
20701除法 难度级别: B: 编程语言:不限:运行时间限制:1000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述 输入正整数n,按从小到大的顺序输出所有 ...
- 刘汝佳黑书 pku等oj题目
原文地址:刘汝佳黑书 pku等oj题目[转]作者:小博博Mr 一.动态规划参考资料:刘汝佳<算法艺术与信息学竞赛><算法导论> 推荐题目:http://acm.pku.edu. ...
- 刘汝佳 算法竞赛-入门经典 第二部分 算法篇 第五章 1(String)
第一题:401 - Palindromes UVA : http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8 ...
- ACM题目推荐(刘汝佳书上出现的一些题目)[非原创]
原地址:http://blog.csdn.net/hncqp/article/details/1758337 推荐一些题目,希望对参与ICPC竞赛的同学有所帮助. POJ上一些题目在http://16 ...
- win7 64 + Ubuntu 14.04.1 64双系统安装,详解UEFI ~ GPT和legacy ~ MBR区别
win7 64 + Ubuntu 14.04.1 64双系统安装 背景:我的笔记本之前的系统是window 7 64 + Ubuntu 14.04.1,用UEFI引导系统.安装过程是先装的win7,再 ...
- HDU 1003 Max Sum【动态规划求最大子序列和详解 】
Max Sum Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Sub ...
随机推荐
- OI 回忆录
时间过得好快,一下子就高三了,一下子就退役了,两年的时间仿佛就在一眨眼的功夫内度过了.不过还是想回忆回忆这两年的雨雪风霜,也就当做个总结吧. 高一其实并没有什么好说的,反正就这么颓到了高二. 高二上学 ...
- crud springmvc
实体类:Student.java package demo.entity; public class Student { private int id; private String name; pr ...
- Unity3D研究院之IOS全自动编辑framework、plist、oc代码
Unity打IOS时会先生成一个Xcode工程,如果你需要增加一些第三方的framework那么需要手动一条一条的添加,这太烦了..而且可能你还需要修改Plist文件,甚至还可能要修改unity自动生 ...
- linux故障救援
今天手贱,误操作将/usr/lib整个文件夹都删了.但是系统内还有很多重要的文件和数据,不敢贸然重装系统. 弄了个u盘启动. 一开始整个系统只有一个分区,原系统装在这个分区里,我用u盘启动后,用分区管 ...
- 360 chrome不能登录Google账户
用过Google chrome,枫树,360 chrome国际版,今天还是回到Google chrome了. 因为除了Chrome自己,其他的都不能正常登录Google账户(偶尔偶尔一次可以,直到今天 ...
- *[codility]Peaks
https://codility.com/demo/take-sample-test/peaks http://blog.csdn.net/caopengcs/article/details/1749 ...
- jdbc的通讯录CRUD
基于JDBC的通讯录练手:项目以MVC模式开发,包名:cn.itcast.txl.domain;cn.itcast.txl.dao;cn.itcast.txl.dao.impl;cn.itcast.t ...
- LR_问题_无法使用LR的Controller,提示缺少license
问题描述 无法使用LR的Controller,提示缺少license 问题解决 使用开始->所有程序->HP LoadRunner->loadrunner,在打开界面的左上角选择co ...
- PCB工艺镀金(电金)和沉金(化金)的区别
1.镀金和沉金的别名分别是什么? 镀金:硬金,电金(镀金也就是电金) 沉金:软金,化金 (沉金也就是化金) 2.别名的由来: 镀金:通过电镀的方式,使金粒子附着到pcb板上,所以叫电金,因为附着 ...
- java:接口实例
接口:打印机接口 interface Printer { public void read(); } 函数一:佳能打印机 class CanPrinter implements Printer { p ...