Алгоритмы на графах / dfs / bicon

{ graphs by boris }
{ depth-first search (biconnected components) }

const
  maxn = 250;

type
  stack = array [0..maxn] of word;      { stack without doubles }
                                        { stack[0] is stack pointer }
procedure init_stack(var s: stack);
begin
  s[0] := 0;
end;

procedure push_to_stack(var s: stack; w: word);
begin
  inc(s[0]);
  s[s[0]] := w;
end;

procedure print_stack(const s: stack);
var
  i: word;
begin
  for i := 1 to s[0] do
    write(s[i], ' ');
end;

var
  a: array [1..maxn, 1..maxn] of boolean;
  pre, low, p: array [1..maxn] of word;
  s: stack;
  pre_c, n: longint;

procedure init;
var
  i, x, y, nn: longint;
begin
  fillchar(a, sizeof(a), false);
  fillchar(pre, sizeof(pre), 0);
  fillchar(p, sizeof(p), 0);
  fillchar(low, sizeof(low), 0);
  pre_c  := 1;

  assign(input, 'graph.in');
  reset(input);

  read(n);
  read(nn);

  for i := 1 to nn do
  begin
    read(x, y);
    a[x, y] := true;
    a[y, x] := true;                    { Ґб«Ё ­Ґ®аЁҐ­вЁа®ў ­­л© Ја д }
  end;
end;

procedure print;
var
  i, j: integer;
begin
  writeln;
  writeln('number of vertex : ', n);
  writeln('adjacency matrix');
  for i := 1 to n do
  begin
    for j := 1 to n do
      write(ord(a[i, j]));
    writeln;
  end;
end;

function min(a, b: longint): longint;
begin
  if a > b then min := b else min := a;
end;

procedure dfs(v: longint);
var
  i: integer;
begin
  pre[v] := pre_c;
  low[v] := pre_c;
  inc(pre_c);

  for i := 1 to n do
    if a[v, i] then
      if pre[i] = 0 then
      begin
        push_to_stack(s, i);
        p[i] := v;
        dfs(i);
        low[v] := min(low[v], low[i]);
        if low[i] >= pre[v] then
        begin
          print_stack(s);
        end;
      end
      else if (p[v] <> i) and (pre[i] < pre[v]) then
      begin
        push_to_stack(s, i);
        low[v] := min(low[v], pre[i]);
      end;
end;

var
  i: longint;

begin
  init;

  writeln('depth-first search (biconnected components)');
  for i := 1 to n do
    if pre[i] = 0 then
      dfs(i);

  writeln;

  write('pre order : ');
  for i := 1 to n do
    write(pre[i], ' ');

  writeln;

  write('low : ');
  for i := 1 to n do
    write(low[i], ' ');

  print;
end.