結果

問題 No.2436 Min Diff Distance
ユーザー KudeKude
提出日時 2023-08-19 01:11:15
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,844 ms / 2,000 ms
コード長 3,610 bytes
コンパイル時間 3,268 ms
コンパイル使用メモリ 247,532 KB
実行使用メモリ 21,388 KB
最終ジャッジ日時 2024-11-28 13:33:06
合計ジャッジ時間 27,173 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 2 ms
6,816 KB
testcase_03 AC 1,830 ms
20,312 KB
testcase_04 AC 1,825 ms
20,196 KB
testcase_05 AC 1,833 ms
20,308 KB
testcase_06 AC 1,832 ms
20,184 KB
testcase_07 AC 1,842 ms
20,312 KB
testcase_08 AC 1,844 ms
20,312 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 2 ms
5,248 KB
testcase_11 AC 970 ms
20,644 KB
testcase_12 AC 967 ms
20,516 KB
testcase_13 AC 1,175 ms
18,092 KB
testcase_14 AC 156 ms
5,504 KB
testcase_15 AC 1,771 ms
21,388 KB
testcase_16 AC 775 ms
12,272 KB
testcase_17 AC 1,522 ms
20,112 KB
testcase_18 AC 1,377 ms
19,300 KB
testcase_19 AC 1,570 ms
20,396 KB
testcase_20 AC 1,036 ms
14,472 KB
testcase_21 AC 417 ms
7,860 KB
testcase_22 AC 159 ms
5,472 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic warning "-Wunused-function"
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;

template<class T>
vector<T> manhattan_dist_max(vector<pair<T, T>> ps) {
  T v1, v2, v3, v4;
  v1 = v3 = numeric_limits<T>::max();
  v2 = v4 = numeric_limits<T>::min();
  for(auto [x, y] : ps) {
    T s = x + y;
    T t = -x + y;
    chmin(v1, s);
    chmax(v2, s);
    chmin(v3, t);
    chmax(v4, t);
  }
  int n = ps.size();
  vector<T> res;
  res.reserve(n);
  for(auto [x, y] : ps) {
    T s = x + y;
    T t = -x + y;
    res.emplace_back(max({s - v1, v2 - s, t - v3, v4 - t}));
  }
  return res;
}

template<class T>
vector<T> manhattan_dist_min(const vector<pair<T, T>>& points) {
  const int n = points.size();
  assert(n >= 2);
  struct S {
    T x, y;
    int i;
  };
  vector<S> ps(n);
  for (int i = 0; i < n; i++) {
    auto [x, y] = points[i];
    ps[i] = {x + y, -x + y, i};
  }
  sort(ps.begin(), ps.end(), [](const S& p1, const S& p2) {
    return p1.x < p2.x;
  });
  vector<T> vals(n);
  for (int i = 0; i < n; i++) vals[i] = ps[i].y;
  sort(vals.begin(), vals.end());
  vals.erase(unique(vals.begin(), vals.end()), vals.end());
  const int sz = vals.size();
  vector<int> ys(n);
  for (int i = 0; i < n; i++) {
    ys[i] = lower_bound(vals.begin(), vals.end(), ps[i].y) - vals.begin();
  }

  vector<T> ls(n, -1), rs(n);
  auto dist = [&](const S& p1, const S& p2) {
    return max(abs(p1.x - p2.x), abs(p1.y - p2.y));
  };
  for (int i = 1; i < n; i++) rs[0] = rs[i] = dist(ps[0], ps[i]);
  struct E {
    T x;
    bool t;
    int i;
    bool operator<(const E& rhs) const {
      return x < rhs.x || (x == rhs.x && t < rhs.t);
    }
  };
  vector<E> evs;
  fenwick_tree<int> ft(sz);
  vector<int> cnt(n), yl(n), yr(n);
  while(true) {
    evs.clear();
    for (int i = 0; i < n; i++) if (rs[i] - ls[i] > 1) {
      T c = (rs[i] + ls[i]) / 2;
      assert(c >= 0);
      yl[i] = lower_bound(vals.begin(), vals.end(), ps[i].y - c) - vals.begin();
      yr[i] = upper_bound(vals.begin(), vals.end(), ps[i].y + c) - vals.begin();
      evs.push_back(E{ps[i].x - c, false, i});
      evs.push_back(E{ps[i].x + c, true, i});
    }
    if (evs.empty()) break;
    sort(evs.begin(), evs.end());
    int ptr = 0;
    for (auto [x, t, i] : evs) {
      while(ptr < n && ((!t && ps[ptr].x < x) || (t && ps[ptr].x <= x))) {
        ft.add(ys[ptr++], 1);
      }
      int sm = ft.sum(yl[i], yr[i]);
      if (!t) {
        cnt[i] = sm;
      } else {
        assert(cnt[i] + 1 <= sm);
        (cnt[i] + 1 == sm ? ls[i] : rs[i]) = (ls[i] + rs[i]) / 2;
      }
    }
  }
  vector<T> res(n);
  for (int i = 0; i < n; i++) res[ps[i].i] = rs[i];
  return res;
}


} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  vector<pair<int, int>> xy(n);
  for(auto& [x, y] : xy) cin >> x >> y;
  VI mx = manhattan_dist_max(xy);
  VI mn = manhattan_dist_min(xy);
  int ans = 2001001001;
  rep(i, n) chmin(ans, mx[i] - mn[i]);
  cout << ans << '\n';
}
0