Submission #319099

#	Time	Username	Problem	Language	Result	Execution time	Memory
319099		VROOM_VARUN	Towns (IOI15_towns)	C++14	25 / 100	23 ms	1152 KiB

towns

/*
ID: varunra2
LANG: C++
TASK: towns
*/

#include <bits/stdc++.h>
#include "towns.h"
using namespace std;

#ifdef DEBUG
#include "lib/debug.h"
#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#define debug_arr(...) \
  cerr << "[" << #__VA_ARGS__ << "]:", debug_arr(__VA_ARGS__)
#pragma GCC diagnostic ignored "-Wsign-compare"
//#pragma GCC diagnostic ignored "-Wunused-parameter"
//#pragma GCC diagnostic ignored "-Wunused-variable"
#else
#define debug(...) 42
#endif

#define EPS 1e-9
#define IN(A, B, C) assert(B <= A && A <= C)
#define INF (int)1e9
#define MEM(a, b) memset(a, (b), sizeof(a))
#define MOD 1000000007
#define MP make_pair
#define PB push_back
#define all(cont) cont.begin(), cont.end()
#define rall(cont) cont.end(), cont.begin()
#define x first
#define y second

const double PI = acos(-1.0);
typedef long long ll;
typedef long double ld;
typedef pair<int, int> PII;
typedef map<int, int> MPII;
typedef multiset<int> MSETI;
typedef set<int> SETI;
typedef set<string> SETS;
typedef vector<int> VI;
typedef vector<PII> VII;
typedef vector<VI> VVI;
typedef vector<string> VS;

#define rep(i, a, b) for (int i = a; i < (b); ++i)
#define trav(a, x) for (auto& a : x)
#define sz(x) (int)(x).size()
typedef pair<int, int> pii;
typedef vector<int> vi;
#pragma GCC diagnostic ignored "-Wsign-compare"
// util functions

// we are rooting the tree at leaf 0

int n;
const int mx = 150;
int vals[mx][mx];
map<PII, int> mp;

// int getDistance(int i, int j) { return vals[i][j]; }

int gt(int i, int j) {
  if (i > j) swap(i, j);
  if (i == j) return 0;
  if (mp[MP(i, j)]) return mp[MP(i, j)];
  mp[MP(i, j)] = getDistance(i, j);
  return mp[MP(i, j)];
}

int rootdist(int d1, int d2, int dtot) {
  // d1 + d2 - 2 * dlca = dtot
  // (d1 + d2 - dtot)/2 = dlca
  return (d1 + d2 - dtot) / 2;
}

int leafdist(int d1, int d2, int dtot) {
  // what is the dist(d1, dlca)
  // dist(d1, dlca) = d1 - rootdist(dlca)
  return d1 - rootdist(d1, d2, dtot);
}

void calc(int from, VI& d) {
  for (int i = 0; i < n; i++) {
    if (i == from) continue;
    d[i] = gt(from, i);
  }
}

int ind1;
int ind2;

VI d0;
VI d1;
VI d2;  // we won't be using this
VI lcadists;

void init() {
  mp.clear();
  d0.clear();
  d1.clear();
  d2.clear();
  lcadists.clear();
  d0.assign(n, 0);
  d1.assign(n, 0);
  d2.assign(n, 0);
  lcadists.assign(n, 0);
}

int find(int lcdist) {
  int val1 = d0[ind1] - lcdist;
  int val2 = d1[ind2] - val1;
  int cur = max(val1, val2);
  return cur;
}

VI cop;

VI findHubs() {
  // we will use 2n - 3 queries to find the hubs/radius/center
  calc(0, d0);
  ind1 = max_element(all(d0)) - d0.begin();
  calc(ind1, d1);
  ind2 = max_element(all(d1)) - d1.begin();
  // diameter is from ind1 -> ind2
  // hub lies on the path from 0 -> ind1 && ind1 -> ind2
  // so since we rooted the tree at 0
  // the hub lies on the path from ind1 -> lca(ind1, ind2)
  for (int i = 0; i < n; i++) {
    lcadists[i] = rootdist(d0[ind1], d0[i], d1[i]);
  }
  int mx = lcadists[ind2];
  VI ret;
  int mn = 2 * INF;
  MPII mp;
  sort(all(lcadists));
  cop = lcadists;
  lcadists.resize(unique(all(lcadists)) - lcadists.begin());

  for (int i = 0; i < sz(lcadists); i++) {
    if (lcadists[i] < mx) continue;
    int cur = find(lcadists[i]);
    if (cur < mn) {
      ret.clear();
      ret.PB(lcadists[i]);
      mn = cur;
    } else if (cur == mn)
      ret.PB(lcadists[i]);
  }
  return ret;
}

