GRL_3_B 橋

Calendar Clock iconCalendar Clock icon

AOJ

# 目次

# 問題

http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=GRL_3_B

# 解説

橋、とはグラフの辺でそれを削除するとグラフが分裂するような辺のこと.
橋を見つける手順は連結点を見つける手順とほとんど同一(その過程で両方見つかる).

DFS木を構築し、その過程でルートから順に到達時間をトップダウンで割り当て、葉から戻る時にボトムアップでLowLinkを計算していく.
到達時間はルートを0とし、未到達の頂点に1つ移動するたびに+1して割り当てていく.
LowLinkとは、以下3つのうち最小のもの.

  1. 自身の到達時間
  2. 子のLowLinkのうち最小のもの
  3. 後退辺の先の頂点への到達時間

後退辺とは、DFS木では辺として使われなかった辺のこと.
DFS木を構築する過程で通った辺を記録しておけば未使用のものがわかる.

DFS木の辺のうちで、親頂点の到達時間よりも子のLowLinkのほうが大きければその子は他に頼れる辺がないということなので橋である.

# 計算量

O(V+E)O(V+E)

# 解答

// C++ 14
#include <iostream>
#include <string>
#include <vector>
#include <list>
#include <algorithm>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <unordered_map>
#include <math.h>

#define ll long long
#define Int int
#define loop(x, start, end) for(Int x = start; x < end; x++)
#define loopdown(x, start, end) for(int x = start; x > end; x--)
#define rep(n) for(int x = 0; x < n; x++)
#define span(a,x,y) a.begin()+x,a.begin()+y
#define span_all(a) a.begin(),a.end()
#define len(x) (x.size())
#define last(x) (*(x.end()-1))

using namespace std;

#define MAX_V 10000

vector<Int> graph[MAX_V];
bool used_v[MAX_V],used_e[MAX_V][MAX_V];
Int ord[MAX_V],lowlink[MAX_V];
Int k=0;
vector<pair<Int, Int> > bridges;

bool compare (pair<Int, Int> a, pair<Int, Int> b) {
  if (a.first != b.first)  return a.first < b.first;
  if (a.second != b.second) return a.second < b.second;
  return false;
}

void dfs(Int v, Int parent) {
  used_v[v] = true;
  ord[v]=lowlink[v]=k++;
  
  for (Int u: graph[v]) {
    if(!used_v[u]) {
      used_e[v][u] = true;
      dfs(u, v);
      lowlink[v] = min(lowlink[v], lowlink[u]);
      if (ord[v] < lowlink[u]) {
        if (v < u) bridges.push_back(make_pair(v, u));
        else bridges.push_back(make_pair(u, v));
      }
    }
    else if(!used_e[u][v]) {
      lowlink[v] = min(lowlink[v], ord[u]);
    }
  }
}

Int main(void) {
  Int n, e, u, v;
  cin >> n >> e;
  loop(i,0,e) {
    cin >> u >> v;
    graph[u].push_back(v);
    graph[v].push_back(u);
  }
  dfs(0, -1);
  sort(bridges.begin(), bridges.end(), compare);
  for (auto b: bridges) cout << b.first << ' ' << b.second << endl;
}

リモートフリーランス。ウェブサービス、スマホアプリエンジニア。
東アジアを拠点に世界を移動しながら活動してます!

お仕事のご依頼・お問い合わせはこちら

コメント