首页 > 代码库 > UESTC 899 方老师和农场 --双连通分量的构造
UESTC 899 方老师和农场 --双连通分量的构造
首先将原图中的连通分量缩点,一定可以将原图缩成一棵树的形式,然后统计这棵树的叶子节点个数,答案就是(leaf+1)/2。这里不再证明,可以画个图看一下。
(简单说明一下,首先把两个最近公共祖先最远的两个叶节点之间连接一条边,这样可以把这两个点到祖先的路径上所有点收缩到一起,因为一个形成的环一定是双连通的。然后再找两个最近公共祖先最远的两个叶节点,这样一对一对找完,恰好是(leaf+1)/2次,把所有点收缩到了一起。 --Byvoid)
怎么统计呢?用并查集缩点,可以知道,缩点后留下的边全部是原图的桥,这是我们可以用Tarjan求出原图的所有桥,然后枚举每条桥,桥两端的点度数分别+1,就可以求出每个点(缩点后的点)的度数了,找出度数为1的即为叶子节点。
怎么用Tarjan求桥呢?根据Tarjan算法性质可知,若low[v]>dfn[u],则边(u,v)为桥(v被封死在子树内)
如图,若low[v]>dfn[u],则v被封死在u的子树内,删除点u,或者删除边(u,v),都将使v与u的祖先w不连通。
关于Tarjan求桥可见:http://hi.baidu.com/lydrainbowcat/item/f8a5ac223e092b52c28d591c
代码:
#include <iostream>#include <cstdio>#include <cstring>#include <cmath>#include <algorithm>#include <string>#include <vector>#include <map>#include <set>#include <time.h>#include <queue>#include <cctype>#include <utility>#include <numeric>#include <cstdlib>#include <functional>#include <iomanip>#include <sstream>#define Mod 1000000007#define SMod 10007#define INT 2147483647#define pi acos(-1.0)#define eps 1e-4#define lll __int64#define ll long longusing namespace std;#define N 10007vector<int> G[N];struct Bridge{ int u,v;}bg[2*N];int vis[N],low[N],dfn[N],Time;int fa[N],deg[N];int n,m,cnt;int findset(int x){ if(x != fa[x]) fa[x] = findset(fa[x]); return fa[x];}void Tarjan(int u,int father){ low[u] = dfn[u] = ++Time; vis[u] = 1; for(int i=0;i<G[u].size();i++) { int v = G[u][i]; if(v == father) continue; if(!vis[v]) { Tarjan(v,u); low[u] = min(low[u],low[v]); if(low[v] > dfn[u]) //u->v为桥 bg[cnt].u = u,bg[cnt++].v = v; else //否则,u,v同属一个连通分量,合并 { int fx = findset(u); int fy = findset(v); if(fx != fy) fa[fx] = fy; } } else low[u] = min(low[u],dfn[v]); }}int main(){ int i,j,u,v; while(scanf("%d%d",&n,&m)!=EOF) { for(i=0;i<=n;i++) G[i].clear(); cnt = Time = 0; for(i=0;i<m;i++) { scanf("%d%d",&u,&v); G[u].push_back(v); G[v].push_back(u); } for(i=1;i<=n;i++) fa[i] = i; memset(dfn,0,sizeof(dfn)); memset(low,0,sizeof(low)); memset(vis,0,sizeof(vis)); Tarjan(1,-1); //统计桥端度数 memset(deg,0,sizeof(deg)); for(i=0;i<cnt;i++) { int fx = findset(bg[i].u); //fx,fy为缩点后的代表点 int fy = findset(bg[i].v); deg[fx]++; deg[fy]++; } int leaf = 0; for(i=1;i<=n;i++) if(deg[i] == 1) leaf++; printf("%d\n",(leaf+1)/2); } return 0;}
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。