-
Notifications
You must be signed in to change notification settings - Fork 0
/
BOTTOM.cpp
125 lines (99 loc) · 1.77 KB
/
BOTTOM.cpp
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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
#include <bits/stdc++.h>
using namespace std;
#define all(x) x.begin(), x.end()
using ll = long long;
int n, m;
vector<vector<int>> g, tg;
int qtdcomp = 0;
vector<int> vis, comp;
vector<int> todo;
void getorder(int u) {
vis[u] = 1;
for (auto v : g[u]) {
if (!vis[v]) getorder(v);
}
todo.emplace_back(u);
}
void markcomp(int u, int compid) {
comp[u] = compid;
for (auto v : tg[u]) {
if(!comp[v]) markcomp(v, compid);
}
}
vector<set<int>> cg;
void gen() {
vis.resize(n);
fill(all(vis), 0);
comp.resize(n);
fill(all(comp), 0);
todo.clear();
for (int i =0 ; i < n; i++) {
if (!vis[i]) getorder(i);
}
reverse(all(todo));
qtdcomp = 0;
for (auto u : todo) {
if (!comp[u]) {
markcomp(u, ++qtdcomp);
}
}
cg.resize(qtdcomp+1);
for (int i = 1; i <= qtdcomp; i++)cg[i].clear();
for (int i = 0; i < n; i++) {
for (auto j : g[i]) {
int ci = comp[i];
int cj = comp[j];
if (ci != cj)
cg[ci].emplace(cj);
}
}
}
void run() {
cin >> m;
g.resize(n);
tg.resize(n);
for (int i = 0; i < n; i++) {
g[i].clear();
tg[i].clear();
}
for (int i = 0; i < m; i++) {
int a, b;
cin >> a >> b;
a--, b--;
g[a].emplace_back(b);
tg[b].emplace_back(a);
}
gen();
/*
cout << "comps:\n";
for (int i = 0; i <n; i++) {
cout << i << ' ' << comp[i] << endl;
}
return;
*/
vector<int> ans;
for (int i = 0; i < n; i++) {
int ci = comp[i];
if (cg[ci].size() == 0) ans.emplace_back(i);
}
for (int i = 0; i < ans.size(); i++) {
cout << ans[i] + 1 << " \n"[(int)i==ans.size()-1];
}
if (ans.empty()) cout << '\n';
}
int main() {
ios_base::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
while (cin >> n) {
if(n)
run();
else break;
}
}
/*
* graphs
* strongly connected components
* kosajaru
* https://vjudge.net/problem/SPOJ-BOTTOM
* */