构造边双连通图（求桥） poj 3352

很难找到裸求割边的题，于是做了这道求割边的拓展。

可先看http://blog.csdn.net/huyuncong/archive/2011/05/12/6415915.aspx

一条边（u，v）是割边，必须满足它是树边，且low[v]>rel[u]（即v是u的孩子，且v无法不通过（u，v）到达u的祖先，这样，一旦（u，v）不存在，v便与u的祖先割开），通过dfs，我们便可求出所有割边（桥）。

求边双连通分量比求点双连通分量简单，将所有桥从图中删去，各个子连通分量便是边双连通分量。

简单提一下求点双连通分量，与求强连通分量类似，我们需要用栈来维护，每dfs点x，将与x相连的树边与回边加入栈，直到遇到rel[u]<=low[v]，即u为割点，将栈中元素出栈至遇到（u，v），出栈元素构成点双连通分量，割点同时属于多个分量，其余点属于同一分量，对于点双连通，我暂时还未实现，大家可以去http://www.byvoid.com/blog/biconnect/zh-hans/，其中对于tarjan描述的很清楚。

至于构造双连通图，是指在有桥的图中，加边使之成为双连通图，具体做法是求出边双连通分量，缩点，使之成为一棵树，要加的边的个数即（叶子数+1）/2，加边集合为不断向lca最远叶子连边，收缩路径上点，直至整个图缩完，恰好（叶子数+1）/2次（具体证明我不知道，如有人知道，欢迎评论）。'

poj3352 只需求出边数，据说与poj3177代码一模一样，我鲜有的1A了一回。

var next,sora,tail:array[1..10000]of longint; f,g:array[1..1000,1..1000]of boolean; b,rel,low:array[1..1000]of longint; v:array[1..1000]of boolean; n,m,ans,ss,time:longint; procedure origin; var i:longint; begin fillchar(next,sizeof(next),0);fillchar(sora,sizeof(sora),0); fillchar(tail,sizeof(tail),0);ss:=n; for i:=1 to n do tail[i]:=i end; procedure link(x,y:longint); begin inc(ss);next[tail[x]]:=ss;tail[x]:=ss;sora[ss]:=y; inc(ss);next[tail[y]]:=ss;tail[y]:=ss;sora[ss]:=x; f[x,y]:=true;f[y,x]:=true end; function min(x,y:longint):longint; begin if x<y then exit(x) else exit(y) end; procedure dfs(x,y:longint); var i,ne:longint; begin inc(time);rel[x]:=time;low[x]:=time; i:=x; while next[i]<>0 do begin i:=next[i];ne:=sora[i]; if (rel[ne]<>0)and(ne<>y) then low[x]:=min(low[x],rel[ne]) else if rel[ne]=0 then begin dfs(ne,x); low[x]:=min(low[x],low[ne]); if low[ne]>rel[x] then begin f[x,ne]:=false;f[ne,x]:=false end; end end end; function find(x:longint):longint; begin if b[x]<>x then b[x]:=find(b[x]); find:=b[x] end; procedure link2(x,y:longint); begin next[tail[x]]:=y;tail[x]:=tail[y] end; procedure dfs2(x,y:longint); var i,ne,fa:longint; begin v[x]:=false;i:=x; while next[i]<>0 do begin i:=next[i];ne:=sora[i]; if ne<>0 then if f[x,ne] and v[ne] then begin fa:=find(ne); if fa<>y then begin b[fa]:=y;g[ne,y]:=false;g[y,ne]:=false; link2(y,fa) end; dfs2(ne,y) end end end; procedure init; var i,j,ne,x,y,sum:longint; begin readln(n,m); origin; fillchar(f,sizeof(f),false); for i:=1 to m do begin readln(x,y);link(x,y) end; fillchar(rel,sizeof(rel),0);fillchar(low,sizeof(low),0); time:=0;dfs(1,1); fillchar(v,sizeof(v),true);fillchar(b,sizeof(b),0); for i:=1 to n do b[i]:=i; fillchar(g,sizeof(g),true); for i:=1 to n do if v[i] then dfs2(i,i); ans:=0; for i:=1 to n do if b[i]=i then begin j:=i;sum:=0;g[i,i]:=false; while next[j]<>0 do begin j:=next[j];ne:=sora[j]; if (ne<>0) then if g[i,ne] then begin inc(sum); if sum>1 then break end end; if sum=1 then inc(ans) end; writeln((ans+1)div 2) end; begin assign(input,'3352.in');reset(input); init; close(input) end.