参考:刘汝佳,训练指南
饶齐博客!
求无向图的割顶与桥代码:http://blog.csdn.net/u013480600/article/category/2335227
//求无向图的割顶和桥
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=100000+10;
int n,m;
int dfs_clock;//时钟,每访问一个节点增1
vector<int> G[maxn];//G[i]表示i节点邻接的所有节点
int pre[maxn];//pre[i]表示i节点被第一次访问到的时间戳,若pre[i]==0表示i还未被访问
int low[maxn];//low[i]表示i节点及其后代能通过反向边连回的最早的祖先的pre值
bool iscut[maxn];//标记i节点是不是一个割点
//求出以u为根节点(u在DFS树中的父节点是fa)的树的所有割顶和桥
//初始调用为dfs(root,-1);
int dfs(int u,int fa)
{
int lowu=pre[u]=++dfs_clock;
int child=0; //子节点数目
for(int i=0; i<G[u].size(); i++)
{
int v=G[u][i];
if(!pre[v])
{
child++;//未访问过的节点才能算是u的孩子
int lowv=dfs(v,u);
lowu=min(lowu,lowv);
if(lowv>=pre[u])
{
iscut[u]=true; //u点是割顶
if(lowv>pre[u]) //(u,v)边是桥
printf("边(%d, %d)是桥\n",u,v);
}
}
else if(pre[v]<pre[u] && v!=fa)//v!=fa确保了(u,v)是从u到v的反向边
{
lowu=min(lowu,pre[v]);
}
}
if(fa<0 && child==1 )
iscut[u]=false;//u若是根且孩子数<=1,那u就不是割顶
return low[u]=lowu;
}
int main()
{
while(scanf("%d%d",&n,&m)==2&&n)
{
dfs_clock=0;//初始化时钟
memset(pre,0,sizeof(pre));
memset(iscut,0,sizeof(iscut));
for(int i=0;i<n;i++) G[i].clear();
for(int i=0;i<m;i++)
{
int u,v;
scanf("%d%d",&u,&v);
G[u].push_back(v);
G[v].push_back(u);
}
dfs(0,-1);//初始调用
for(int i=0;i<n;i++)if(iscut[i]==true)
printf("割顶是:%d\n",i);
}
return 0;
}
求无向图的双连通分量代码:
http://blog.csdn.net/u013480600/article/details/44835827
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
const int maxn=1000+10;
int n,m;
int bcc_cnt;
int dfs_clock;//bcc_cnt计数一共有多少个点-双连通分量
int pre[maxn];
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=0;
for(int i=0;i<G[u].size();i++)
{
int v=G[u][i];
Edge e = Edge(u,v);
if(!pre[v])
{
S.push(e);
child++;
int lowv=dfs(v,u);
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();
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) break;
}
}
}
else if(pre[v]<pre[u] && v!=fa) //这个判断条件如果少了,就是WA,可修改POJ2942代码
{
S.push(e);
lowu=min(lowu,pre[v]);
}
}
if(fa<0 && child==1) iscut[u]=false;
return lowu;
}
void find_bcc(int n)
{
memset(pre,0,sizeof(pre));
memset(iscut,0,sizeof(iscut));
memset(bccno,0,sizeof(bccno));
dfs_clock = bcc_cnt = 0;
for(int i=0;i<n;i++)
if(!pre[i]) dfs(i,-1);
}
int main()
{
while(scanf("%d%d",&n,&m)==2&&n)
{
for(int i=0;i<n;i++) G[i].clear();
for(int i=0;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=1;i<=bcc_cnt;i++)
{
printf("第%d个点-双连通分量包含以下点:\n",i);
sort(&bcc[i][0],&bcc[i][0]+bcc[i].size()); //对vector排序,使输出的点从小到大
for(int j=0;j<bcc[i].size();j++)
{
printf("%d ",bcc[i][j]);
}
printf("\n");
}
}
return 0;
}
/*
示例输入:
6 7
1 2
2 3
1 3
3 4
4 5
3 5
5 6
输出:
点-双连通分量一共3个
第1个点-双连通分量包含以下点:
5 6
第2个点-双连通分量包含以下点:
3 4 5
第3个点-双连通分量包含以下点:
1 2 3
*/