Pagini recente » Cod sursa (job #1375073) | Cod sursa (job #611826) | Cod sursa (job #2376660) | Cod sursa (job #1163424) | Cod sursa (job #2194665)
// Copyright 2018 Popescu Alexandru Gabriel <[email protected]>
// Copyright 2018 Darius Neatu <[email protected]>
// Componente Tare Conexe (CTC) - Strongly Connected Components (SCC)
// Kosaraju - O(n + m)
// https://infoarena.ro/problema/ctc
#include <bits/stdc++.h>
#define NMAX 100010
using namespace std;
ifstream f("ctc.in");
ofstream g("ctc.out");
class Task {
public:
void solve() {
read_input();
kosaraju();
print_output();
}
private:
int n, m, sol = 0;
vector<int> adj[NMAX]; // graful initial
vector<int> adj_t[NMAX]; // graful transpus
vector<vector<int> > CTC; // retin solutia
vector<int> visited;
vector<int> Tsort; // retin nodurile crescator dupa timpul de finalizare, nu folosim stiva
void read_input() {
f >> n >> m;
visited = vector<int>(n + 1, 0);
Tsort = vector<int>(n + 1, 0);
int x, y;
for (int i = 1; i <= m; ++i) {
f >> x >> y;
adj[x].push_back(y);
adj_t[y].push_back(x);
}
}
void kosaraju() {
for (int i = 1; i <= n; ++i) {
if (!visited[i]) {
dfs(i);
}
}
for (int i = n; i >= 1; --i) {
if (visited[Tsort[i]]) {
vector<int> current_ctc;
dfs_t(Tsort[i], sol, current_ctc);
CTC.push_back(current_ctc);
++sol;
}
}
}
void dfs(int node) {
visited[node] = 1;
for (auto &vecin : adj[node]) {
if (!visited[vecin]) {
dfs(vecin);
}
}
Tsort[++Tsort[0]] = node; //Tsort[0] functioneaza pe post de counter
}
void dfs_t(int node, int &sol, vector<int> ¤t_ctc) {
visited[node] = 0;
current_ctc.push_back(node);
for (auto &vecin : adj_t[node]) {
if (visited[vecin]) {
dfs_t(vecin, sol, current_ctc);
}
}
}
void print_output() {
// pentru verificare pe infoarena trebuie doar modificata afisarea
//std::cout << "Graful are " << sol << " componente tare conexe\n";
g << sol << '\n';
for (int i = 0; i < sol; ++i) {
//std::cout << "Componenta tare conexa [" << i + 1 << "]: ";
for (auto &node : CTC[i]) {
g << node << " ";
}
g << '\n';
}
}
};
int main() {
Task *task = new Task();
task->solve();
delete task;
return 0;
}