先看如下的两个定义:

  • 点-双连通图:一个连通的无向图内部没有割点,那么该图是点-双连通图。

        注意:孤立点,以及两点一边这两种图都是点-双连通的,因为它们都是内部无割点。

  • 边-双连通图:一个连通的无向图内部没有桥,那么该图就是边-双连通图。

        注意:孤立点是边-双连通的,但是两点一边不是边-双连通的。

由上面定义可以知道:点-双连通图不一定是边-双连通的

对于一张无向图,点-双连通的极大子图称为双连通分量。不难发现,每条边恰好属于一个双连通分量(所以两点一边是一个点-双连通分量)。但不同双连通分量可能会有公共点,可以证明不同双连通分量最多只有一个公共点,且它一定是割顶。另一方面任意割顶都是至少两个不同的点-双连通分量的公共点。

边-双连通的极大子图称为边-双连通分量。除了桥不属于任何边-双连通分量外,其他每条边恰好属于一个边-双连通分量,而且把所有桥删除之后,每个连通分量对应原图中的一个边-双连通分量。

总之:

判断一个图是不是点-双连通的只要看图中是否有割点。(比寻找割点多了个栈而已)

判断一个图是不是边-双连通的只要看图中是否有桥。

如果看到这里还有很多困惑,先画图思考矛盾的地方,然后看代码可以更清楚,前提是能看懂求割点的算法。

 #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)的更多相关文章

  1. [置顶] 刘汝佳《训练指南》动态规划::Beginner (25题)解题报告汇总

    本文出自   http://blog.csdn.net/shuangde800 刘汝佳<算法竞赛入门经典-训练指南>的动态规划部分的习题Beginner  打开 这个专题一共有25题,刷完 ...

  2. Windows XP硬盘安装Ubuntu 12.04双系统图文详解

    Windows XP硬盘安装Ubuntu 12.04双系统图文详解 Ubuntu 12.04 LTS版本于2012年4月26日发布,趁着五一放假,赶紧在自己的Windows XP的电脑上安装下Ubun ...

  3. Windows XP硬盘安装Ubuntu 16.04双系统图文详解

    需要下载的东西有两个,一个是grub4dos,另一个是Ubuntu 16.04 LTS的镜像文件,具体下载地址如下:   1 2 3       1.grub4dos  点击下载 grub4dos 2 ...

  4. c++20701除法(刘汝佳1、2册第七章,暴搜解决)

    20701除法 难度级别: B: 编程语言:不限:运行时间限制:1000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述     输入正整数n,按从小到大的顺序输出所有 ...

  5. 刘汝佳黑书 pku等oj题目

    原文地址:刘汝佳黑书 pku等oj题目[转]作者:小博博Mr 一.动态规划参考资料:刘汝佳<算法艺术与信息学竞赛><算法导论> 推荐题目:http://acm.pku.edu. ...

  6. 刘汝佳 算法竞赛-入门经典 第二部分 算法篇 第五章 1(String)

    第一题:401 - Palindromes UVA : http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8 ...

  7. ACM题目推荐(刘汝佳书上出现的一些题目)[非原创]

    原地址:http://blog.csdn.net/hncqp/article/details/1758337 推荐一些题目,希望对参与ICPC竞赛的同学有所帮助. POJ上一些题目在http://16 ...

  8. 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,再 ...

  9. HDU 1003 Max Sum【动态规划求最大子序列和详解 】

    Max Sum Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Sub ...

随机推荐

  1. Unix守护进程

    问题描述:         Unix守护进程 问题解决:     Unix守护进程没有控制终端,终端名设置为问号(?),终端前台进程组ID设置(TPGID)为-1 守护进程编写规则:      (1) ...

  2. android 解析XML方式(三)

    上一节中,我们使用SAX方式解析xml文档, SAX方式是基于事件驱动的.当然android的事件机制是基于回调函数的.在这一节中,我们用另外一种方式解析xml文档,这种方式也是基于事件驱动的,与SA ...

  3. HDU4758 Walk Through Squares AC自动机&&dp

    这道题当时做的时候觉得是数论题,包含两个01串什么的,但是算重复的时候又很蛋疼,赛后听说是字符串,然后就觉得很有可能.昨天队友问到这一题,在学了AC自动机之后就觉得简单了许多.那个时候不懂AC自动机, ...

  4. [STL]单词转换

    如果单词转换文件的内容是: 'em         themcuz         becausegratz      grateful i             Inah        nopos ...

  5. [shell编程]一个简单的脚本

    首先,为什么要学习shell呢?哈哈,当然不是shell能够怎样怎样然后100字. 最近看到一篇博文<开阔自己的视野,勇敢的接触新知识>,读完反思良久.常常感慨自己所会不多,对新知识又有畏 ...

  6. 空格的URL编码

    Q: 为什么我看的教材一会说是“+” 一会说是“%20” A: urlencode(" ") '返回+encodeURI(" ") '返回%20是有区别的

  7. FileOutputStream和FileInputStream

    package one.string; import java.io.File; import java.io.FileInputStream; import java.io.FileNotFound ...

  8. Android中自定义Checkbox

    custom_checkbox.xml文件: <?xml version="1.0" encoding="utf-8"?> <selector ...

  9. MySQL 建表字段长度的限制

    脑补,varchar(N),N指的是最大字符数,不是字节数. 先上测试说明:   在MySQL建表时,遇到一个奇怪的现象: root@localhost : test 10:30:54>CREA ...

  10. Android:控件AutoCompleteTextView 客户端保存搜索历史自动提示

    <?xml version="1.0" encoding="utf-8"?> <LinearLayout xmlns:android=&quo ...