-
Notifications
You must be signed in to change notification settings - Fork 5
Expand file tree
/
Copy pathkosaraju algorithm for SCC.cpp
More file actions
68 lines (52 loc) · 1.06 KB
/
kosaraju algorithm for SCC.cpp
File metadata and controls
68 lines (52 loc) · 1.06 KB
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
#include<bits/stdc++.h>
using namespace std;
const int N = 1e5;
vector<int> gr[N], grr[N];
int vis[N], col[N];
vector<int> order;
vector<int> components[N];
// topological sorting
void dfs1(int cur) {
vis[cur] = 1;
for (auto x : gr[cur]) {
if (!vis[x]) dfs1(x);
}
order.push_back(cur);
}
void dfs2(int cur, int comp) {
vis[cur] = 1;
col[cur] = comp;
components[comp].push_back(cur);
for (auto x : grr[cur]) {
if (!vis[x]) dfs2(x, comp);
}
}
int main()
{
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
int n, m;
cin >> n >> m;
for (int i = 0; i < m; i++) {
int x, y;
cin >> x >> y;
gr[x].push_back(y);
grr[y].push_back(x);
}
for (int i = 1; i <= n; i++) {
if (!vis[i]) dfs1(i);
}
reverse(order.begin(), order.end());
memset(vis, 0, sizeof(vis));
int comp = 1;
for (auto x : order) {
// cout << x << " ";
if (!vis[x]) dfs2(x, comp++);
}
for (int i = 1; i <= n; i++) {
cout << i << " " << col[i] << '\n';
}
cout << "Total strongly components are -> " << comp - 1 << '\n';
// complexity O(n+m)
return 0;
}