Cod sursa(job #2569428)

Utilizator radustn92Radu Stancu radustn92 Data 4 martie 2020 12:07:40
Problema Componente biconexe Scor 90
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.22 kb
#include <cstdio>
#include <vector>
#include <stack>
#include <cstring>
using namespace std;
const int NMAX = 100505;

int N, M, parent[NMAX];
vector<vector<int>> graph;

int low[NMAX], lvl[NMAX];
bool visited[NMAX];
stack<pair<int, int>> edges;
int nodeToComp[NMAX], totalComps;

void read() {
	scanf("%d%d", &N, &M);
	graph.assign(N + 1, vector<int>());
	int x, y;
	for (int edgeIdx = 0; edgeIdx < M; edgeIdx++) {
		scanf("%d%d", &x, &y);
		graph[x].push_back(y);
		graph[y].push_back(x);
	}
}

void addNodeToLatestComp(int node, int compsNo, vector<int>& newComp) {
	if (nodeToComp[node] != compsNo) {
		nodeToComp[node] = compsNo;
		newComp.push_back(node);
	}
}

void addEdgeToLatestComp(pair<int, int>& edge, int compsNo, vector<int>& newComp) {
	addNodeToLatestComp(edge.first, compsNo, newComp);
	addNodeToLatestComp(edge.second, compsNo, newComp);
}
 
void dfs(int node) {
	visited[node] = true;
	low[node] = lvl[node];
	for (auto& neighbour : graph[node]) {
		if (!visited[neighbour]) {
			parent[neighbour] = node;
			lvl[neighbour] = lvl[node] + 1;
			dfs(neighbour);

			if (low[neighbour] >= lvl[node]) {
				totalComps++;
			}
			low[node] = min(low[node], low[neighbour]);
		} else if (parent[node] != neighbour) {
			low[node] = min(low[node], low[neighbour]);
		}
	}
}

void dfsRecons(int node) {
	visited[node] = true;
	for (auto& neighbour : graph[node]) {
		if (!visited[neighbour]) {
			edges.push({node, neighbour});
			dfsRecons(neighbour);
			if (low[neighbour] >= lvl[node]) {
				vector<int> newComp;
				totalComps++;
				while (!edges.empty()) {
					pair<int, int> edge = edges.top();
					edges.pop();
					addEdgeToLatestComp(edge, totalComps, newComp);
					if (edge.first == node && edge.second == neighbour) {
						break;
					}
				}
				for (auto node : newComp) {
					printf("%d ", node);
				}
				printf("\n");
			}
		}
	}
}

void solve() {
	for (int node = 1; node <= N; node++) {
		if (!visited[node]) {
			dfs(node);
		}
	}

	printf("%d\n", totalComps);
	memset(visited, false, sizeof(visited));
	totalComps = 0;
	for (int node = 1; node <= N; node++) {
		if (!visited[node]) {
			dfsRecons(node);
		}
	}
}

int main() {
	freopen("biconex.in", "r", stdin);
	freopen("biconex.out", "w", stdout);

	read();
	solve();
	return 0;
}