summaryrefslogtreecommitdiff
path: root/graph/articulationPoints.cpp
blob: 0df1eeb693a4f1327c277cb30dc390aad749393e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
vector<vector<edge>> adjlist;
vector<int> num;
int counter, rootCount, root;
vector<bool> isArt;
vector<edge> bridges, st;
vector<vector<edge>> bcc;

int dfs(int v, int parent = -1) {
	int me = num[v] = ++counter, top = me;
	for (edge& e : adjlist[v]) {
		if (e.id == parent){}
		else if (num[e.to]) {
			top = min(top, num[e.to]);
			if (num[e.to] < me) st.push_back(e);
		} else {
			if (v == root) rootCount++;
			int si = sz(st);
			int up = dfs(e.to, e.id);
			top = min(top, up);
			if (up >= me) isArt[v] = true;
			if (up > me) bridges.push_back(e);
			if (up <= me) st.push_back(e);
			if (up == me) {
				bcc.emplace_back();
				while (sz(st) > si) {
					bcc.back().push_back(st.back());
					st.pop_back();
	}}}}
	return top;
}

void findArticulationPoints() {
	counter = 0;
	num.assign(sz(adjlist), 0);
	isArt.assign(sz(adjlist), false);
	bridges.clear();
	st.clear();
	bcc.clear();
	for (int v = 0; v < sz(adjlist); v++) {
		if (!num[v]) {
			root = v;
			rootCount = 0;
			dfs(v);
			isArt[v] = rootCount > 1;
}}}