Submission #1055505

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

islands

#include <bits/stdc++.h>
 
using namespace std;
 
const int N = 2000002;
 
int a[N], l[N], cyc[N];
 
long long d[N], dp[N], mx[N];
 
vector<int> g[N];
 
long long Solve(int v, int pv) {
  long long mx = 0;
  for (int u : g[v]) {
    if (d[u] != -2 || u == pv) {
      continue;
    }
    mx = max(mx, Solve(u, v));
    mx = max(mx, dp[v] + dp[u] + l[u]);
    dp[v] = max(dp[v], dp[u] + l[u]);
  }
  vector<int>().swap(g[v]);
  return mx;
}
 
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;
  }
  long long c, t;
  int low, high, pos, mid;
  for (int i = 0; i < n; i++) {
    if (d[i] != -1) {
      continue;
    }
    int sz = 0;
    cyc[sz++] = i;
    d[i] = 0;
    for (int b = 0; b < sz; b++) {
      int x = cyc[b];
      if (d[a[x]] == -1) {
        d[a[x]] = 0;
        cyc[sz++] = a[x];
      }
      for (int y : g[x]) {
        if (d[y] == -1) {
          d[y] = 0;
          cyc[sz++] = y;
        }
      }
    }
    for (int x = 0; x < sz; x++) {
      d[cyc[x]] = -1;
    }
    int x = i;
    while (d[x] == -1) {
      d[x] = 0;
      x = a[x];
    }
    for (int x = 0; x < sz; x++) {
      d[cyc[x]] = -2;
    }
    sz = 0;
    cyc[sz++] = x;
    int y = a[x];
    while (y != x) {
      cyc[sz++] = y;
      y = a[y];
    }
    d[cyc[0]] = 0;
    for (int i = 1; i < sz; i++) {
      int x = cyc[i - 1], y = cyc[i];
      d[y] = d[x] + l[x];
    }
    c = 0;
    for (int i = 0; i < sz; i++) {
      c = max(c, Solve(cyc[i], cyc[i]));
      c = max(c, dp[cyc[i]]);
    }
    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]]);
    }
    t = d[cyc[sz - 1]] + l[cyc[sz - 1]];
    for (int i = 0; i < sz; i++) {
      low = i + 1, high = sz - 1, pos = i;
      while (low <= high) {
        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++) {
      low = i + 1, high = sz - 1, pos = i;
      while (low <= high) {
        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]);
      }
    }
    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