結果

問題 No.2988 Min-Plus Convolution Query
ユーザー Kude
提出日時 2024-12-13 12:03:34
言語 C++23
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

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