int ccnt(int dep) {
  int ret = 0;
  for (auto& x : cop) {
    if (x >= dep) ret++;
  }
  return ret;
}

int chk;

bool same(int i, int j) {
  // we need to figure out the distance of the lca of i, j from the root
  return rootdist(d0[i], d0[j], gt(i, j)) > chk;
  // if its greater than check then they in the subtree of check's children
  // if its equal, they in subtree of check, but different subtree of his
  // children if its smaller, they aren't in subtree of check
}

bool balanced(int depth) {
  chk = depth;
  // we need to check if an array contains a majority using (3n/2 - 2)
  // operations http://www.cs.yale.edu/publications/techreports/tr252.pdf ^paper
  // highlights how to solve it
  VI bucket;
  VI list;

  // phase 1

  for (int i = 0; i < n; i++) {
    if (list.empty())
      list.PB(i);
    else if (same(list.back(), i))
      bucket.PB(i);
    else {
      list.PB(i);
      if (!bucket.empty()) {
        list.PB(bucket.back());
        bucket.pop_back();
      }
    }
  }

  // debug(bucket);
  // debug(list);

  // phase 2

  int T = list.back();

  while (!list.empty()) {
    if (same(list.back(), T)) {
      if (sz(list) >= 2) {
        list.pop_back();
        list.pop_back();
      } else {
        bucket.PB(list.back());
        list.pop_back();
      }
    } else {
      list.pop_back();
      if (bucket.empty()) return true;
      bucket.pop_back();
    }
  }

  // debug(bucket);
  return bucket.empty();
}

int hubDistance(int n, int sub) {
  ::n = n;
  init();
  auto use = findHubs();
  // now we need to check if its balanced or not
  int radius = find(use[0]);
  if (sz(use) == 1) {
    return (balanced(use[0]) ? radius : -radius);
  }
  if (sz(use) == 2) {
    // if hub[0] has < n/2 leaves
    int cnt = ccnt(use[0]);
    if (cnt > n / 2)
      return (balanced(use[1]) ? radius : -radius);
    else
      return (balanced(use[0]) ? radius : -radius);
  }
}

// int main() {
// #ifndef ONLINE_JUDGE
//   freopen("towns.in", "r", stdin);
//   freopen("towns.out", "w", stdout);
// #endif
//   cin.sync_with_stdio(0);
//   cin.tie(0);
//   int subbb;
//   cin >> subbb;
//   int t;
//   cin >> t;
//   while (t--) {
//     int n;
//     cin >> n;
//     for (int i = 0; i < n; i++) {
//       for (int j = 0; j < n; j++) {
//         cin >> vals[i][j];
//       }
//     }

//     cout << hubDistance(n, 2) << '\n';
//   }
//   return 0;
// }

Compilation message (stderr)

towns.cpp: In function 'VI findHubs()':
towns.cpp:124:31: warning: conversion from '__gnu_cxx::__normal_iterator<int*, std::vector<int> >::difference_type' {aka 'long int'} to 'int' may change value [-Wconversion]
  124 |   ind1 = max_element(all(d0)) - d0.begin();
      |          ~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~
towns.cpp:126:31: warning: conversion from '__gnu_cxx::__normal_iterator<int*, std::vector<int> >::difference_type' {aka 'long int'} to 'int' may change value [-Wconversion]
  126 |   ind2 = max_element(all(d1)) - d1.begin();
      |          ~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~
towns.cpp:134:7: warning: declaration of 'mx' shadows a global declaration [-Wshadow]
  134 |   int mx = lcadists[ind2];
      |       ^~
towns.cpp:59:11: note: shadowed declaration is here
   59 | const int mx = 150;
      |           ^~
towns.cpp:137:8: warning: declaration of 'mp' shadows a global declaration [-Wshadow]
  137 |   MPII mp;
      |        ^~
towns.cpp:61:15: note: shadowed declaration is here
   61 | map<PII, int> mp;
      |               ^~
towns.cpp: In function 'int hubDistance(int, int)':
towns.cpp:224:31: warning: declaration of 'n' shadows a global declaration [-Wshadow]
  224 | int hubDistance(int n, int sub) {
      |                               ^
towns.cpp:58:5: note: shadowed declaration is here
   58 | int n;
      |     ^
towns.cpp:224:28: warning: unused parameter 'sub' [-Wunused-parameter]
  224 | int hubDistance(int n, int sub) {
      |                        ~~~~^~~
towns.cpp:227:23: warning: control reaches end of non-void function [-Wreturn-type]
  227 |   auto use = findHubs();
      |                       ^

• Subtask 1: 13 / 13
• Subtask 2: 12 / 12
• Subtask 3: 0 / 13
• Subtask 4: 0 / 10
• Subtask 5: 0 / 13
• Subtask 6: 0 / 39