Submission #852809

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
852809	2023-09-22T18:16:07 Z	lorenzoferrari	Islands (IOI08_islands)	C++17	90 / 100	1233 ms	131072 KB

islands

#include "bits/stdc++.h"
using namespace std;
using LL = long long;

static constexpr int N = 1e6+5;

vector<array<int, 2>> edges;
vector<vector<int>> adj;
vector<array<int, 2>> par;
bitset<N> vis, cyc;
int cycle_node = -1;
int back_edge = -1;

void dfs(int v, int pi) {
    vis[v] = true;
    for (auto i : adj[v]) {
        auto [uv, w] = edges[i];
        int u = uv^v;
        if (i == pi) continue;   // deal with cycles of length 2
        if (!vis[u]) {
            par[u] = {v, w};
            dfs(u, i);
        } else if (vis[u] && cycle_node == -1) {
            cycle_node = v;
            par[u] = {v, w};
            back_edge = i;
        }
    }
}

void find_diam(int v, int p, LL dist, array<LL, 2>& max_dist) {
    max_dist = max(max_dist, array<LL, 2>{dist, v});
    for (auto i : adj[v]) {
        auto [uv, w] = edges[i];
        int u = uv^v;
        if (u == p || i == back_edge) continue;
        find_diam(u, v, dist + w, max_dist);
    }
}

LL height(int v, int p) {
    LL ans = 0;
    for (auto i : adj[v]) {
        auto [uv, w] = edges[i];
        int u = uv^v;
        if (u == p || i == back_edge) continue;
        if (cyc[u]) continue;
        ans = max(ans, w + height(u, v));
    }
    return ans;
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);

    int n; cin >> n;
    adj.resize(n+1);
    par.resize(n+1);
    edges.resize(n+1);

    for (int v = 1, u, w; v <= n; ++v) {
        cin >> u >> w;
        edges[v] = {u^v, w};
        adj[v].push_back(v);
        adj[u].push_back(v);
    }

    LL ans = 0;
    for (int v = 1; v <= n; ++v) {
        if (!vis[v]) {
            // first dfs
            dfs(v, -1);

            // diameter
            array<LL, 2> max_dist{0, v};
            find_diam(v, -1, 0, max_dist);
            int u = max_dist[1];
            max_dist = {0, u};
            find_diam(u, -1, 0, max_dist);
            LL diameter = max_dist[0];

            // path traversing the cycle
            LL a = cycle_node;
            LL cycle_length = 0;
            while (par[a][0] != cycle_node) {
                cyc[a] = true;
                cycle_length += par[a][1];
                a = par[a][0];
            }
            cyc[a] = true;
            cycle_length += par[a][1];

            a = cycle_node;
            LL max1 = -1e9, max2 = -1e9;
            LL prv1 = height(a, -1);
            LL prv2 = prv1;
            while (par[a][0] != cycle_node) {
                auto [b, w] = par[a];
                prv1 += w;
                prv2 -= w;
                LL h = height(b, -1);
                max1 = max(max1, prv1 + h);
                max2 = max(max2, prv2 + h);
                prv1 = max(prv1, h);
                prv2 = max(prv2, h);
                a = b;
            }
            LL cycle_path = max(max1, cycle_length + max2);

            // increase the answer by the maximum path in the current cc
            ans += max(diameter, cycle_path);

            cycle_node = -1;    // reset
            back_edge = -1;
        }
    }

    cout << ans << "\n";
}

Subtask #1
33.0 / 33.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	0 ms	348 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	348 KB	Output is correct
7	Correct	0 ms	348 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct
11	Correct	1 ms	348 KB	Output is correct

Subtask #2
3.0 / 3.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	1 ms	344 KB	Output is correct

Subtask #3
4.0 / 4.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	604 KB	Output is correct
2	Correct	2 ms	860 KB	Output is correct

Subtask #4
10.0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	6 ms	1628 KB	Output is correct
2	Correct	13 ms	4184 KB	Output is correct
3	Correct	8 ms	1832 KB	Output is correct
4	Correct	4 ms	1192 KB	Output is correct

Subtask #5
10.0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	17 ms	5720 KB	Output is correct
2	Correct	25 ms	9148 KB	Output is correct

Subtask #6
10.0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	62 ms	16492 KB	Output is correct
2	Correct	68 ms	23792 KB	Output is correct
3	Correct	82 ms	32196 KB	Output is correct

Subtask #7
10.0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	98 ms	30292 KB	Output is correct
2	Correct	141 ms	54604 KB	Output is correct
3	Correct	169 ms	67296 KB	Output is correct
4	Correct	200 ms	81916 KB	Output is correct

Subtask #8
10.0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	199 ms	74836 KB	Output is correct
2	Correct	824 ms	118684 KB	Output is correct
3	Correct	204 ms	58984 KB	Output is correct
4	Correct	278 ms	99668 KB	Output is correct
5	Correct	282 ms	99244 KB	Output is correct
6	Correct	1233 ms	69200 KB	Output is correct
7	Correct	311 ms	131072 KB	Output is correct

Subtask #9
0 / 10.0

#	Verdict	Execution time	Memory	Grader output
1	Runtime error	243 ms	131072 KB	Execution killed with signal 9
2	Halted	0 ms	0 KB	-

Submission #852809

Compilation message