The OPTGRAPH Procedure

Biconnected Components and Articulation Points

A biconnected component of a graph $G = (N,A)$ is a connected subgraph that cannot be broken into disconnected pieces by deleting any single node (and its incident links). An articulation point is a node of a graph whose removal would cause an increase in the number of connected components. Articulation points can be important when you analyze any graph that represents a communications network. Consider an articulation point $i \in N$ which, if removed, disconnects the graph into two components $C^1$ and $C^2$. All paths in G between some nodes in $C^1$ and some nodes in $C^2$ must pass through node i. In this sense, articulation points are critical to communication. Examples of where articulation points are important are airline hubs, electric circuits, network wires, protein bonds, traffic routers, and numerous other industrial applications.

In PROC OPTGRAPH, you can find biconnected components and articulation points of an input graph by invoking the BICONCOMP statement. This algorithm works only with undirected graphs.

The results for the biconnected components algorithm are written to the output links data set that is specified in the OUT_LINKS= option in the PROC OPTGRAPH statement. For each link in the links data set, the variable biconcomp identifies its component. The component identifiers are numbered sequentially starting from 1. The results for the articulation points are written to the output nodes data set that is specified in the OUT_NODES= option in the PROC OPTGRAPH statement. For each node in the nodes data set, the variable artpoint is either 1 (if the node is an articulation point) or 0 (otherwise).

The biconnected components algorithm reports status information in a macro variable called _OPTGRAPH_BICONCOMP_. See the section Macro Variable _OPTGRAPH_BICONCOMP_ for more information about this macro variable.

The algorithm that PROC OPTGRAPH uses to compute biconnected components is a variant of depth-first search (Tarjan 1972). This algorithm runs in time $O(|N|+|A|)$ and therefore should scale to very large graphs.

Biconnected Components of a Simple Undirected Graph

This section illustrates the use of the biconnected components algorithm on the simple undirected graph G that is shown in Figure 1.17.

Figure 1.17: A Simple Undirected Graph G

A Simple Undirected Graph


The undirected graph G can be represented by the following links data set LinkSetInBiCC:

data LinkSetInBiCC;
   input from $ to $ @@;
   datalines;
A B  A F  A G  B C  B D
B E  C D  E F  G I  G H
H I
;

The following statements calculate the biconnected components and articulation points and output the results in the data sets LinkSetOut and NodeSetOut:

proc optgraph
   data_links = LinkSetInBiCC
   out_links  = LinkSetOut
   out_nodes  = NodeSetOut;
   biconcomp;
run;

The data set LinkSetOut now contains the biconnected components of the input graph, as shown in Figure 1.18.

Figure 1.18: Biconnected Components of a Simple Undirected Graph

from to biconcomp
A B 2
A F 2
A G 4
B C 1
B D 1
B E 2
C D 1
E F 2
G I 3
G H 3
H I 3



In addition, the data set NodeSetOut contains the articulation points of the input graph, as shown in Figure 1.19.

Figure 1.19: Articulation Points of a Simple Undirected Graph

node artpoint
A 1
B 1
F 0
G 1
C 0
D 0
E 0
I 0
H 0



The biconnected components are shown graphically in Figure 1.20 and Figure 1.21.

Figure 1.20: Biconnected Components $C^1$ and $C^2$

$C^1 = \{ B,C,D\} $

$C^2 = \{ A,B,E,F\} $

$\quad $

biconcomp1_1

biconcomp1_2


Figure 1.21: Biconnected Components $C^3$ and $C^4$

$C^3 = \{ G,H,I\} $

$C^4 = \{ A,G\} $

$\quad $

biconcomp1_4

biconcomp1_3


For a more detailed example, see Articulation Points in a Terrorist Network.