oj.uz

Submission #489250

# Submission time Handle Problem Language Result Execution time Memory
489250 2021-11-21T18:15:35 Z dxz05 Islands (IOI08_islands) C++14
90 / 100
 1207 ms 131076 KB 
//#pragma GCC optimize("Ofast,O2,O3,unroll-loops")
//#pragma GCC target("avx2")

#include <iostream>
#include <vector>
#include <bitset>
#include <stack>
#include <algorithm>
#include <chrono>
#include <random>
using namespace std;

#define SZ(s) ((int)s.size())
#define all(x) (x).begin(), (x).end()
#define lla(x) (x).rbegin(), (x).rend()

#define MP make_pair

typedef long long ll;
const double eps = 0.000001;
const int MOD = 998244353;
const int INF = 1000000101;
const long long LLINF = 1223372000000000555;
const int N = 1e6 + 1;
const int M = 222;

mt19937 rng(chrono::high_resolution_clock::now().time_since_epoch().count());

int f[N], w[N];

vector<int> cycle;
stack<int> st;
bitset<N> used;

inline void find_cycle(int x){
    while (!st.empty()) st.pop();
    cycle.clear();
    while (true){
        st.push(x);
        used[x] = true;
        if (used[f[x]]){
            while (!st.empty()){
                int y = st.top(); st.pop();
                cycle.push_back(y);
                if (y == f[x]) break;
            }
            return;
        }
        x = f[x];
    }
}

int parent[N];
int find(const int x){
    return (x == parent[x] ? x : parent[x] = find(parent[x]));
}

vector<int> g[N];

inline int A(const int &u, const int &v){
    if (f[u] == v && f[v] != u) return w[u];
    if (f[u] != v && f[v] == u) return w[v];
    return max(w[u], w[v]);
}

int X; ll D;

void dfs(const int v, const int p, const ll d, const int c1, const int c2, const bool flag){
    if (d > D) D = d, X = v;
    for (int u : g[v]){
        if (u != p && u != c1 && u != c2) dfs(u, v, d + A(u, v), c1, c2, flag);
    }
    if (flag) g[v].clear();
}

ll arr[N];

inline ll component(const int root){
    find_cycle(root);

    const int n = cycle.size();

    ll ans = 0;
    for (int i = 0; i < n; i++){
        int v = cycle[i];
        D = -1;
        dfs(v, v, 0, cycle[(i + n - 1) % n], cycle[(i + 1) % n], false);

        ans = max(ans, D);
        arr[v] = D;

        D = -1;
        dfs(X, X, 0, cycle[(i + n - 1) % n], cycle[(i + 1) % n], true);

        ans = max(ans, D);
    }

    if (n == 2){
        int u = cycle[0], v = cycle[1];
        ans = max(ans, arr[u] + arr[v] + A(u, v));
        return ans;
    }

    vector<ll> a(n);
    a[0] = {arr[cycle[0]]};
    ll sum = 0;
    for (int i = 1; i < n; i++){
        sum += A(cycle[i - 1], cycle[i]);
        a[i] = arr[cycle[i]] + sum;
    }
    sum += A(cycle[0], cycle.back());

    ans = max(ans, arr[cycle[0]] + *max_element(a.begin() + 1, a.end()));

    vector<ll> sf(n, a.back());
    for (int i = n - 2; i >= 0; i--) sf[i] = max(sf[i + 1], a[i]);

    ll add = 0;
    int v, u;
    ll mx = 0;
    for (int i = 1; i < n; i++){
        v = cycle[i], u = cycle[i - 1];
        add -= A(v, u);
        mx = max(mx, a[i - 1] + sum);

        ans = max(ans, arr[v] + add + mx);
        if (i + 1 < n) ans = max(ans, arr[v] + add + sf[i + 1]);
        mx = max(mx, a[i]);

    }

    return ans;
}

int main() {
    ios_base::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr);
    int n;
    cin >> n;

    for (int i = 1; i <= n; i++) parent[i] = i;

    for (int i = 1; i <= n; i++){
        cin >> f[i] >> w[i];
        if (rng() & 1) parent[find(i)] = find(f[i]); else
            parent[find(f[i])] = find(i);
        if (f[f[i]] != i){
            g[i].push_back(f[i]);
            g[f[i]].push_back(i);
        }
    }

    bitset<N> comp;

    ll ans = 0;
    for (int i = 1; i <= n; i++){
        if (!comp[find(i)]){
            comp[find(i)] = true;
            ans += component(i);
        }
    }

    cout << ans;


    return 0;
}

Subtask #1
33.0 / 33.0

# Verdict Execution time Memory Grader output
1 Correct 12 ms 23884 KB Output is correct
2 Correct 12 ms 23884 KB Output is correct
3 Correct 13 ms 23884 KB Output is correct
4 Correct 13 ms 23884 KB Output is correct
5 Correct 12 ms 23920 KB Output is correct
6 Correct 12 ms 23884 KB Output is correct
7 Correct 12 ms 23884 KB Output is correct
8 Correct 13 ms 24048 KB Output is correct
9 Correct 12 ms 23948 KB Output is correct
10 Correct 12 ms 23948 KB Output is correct
11 Correct 11 ms 23976 KB Output is correct

Subtask #2
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 12 ms 24012 KB Output is correct
2 Correct 16 ms 23980 KB Output is correct

Subtask #3
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 14 ms 24060 KB Output is correct
2 Correct 14 ms 24152 KB Output is correct

Subtask #4
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 19 ms 24792 KB Output is correct
2 Correct 26 ms 26184 KB Output is correct
3 Correct 21 ms 25060 KB Output is correct
4 Correct 16 ms 24396 KB Output is correct

Subtask #5
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 30 ms 27212 KB Output is correct
2 Correct 47 ms 29460 KB Output is correct

Subtask #6
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 78 ms 35200 KB Output is correct
2 Correct 94 ms 37400 KB Output is correct
3 Correct 101 ms 41648 KB Output is correct

Subtask #7
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 134 ms 44804 KB Output is correct
2 Correct 188 ms 56824 KB Output is correct
3 Correct 183 ms 61736 KB Output is correct
4 Correct 224 ms 69124 KB Output is correct

Subtask #8
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 292 ms 74912 KB Output is correct
2 Correct 803 ms 86280 KB Output is correct
3 Correct 255 ms 65384 KB Output is correct
4 Correct 333 ms 81824 KB Output is correct
5 Correct 333 ms 83548 KB Output is correct
6 Correct 1207 ms 68004 KB Output is correct
7 Correct 352 ms 94752 KB Output is correct

Subtask #9
0 / 10.0

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