結果

問題 No.2988 Min-Plus Convolution Query
ユーザー KudeKude
提出日時 2024-12-13 12:03:34
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,589 ms / 2,000 ms
コード長 3,032 bytes
コンパイル時間 3,596 ms
コンパイル使用メモリ 280,512 KB
実行使用メモリ 98,604 KB
最終ジャッジ日時 2024-12-13 12:04:09
合計ジャッジ時間 32,944 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 1,550 ms
97,672 KB
testcase_03 AC 1,589 ms
96,932 KB
testcase_04 AC 1,507 ms
96,908 KB
testcase_05 AC 89 ms
25,924 KB
testcase_06 AC 2 ms
6,820 KB
testcase_07 AC 2 ms
6,816 KB
testcase_08 AC 2 ms
6,820 KB
testcase_09 AC 66 ms
15,644 KB
testcase_10 AC 77 ms
16,704 KB
testcase_11 AC 17 ms
6,816 KB
testcase_12 AC 27 ms
6,816 KB
testcase_13 AC 323 ms
25,240 KB
testcase_14 AC 36 ms
6,912 KB
testcase_15 AC 555 ms
41,308 KB
testcase_16 AC 1,413 ms
92,980 KB
testcase_17 AC 307 ms
35,280 KB
testcase_18 AC 289 ms
34,416 KB
testcase_19 AC 206 ms
31,448 KB
testcase_20 AC 1,220 ms
97,572 KB
testcase_21 AC 1,247 ms
97,148 KB
testcase_22 AC 1,414 ms
97,784 KB
testcase_23 AC 1,198 ms
97,272 KB
testcase_24 AC 1,193 ms
97,724 KB
testcase_25 AC 1,222 ms
97,152 KB
testcase_26 AC 1,175 ms
97,224 KB
testcase_27 AC 1,340 ms
97,284 KB
testcase_28 AC 1,217 ms
96,872 KB
testcase_29 AC 1,420 ms
96,912 KB
testcase_30 AC 1,217 ms
98,108 KB
testcase_31 AC 1,477 ms
98,604 KB
testcase_32 AC 1 ms
6,820 KB
testcase_33 AC 2 ms
6,820 KB
testcase_34 AC 2 ms
6,816 KB
testcase_35 AC 2 ms
6,820 KB
testcase_36 AC 2 ms
6,816 KB
testcase_37 AC 2 ms
6,816 KB
testcase_38 AC 2 ms
6,816 KB
testcase_39 AC 2 ms
6,816 KB
testcase_40 AC 2 ms
6,820 KB
testcase_41 AC 2 ms
6,820 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>;
constexpr int INF = 1001001001;

} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n, q;
  cin >> n >> q;
  VI a(n), b(n);
  rep(i, n) cin >> a[i];
  rep(i, n) cin >> b[i];
  VI t(n), k(q);
  int sz1 = bit_ceil(0U + q);
  vector<vector<P>> add_queries(2 * sz1);
  auto reg = [&](int i, int x, int l, int r) {
    l += sz1, r += sz1;
    while (l < r) {
      if (l & 1) add_queries[l++].emplace_back(i, x);
      if (r & 1) add_queries[--r].emplace_back(i, x);
      l >>= 1;
      r >>= 1;
    }
  };
  rep(i, q) {
    int p, x;
    cin >> p >> x >> k[i];
    p--, k[i] -= 2;
    reg(p, a[p], t[p], i);
    a[p] = x, t[p] = i;
  }
  rep(i, n) reg(i, a[i], t[i], q);
  int sz2 = bit_ceil(0U + 2 * n - 1);
  int log2 = countr_zero(0U + sz2);
  VI d(2 * sz2, -1);
  vector<P> hist;
  auto cmp = [&](int p, int q, int k) {
    int pb = k - p, qb = k - q;
    if (pb >= n) return false;
    if (qb < 0) return true;
    return a[p] + b[pb] < a[q] + b[qb];
  };
  auto get = [&](int k) {
    int k2 = sz2 + k;
    int res = 2 * INF;
    for (int i = log2; i >= 0; i--) {
      int p = d[k2 >> i];
      int pb = k - p;
      if (p != -1 && 0 <= pb && pb < n) chmin(res, a[p] + b[pb]);
    }
    return res;
  };
  auto dfs1 = [&](auto&& self, int u, int l, int r) -> void {
    int now = ssize(hist);
    for (auto [p, x] : add_queries[u]) {
      a[p] = x;
      auto dfs2  = [&](auto&& self, int u, int l, int r, int p) -> void {
        int q = d[u];
        hist.emplace_back(u, q);
        if (q == -1) {
          d[u] = p;
          return;
        }
        if (p > q) swap(p, q);
        if (l + 1 == r) {
          d[u] = cmp(p, q, l) ? p : q;
        } else {
          int c = (l + r) / 2;
          if (cmp(p, q, c)) {
            d[u] = p;
            self(self, 2 * u + 1, c, r, q);
          } else {
            d[u] = q;
            self(self, 2 * u, l, c, p);
          }
        }
      };
      dfs2(dfs2, 1, 0, sz2, p);
    }
    if (l + 1 == r) {
      int iq = u - sz1;
      if (iq < q) cout << get(k[iq]) << '\n';
    } else {
      int c = (l + r) / 2;
      self(self, 2 * u, l, c);
      self(self, 2 * u + 1, c, r);
    }
    while (ssize(hist) != now) {
      auto [v, p] = hist.back(); hist.pop_back();
      d[v] = p;
    }
  };
  dfs1(dfs1, 1, 0, sz1);
}
0