Submission #1055447

#	Time	Username	Problem	Language	Result	Execution time	Memory
1055447		MilosMilutinovic	Islands (IOI08_islands)	C++14	90 / 100	428 ms	131072 KiB

islands

#include <bits/stdc++.h>

using namespace std;

const int N = 1000000;

int a[N], l[N];

long long d[N], dp[N], mx[N];

vector<int> g[N];

int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  for (int i = 0; i < n; i++) {
    cin >> a[i] >> l[i];
    --a[i];
  }
  for (int i = 0; i < n; i++) {
    g[a[i]].push_back(i);
  }
  long long res = 0;
  for (int i = 0; i < n; i++) {
    d[i] = -1;
  }
  for (int i = 0; i < n; i++) {
    if (d[i] != -1) {
      continue;
    }
    vector<int> que(1, i);
    d[i] = 0;
    for (int b = 0; b < (int) que.size(); b++) {
      int x = que[b];
      if (d[a[x]] == -1) {
        d[a[x]] = 0;
        que.push_back(a[x]);
      }
      for (int y : g[x]) {
        if (d[y] == -1) {
          d[y] = 0;
          que.push_back(y);
        }
      }
    }
    for (int x : que) {
      d[x] = -1;
    }
    int x = i;
    while (d[x] == -1) {
      d[x] = 0;
      x = a[x];
    }
    for (int x : que) {
      d[x] = -1;
    }
    int y = a[x];
    vector<int> cyc(1, x);
    while (y != x) {
      cyc.push_back(y);
      y = a[y];
    }
    d[cyc[0]] = 0;
    for (int i = 1; i < (int) cyc.size(); i++) {
      int x = cyc[i - 1], y = cyc[i];
      d[y] = d[x] + l[x];
    }
    for (int x : que) {
      if (d[x] == -1) {
        d[x] = -2;
      }
    }
    vector<int>().swap(que);
    long long c = 0;
    function<void(int, int)> Solve = [&](int v, int pv) {
      for (int u : g[v]) {
        if (d[u] != -2 || u == pv) {
          continue;
        }
        Solve(u, v);
        c = max(c, dp[v] + dp[u] + l[u]);
        dp[v] = max(dp[v], dp[u] + l[u]);
      }
    };
    for (int i : cyc) {
      Solve(i, i);
      c = max(c, dp[i]);
    }
    int sz = (int) cyc.size();
    for (int i = sz - 1; i >= 0; i--) {
      if (i + 1 < sz) {
        mx[i] = mx[i + 1];
      } else {
        mx[i] = d[cyc[i]] + dp[cyc[i]];
      }
      mx[i] = max(mx[i], d[cyc[i]] + dp[cyc[i]]);
    }
    long long t = d[cyc.back()] + l[cyc.back()];
    for (int i = 0; i < sz; i++) {
      int low = i + 1, high = sz - 1, pos = i;
      while (low <= high) {
        int mid = (low + high) >> 1;
        if (t - (d[cyc[mid]] - d[cyc[i]]) >= (d[cyc[mid]] - d[cyc[i]])) {
          pos = mid;
          low = mid + 1;
        } else {
          high = mid - 1;
        }
      }
      if (pos + 1 < sz) {
        c = max(c, -d[cyc[i]] + dp[cyc[i]] + mx[pos + 1]);
      }
    }
    for (int i = sz - 1; i >= 0; i--) {
      if (i + 1 < sz) {
        mx[i] = mx[i + 1];
      } else {
        mx[i] = -d[cyc[i]] + dp[cyc[i]];
      }
      mx[i] = max(mx[i], -d[cyc[i]] + dp[cyc[i]]);
    }
    for (int i = 0; i < sz; i++) {
      int low = i + 1, high = sz - 1, pos = i;
      while (low <= high) {
        int mid = (low + high) >> 1;
        if (t - (d[cyc[mid]] - d[cyc[i]]) >= (d[cyc[mid]] - d[cyc[i]])) {
          pos = mid;
          low = mid + 1;
        } else {
          high = mid - 1;
        }
      }
      if (pos > i) {
        c = max(c, t + d[cyc[i]] + dp[cyc[i]] + mx[i + 1]);
      }
    }
    vector<int>().swap(cyc);
    res += c;
  }
  cout << res << '\n';
  return 0;
}

• Subtask 1: 33 / 33
• Subtask 2: 3 / 3
• Subtask 3: 4 / 4
• Subtask 4: 10 / 10
• Subtask 5: 10 / 10
• Subtask 6: 10 / 10
• Subtask 7: 10 / 10
• Subtask 8: 10 / 10
• Subtask 9: 0 / 10