洛谷P1656 炸铁路 （求割边）

用tarjan变种求割边的模板题

其实还可以求出所有的边双（用栈），但本题不需要求。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 const int N=1e5+10;
 4 int head[N],nxt[N<<1],to[N<<1],tot;
 5 pair<int,int> edge[N];
 6 int dfn[N],low[N],top/*st[N]*/;
 7 int cnt,idx,anscnt=0/*bel[N]*/;
 8 //st是用来提取连通块的栈
 9 //bel表示每个点所属连通块（这道题不需要用）
10
11 void add(int a,int b){
12     nxt[++tot]=head[a];head[a]=tot;to[tot]=b;
13     nxt[++tot]=head[b];head[b]=tot;to[tot]=a;
14 }
15
16 bool cmp(pair<int,int>x,pair<int,int>y){
17     if(x.first==y.first) return x.second<y.second;
18     else return x.first<y.first;
19 }
20
21 void dfs(int u,int fa){
22     low[u]=dfn[u]=++cnt;
23     //st[++top]=u;
24     for(int i=head[u];i;i=nxt[i]){
25         int v=to[i];
26         if(v!=fa){
27             if(!dfn[v]){//v没有被访问过
28                 dfs(v,u);
29                 low[u]=min(low[u],low[v]);
30                 if(low[v]>dfn[u])//是割边
31                     edge[anscnt++]=make_pair(min(u,v),max(u,v));//存入割边
32             }
33             else low[u]=min(low[u],dfn[v]);
34             //v被访问过了，更新u的low值
35         }
36         /*if(low[u]==dfn[u]){        用来求边双连通分量
37             int v;++idx;
38             do{
39                 v=st[top--];
40                 bel[v]=idx;
41             }while(v!=u);
42         }*/
43     }
44 }
45
46 int main(){
47     int n,m;
48     scanf("%d%d",&n,&m);
49     for(int i=1;i<=m;i++){
50         int a,b;
51         scanf("%d%d",&a,&b);
52         add(a,b);
53     }
54     for(int i=1;i<=n;i++)
55         if(!dfn[i]) dfs(i,0);
56     sort(edge,edge+anscnt,cmp);//按照题目要求给答案排序
57     for(int i=0;i<anscnt;i++){
58         int a=edge[i].first,b=edge[i].second;
59         printf("%d %d\n",a,b);
60     }
61 }