Submission #501734

# Submission time Handle Problem Language Result Execution time Memory
501734 2022-01-04T12:16:46 Z 600Mihnea Islands (IOI08_islands) C++17
90 / 100
1355 ms 131076 KB
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

struct Edge {
  int a, b, cost;
  bool seen = false;
};

const int N = (int) 1e6 + 7;
const ll INF = (ll) 1e18;
int n;
int parrent[N];
Edge graph_edges[N];
vector<int> g[N], verts;
bool vis[N];
int dep[N];
int epar[N];

vector<int> geallonp(int a, int b) {
  vector<int> p1, p2;
  while (a != b) {
    if (dep[a] > dep[b]) {
      p1.push_back(a);
      a = parrent[a];
    } else {
      p2.push_back(b);
      b = parrent[b];
    }
  }
  p1.push_back(a);
  reverse(p2.begin(), p2.end());
  for (auto &x : p2) {
    p1.push_back(x);
  }
  return p1;
}

vector<int> geallonpedges(int a, int b) {
  vector<int> p1, p2;
  while (a != b) {
    if (dep[a] > dep[b]) {
      p1.push_back(epar[a]);
      a = parrent[a];
    } else {
      p2.push_back(epar[b]);
      b = parrent[b];
    }
  }
  reverse(p2.begin(), p2.end());
  for (auto &x : p2) {
    p1.push_back(x);
  }
  return p1;
}

bool blocked[N];
ll under[N];
ll bbpath[N];

void dfs(int a, int p = -1) {
  for (auto &i : g[a]) {
    int b = graph_edges[i].a + graph_edges[i].b - a;
    assert(p == parrent[a]);
    ///cout << p << " vs " << parrent[a] << "\n";
    if (blocked[b] || b == p) {
      continue;
    }
    dfs(b, a);
    bbpath[a] = max(bbpath[a], bbpath[b]);
    bbpath[a] = max(bbpath[a], under[b] + graph_edges[i].cost + under[a]);
    under[a] = max(under[a], under[b] + graph_edges[i].cost);
  }

}

ll sum[N];
bool v2[N];

ll solveComponent(int rootVertex) {
  ll sol = 0;
  verts.clear();

  queue<int> q;
  q.push(rootVertex);
  vis[rootVertex] = true;

  while (!q.empty()) {
    int a = q.front();
    q.pop();
    verts.push_back(a);
    for (auto &i : g[a]) {
      int b = graph_edges[i].a + graph_edges[i].b - a;
      if (!vis[b]) {
        dep[b] = 1 + dep[a];
        epar[b] = i;
        vis[b] = true;
        q.push(b);
        parrent[b] = a;
        graph_edges[i].seen = true;
      }
    }
  }
  vector<int> nnseeng;
  for (auto &v : verts) {
    for (auto &i : g[v]) {
      if (!graph_edges[i].seen) {
        graph_edges[i].seen = true;
        nnseeng.push_back(i);
      }
    }
  }
  {
    int index = nnseeng[0];
    int a = graph_edges[index].a;
    int b = graph_edges[index].b;
    int cost = graph_edges[index].cost;
    vector<int> path = geallonp(a, b);
    vector<int> edges = geallonpedges(a, b);
    vector<ll> value;
    edges.push_back(index);

    vector<int> ord;

    for (auto &v : path) {
      blocked[v] = true;
      q.push(v);
      parrent[v] = -1;
      v2[v] = true;
    }


    while (!q.empty()) {
      int a = q.front();
      q.pop();
      ord.push_back(a);
      for (auto &i : g[a]) {
        int b = graph_edges[i].a + graph_edges[i].b - a;
        if (!v2[b]) {
          parrent[b] = a;
          q.push(b);
          v2[b] = true;
        }
      }
    }

    for (auto &v : path) {
      dfs(v);
      sol = max(sol, bbpath[v]);
      value.push_back(under[v]);
    }

    ll total = 0;
    for (auto &i : edges) {
      total += graph_edges[i].cost;
    }

    ll cur = 0;
    for (int i = 1; i <= (int) edges.size(); i++) {
      cur += graph_edges[edges[i - 1]].cost;
      sum[i] = cur;
    }

    int y = (int) value.size();

    ll c1 = -INF, c2 = -INF;


    for (int r = 0; r < y; r++) {

      ll b1 = c1 + sum[r] + value[r];
      ll b2 = c2 + total + value[r] - sum[r];
      c1 = max(c1, value[r] - sum[r]);
      c2 = max(c2, sum[r] + value[r]);
      sol = max(sol, b1);
      sol = max(sol, b2);
    }
  }
  return sol;
}

signed main() {
  /// L < P

  ios::sync_with_stdio(0); cin.tie(0);

  ///freopen ("input2", "r", stdin);

  cin >> n;
  for (int i = 1; i <= n; i++) {
    int j, c;
    cin >> j >> c;
    graph_edges[i].a = i;
    graph_edges[i].b = j;
    graph_edges[i].cost = c;
  }
  for (int i = 1; i <= n; i++) {
    g[graph_edges[i].a].push_back(i);
    g[graph_edges[i].b].push_back(i);
  }

  ll sol = 0;

  for (int i = 1; i <= n; i++) {
    if (!vis[i]) {
      sol += solveComponent(i);
    }
  }

  cout << sol << "\n";

  return 0;
}

Compilation message

islands.cpp: In function 'll solveComponent(int)':
islands.cpp:119:9: warning: unused variable 'cost' [-Wunused-variable]
  119 |     int cost = graph_edges[index].cost;
      |         ^~~~
# Verdict Execution time Memory Grader output
1 Correct 19 ms 39500 KB Output is correct
2 Correct 20 ms 39500 KB Output is correct
3 Correct 20 ms 39444 KB Output is correct
4 Correct 20 ms 39500 KB Output is correct
5 Correct 19 ms 39456 KB Output is correct
6 Correct 18 ms 39500 KB Output is correct
7 Correct 19 ms 39464 KB Output is correct
8 Correct 19 ms 39452 KB Output is correct
9 Correct 19 ms 39484 KB Output is correct
10 Correct 18 ms 39500 KB Output is correct
11 Correct 19 ms 39500 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 18 ms 39500 KB Output is correct
2 Correct 19 ms 39460 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 19 ms 39628 KB Output is correct
2 Correct 19 ms 39864 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 24 ms 40600 KB Output is correct
2 Correct 34 ms 42520 KB Output is correct
3 Correct 29 ms 41120 KB Output is correct
4 Correct 23 ms 40296 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 45 ms 43692 KB Output is correct
2 Correct 48 ms 45716 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 101 ms 54328 KB Output is correct
2 Correct 98 ms 55508 KB Output is correct
3 Correct 113 ms 61544 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 149 ms 63328 KB Output is correct
2 Correct 178 ms 81812 KB Output is correct
3 Correct 213 ms 80132 KB Output is correct
4 Correct 232 ms 98940 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 292 ms 102524 KB Output is correct
2 Correct 894 ms 119416 KB Output is correct
3 Correct 260 ms 88428 KB Output is correct
4 Correct 354 ms 111956 KB Output is correct
5 Correct 350 ms 111068 KB Output is correct
6 Correct 1355 ms 104360 KB Output is correct
7 Correct 368 ms 129604 KB Output is correct
# Verdict Execution time Memory Grader output
1 Runtime error 324 ms 131076 KB Execution killed with signal 9
2 Halted 0 ms 0 KB -