/* There's someone in my head but it's not me */
#include <bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef pair<int, int> pii;
#define SZ(x) (int) x.size()
#define F first
#define S second
const int N = 1e6 + 1;
int n, rt[N], P[N], C[N], M[N], out[N], w[N], R; ll A[N], ret[N], H[N], S[N];
vector<int> adj[N], vec; vector<vector<int>> comps;
int Find(int v) {
return !rt[v] ? v : rt[v] = Find(rt[v]);
}
int Union(int v, int u) {
u = Find(u), v = Find(v);
if (u == v) return 1;
rt[u] = v;
return 0;
}
void preDFS(int v, int tr, int p = -1) {
vec.push_back(v);
for (int id : adj[v]) {
int u = out[id] ^ id ^ v;
if ((!C[u] || u == tr) && u != p)
H[u] = H[v] + w[id], preDFS(u, tr, v);
}
}
void cycDFS(int v, int tr, int p = -1) {
P[v] = p; vec.push_back(v);
for (int id : adj[v]) {
if (id != tr) {
int u = out[id] ^ id ^ v;
if (u != p) cycDFS(u, tr, v);
}
}
}
ll CYC(int v, int p = -1) {
ll ans = -1e18;
for (int id : adj[v]) {
int u = out[id] ^ id ^ v;
if (C[u] && u != p && u != R) {
H[u] = H[v] + w[id], ans = CYC(u, v);
break;
}
}
ret[Find(v)] = max(ret[Find(v)], ans + A[v] - H[v]);
return max(ans, A[v] + H[v]);
}
ll CYC2(int v, int p = -1) {
ll ans = -1e18;
for (int id : adj[v]) {
int u = out[id] ^ id ^ v;
if (C[u] && u != p && u != R) {
ans = CYC2(u, v);
break;
}
}
ret[Find(v)] = max(ret[Find(v)], ans + A[v] + H[v]);
return max(ans, S[Find(v)] - H[v] + A[v]);
}
int main() {
scanf("%d", &n);
vector<int> bads;
for (int i = 1; i <= n; i++) {
scanf("%d%d", &out[i], &w[i]);
adj[i].push_back(i);
adj[out[i]].push_back(i);
if (Union(i, out[i]))
bads.push_back(i);
}
for (int id : bads) {
vec = {};
cycDFS(out[id], id);
comps.push_back(vec);
int v = id;
while (~v)
C[v] = 1, v = P[v];
}
bads = {}, bads.shrink_to_fit();
for (int i = 1; i <= n; i++) {
if (!C[i]) continue;
vec = {}, preDFS(i, i);
for (int id : adj[i]) {
int u = out[id] ^ id ^ i;
if (C[u]) S[Find(i)] += w[id];
}
int b = -1;
for (int v : vec)
if (!~b || H[b] < H[v]) b = v;
A[i] = H[b]; vec = {};
H[b] = 0, preDFS(b, i);
b = -1;
for (int v : vec)
if (!~b || H[b] < H[v]) b = v;
ret[Find(i)] = max(ret[Find(i)], H[b]);
}
ll res = 0;
for (auto &c : comps) {
S[Find(c.back())] >>= 1;
for (int &v : c)
if (C[v]) R = v;
H[R] = 0;
CYC(R), CYC2(R);
res += ret[Find(c.back())];
}
printf("%lld\n", res);
return 0;
}
Compilation message
islands.cpp: In function 'int main()':
islands.cpp:75:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
75 | scanf("%d", &n);
| ~~~~~^~~~~~~~~~
islands.cpp:78:14: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
78 | scanf("%d%d", &out[i], &w[i]);
| ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
15 ms |
23884 KB |
Output is correct |
2 |
Correct |
17 ms |
23756 KB |
Output is correct |
3 |
Correct |
16 ms |
23824 KB |
Output is correct |
4 |
Correct |
15 ms |
23852 KB |
Output is correct |
5 |
Correct |
15 ms |
23780 KB |
Output is correct |
6 |
Correct |
15 ms |
23756 KB |
Output is correct |
7 |
Correct |
17 ms |
23756 KB |
Output is correct |
8 |
Correct |
18 ms |
23796 KB |
Output is correct |
9 |
Correct |
18 ms |
23756 KB |
Output is correct |
10 |
Correct |
16 ms |
23872 KB |
Output is correct |
11 |
Correct |
15 ms |
23872 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
16 ms |
23944 KB |
Output is correct |
2 |
Correct |
16 ms |
23960 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
18 ms |
24012 KB |
Output is correct |
2 |
Correct |
18 ms |
24304 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
30 ms |
25028 KB |
Output is correct |
2 |
Correct |
41 ms |
27464 KB |
Output is correct |
3 |
Correct |
31 ms |
25496 KB |
Output is correct |
4 |
Correct |
22 ms |
24668 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
49 ms |
28948 KB |
Output is correct |
2 |
Correct |
69 ms |
32200 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
152 ms |
40248 KB |
Output is correct |
2 |
Correct |
133 ms |
44448 KB |
Output is correct |
3 |
Correct |
166 ms |
51928 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
215 ms |
53540 KB |
Output is correct |
2 |
Correct |
270 ms |
73636 KB |
Output is correct |
3 |
Correct |
313 ms |
82736 KB |
Output is correct |
4 |
Correct |
366 ms |
95656 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
433 ms |
100884 KB |
Output is correct |
2 |
Correct |
1247 ms |
124460 KB |
Output is correct |
3 |
Correct |
393 ms |
82132 KB |
Output is correct |
4 |
Correct |
523 ms |
115408 KB |
Output is correct |
5 |
Correct |
530 ms |
115984 KB |
Output is correct |
6 |
Correct |
1935 ms |
87156 KB |
Output is correct |
7 |
Runtime error |
513 ms |
131076 KB |
Execution killed with signal 9 |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
426 ms |
131076 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |