結果

問題 No.1698 Face to Face
ユーザー noshi91noshi91
提出日時 2021-10-01 23:36:43
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 718 ms / 5,000 ms
コード長 6,053 bytes
コンパイル時間 1,469 ms
コンパイル使用メモリ 101,100 KB
実行使用メモリ 18,976 KB
最終ジャッジ日時 2023-09-26 22:51:13
合計ジャッジ時間 18,228 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,384 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 90 ms
17,532 KB
testcase_04 AC 354 ms
18,976 KB
testcase_05 AC 484 ms
17,956 KB
testcase_06 AC 449 ms
18,152 KB
testcase_07 AC 316 ms
18,316 KB
testcase_08 AC 343 ms
18,260 KB
testcase_09 AC 344 ms
18,356 KB
testcase_10 AC 275 ms
18,284 KB
testcase_11 AC 584 ms
18,108 KB
testcase_12 AC 375 ms
18,164 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 2 ms
4,380 KB
testcase_15 AC 718 ms
18,192 KB
testcase_16 AC 694 ms
17,780 KB
testcase_17 AC 715 ms
18,012 KB
testcase_18 AC 624 ms
16,668 KB
testcase_19 AC 598 ms
16,720 KB
testcase_20 AC 563 ms
16,216 KB
testcase_21 AC 256 ms
10,084 KB
testcase_22 AC 184 ms
8,800 KB
testcase_23 AC 36 ms
4,592 KB
testcase_24 AC 305 ms
10,636 KB
testcase_25 AC 657 ms
17,168 KB
testcase_26 AC 20 ms
4,380 KB
testcase_27 AC 118 ms
6,868 KB
testcase_28 AC 137 ms
7,108 KB
testcase_29 AC 91 ms
6,268 KB
testcase_30 AC 370 ms
11,556 KB
testcase_31 AC 259 ms
9,988 KB
testcase_32 AC 473 ms
14,856 KB
testcase_33 AC 475 ms
14,936 KB
testcase_34 AC 657 ms
17,228 KB
testcase_35 AC 2 ms
4,384 KB
testcase_36 AC 1 ms
4,384 KB
testcase_37 AC 2 ms
4,384 KB
testcase_38 AC 2 ms
4,384 KB
testcase_39 AC 2 ms
4,384 KB
testcase_40 AC 1 ms
4,380 KB
testcase_41 AC 1 ms
4,384 KB
testcase_42 AC 2 ms
4,380 KB
testcase_43 AC 2 ms
4,380 KB
testcase_44 AC 1 ms
4,384 KB
testcase_45 AC 453 ms
18,808 KB
testcase_46 AC 431 ms
18,636 KB
testcase_47 AC 173 ms
18,616 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#define NDEBUG
#pragma warning(disable : 4146)

#include <algorithm>
#include <cassert>
#include <cstddef>
#include <cstdint>
#include <iomanip>
#include <iostream>
#include <utility>
#include <vector>

namespace n91 {

using i32 = std::int32_t;
using i64 = std::int64_t;
using u32 = std::uint32_t;
using u64 = std::uint64_t;
using isize = std::ptrdiff_t;
using usize = std::size_t;
using f64 = double;

struct rep {
  struct itr {
    usize i;
    constexpr itr(const usize i) noexcept : i(i) {}
    void operator++() noexcept { ++i; }
    constexpr usize operator*() const noexcept { return i; }
    constexpr bool operator!=(const itr x) const noexcept { return i != x.i; }
  };
  const itr f, l;
  constexpr rep(const usize f, const usize l) noexcept
      : f(std::min(f, l)), l(l) {}
  constexpr auto begin() const noexcept { return f; }
  constexpr auto end() const noexcept { return l; }
};
struct revrep {
  struct itr {
    usize i;
    constexpr itr(const usize i) noexcept : i(i) {}
    void operator++() noexcept { --i; }
    constexpr usize operator*() const noexcept { return i; }
    constexpr bool operator!=(const itr x) const noexcept { return i != x.i; }
  };
  const itr f, l;
  constexpr revrep(const usize f, const usize l) noexcept
      : f(l - 1), l(std::min(f, l) - 1) {}
  constexpr auto begin() const noexcept { return f; }
  constexpr auto end() const noexcept { return l; }
};
template <class T> auto md_vec(const usize n, const T &value) {
  return std::vector<T>(n, value);
}
template <class... Args> auto md_vec(const usize n, Args... args) {
  return std::vector<decltype(md_vec(args...))>(n, md_vec(args...));
}
template <class T> constexpr T difference(const T &a, const T &b) noexcept {
  return a < b ? b - a : a - b;
}
template <class T> void chmin(T &a, const T &b) noexcept {
  if (b < a)
    a = b;
}
template <class T> void chmax(T &a, const T &b) noexcept {
  if (a < b)
    a = b;
}
template <class F> class rec_lambda {
  F f;

public:
  rec_lambda(F &&f_) : f(std::forward<F>(f_)) {}
  template <class... Args> auto operator()(Args &&... args) const {
    return f(*this, std::forward<Args>(args)...);
  }
};
template <class T> T scan() {
  T ret;
  std::cin >> ret;
  return ret;
}
constexpr char eoln = '\n';

i64 floor_div(const i64 n, const i64 d) {
  assert(d != 0);
  return n / d - static_cast<i64>((n ^ d) < 0 && n % d != 0);
}

i64 ceil_div(const i64 n, const i64 d) {
  assert(d != 0);
  return n / d + static_cast<i64>((n ^ d) >= 0 && n % d != 0);
}

#ifdef N91_LOCAL
#define OJ_LOCAL(a, b) b
#else
#define OJ_LOCAL(a, b) a
#endif

} // namespace n91

