双连通分量是图论中的一个重要模型,它的求解方式和割顶的差不多。
代码里用一个栈来保存双连通分量中的点。
注意,这里所说的双连通分量是点-双连通分量;
下面是核心代码:
class edge
{
public:
int u;
int v;
};
bool iscut[maxn];
vector<int>G[maxn];
vector<int>bcc[maxn];
stack<edge>S;
int n,m,low[maxn],pre[maxn],dfs_clock=0,bcc_cnt=0;
int bccno[maxn];
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[bcc_cnt].clear();
for(;;)
{
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)
{
S.push(e);
lowu = min(lowu,pre[v]);
}
}
if(fa < 0&&child==1)
iscut[u] = false;
low[u] = lowu;
return lowu;
}
void find_bcc(int n)
{
memset(pre,0,sizeof(pre));
memset(bccno,0,sizeof(bccno));
memset(iscut,0,sizeof(iscut));
bcc_cnt = dfs_clock = 0;
for(int i = 0 ; i < n ; i++ )
if(!pre[i])
dfs(i,-1);
}
代码有点长,其中很大一部分是空格、回车之类的,为的是让代码更美观。
PS:博主有强迫症,不喜欢看难看的代码;