#include <cassert>
#include <cstdint>

namespace binary_search_impl {

using i64 = std::int64_t;

template <class F> i64 binary_search(F f, i64 left, i64 right) {
  assert(left <= right);

  while (left != right) {
    const i64 mid = left + (right - left) / 2;
    if (f(mid)) {
      left = mid + 1;
    } else {
      right = mid;
    }
  }
  return left;
}

} // namespace binary_search_impl

// first false
using binary_search_impl::binary_search;

#include <atcoder/dsu>
#include <atcoder/segtree>

namespace n91 {

std::vector<usize> inverse(const std::vector<usize> &p) {
  const usize n = p.size();
  std::vector<usize> q(n, -1);
  for (usize i = 0; i != n; i += 1) {
    assert(p[i] < n);
    assert(q[p[i]] == -1);
    q[p[i]] = i;
  }
  return q;
}

usize op(usize x, usize y) { return std::max(x, y); }

usize id() { return 0; }

void main_() {
  const usize n = scan<usize>();

  std::vector<usize> a(n), b(n), z(n);
  for (const usize i : rep(0, n)) {
    std::cin >> a[i];
    a[i] -= 1;
  }
  for (const usize i : rep(0, n)) {
    std::cin >> b[i];
    b[i] -= 1;
  }
  for (const usize i : rep(0, n)) {
    std::cin >> z[i];
    z[i] -= 1;
  }

  const auto ra = inverse(a), rb = inverse(b), rz = inverse(z);

  atcoder::segtree<usize, op, id> sega(a), segb(b);
  const auto get_left = [&](auto &seg, const usize i, const usize t) -> usize {
    if (seg.get(i) > t) {
      return i;
    } else {
      return seg.min_left(i, [&](const usize x) { return x <= t; });
    }
  };
  const auto get_right = [&](auto &seg, const usize i, const usize t) -> usize {
    if (seg.get(i) > t) {
      return i + 1;
    } else {
      return seg.max_right(i + 1, [&](const usize x) { return x <= t; });
    }
  };

  std::vector<std::vector<usize>> raft(n);
  {
    for (const usize i : rep(0, n)) {
      raft[binary_search(
               [&](const usize t) {
                 const usize l = std::max(get_left(sega, ra[i], t),
                                          get_left(segb, rb[z[i]], t));
                 const usize r = std::min(get_right(sega, ra[i], t),
                                          get_right(segb, rb[z[i]], t));
                 return !(l < r);
               },
               0, n - 1)]
          .push_back(i);
    }
  }

  usize ans = 0;

  std::vector<usize> count(n + 1, 4);
  std::vector<std::pair<usize, usize>> data(n, {1, 0});
  atcoder::dsu uf(n);

  const auto cost = [&](const usize i) {
    return std::min(data[i].first, data[i].second);
  };

  const auto dec = [&](const usize i) {
    count[i] -= 1;
    if (count[i] != 0) {
      return;
    }
    usize x = i - 1;
    usize y = i;
    x = uf.leader(x);
    y = uf.leader(y);
    ans -= cost(x);
    ans -= cost(y);
    uf.merge(x, y);
    if (uf.leader(x) != x) {
      std::swap(x, y);
    }
    data[x].first += data[y].first;
    data[x].second += data[y].second;
    ans += cost(x);
  };

  for (const usize t : rep(0, n)) {
    dec(ra[t]);
    dec(ra[t] + 1);
    dec(rb[t]);
    dec(rb[t] + 1);
    for (const usize i : raft[t]) {
      usize l = std::max(get_left(sega, ra[i], t), get_left(segb, rb[z[i]], t));
      l = uf.leader(l);
      ans -= cost(l);
      data[l].second += 1;
      ans += cost(l);
    }
    std::cout << ans << eoln;
  }
}

} // namespace n91

int main() {
  //*
  std::ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  //*/
  std::cout << std::setprecision(20);
  n91::main_();
  return 0;
}
0