結果
問題 | No.3116 More and more teleporter |
ユーザー |
|
提出日時 | 2025-04-19 02:51:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 281 ms / 2,000 ms |
コード長 | 4,042 bytes |
コンパイル時間 | 2,056 ms |
コンパイル使用メモリ | 199,208 KB |
実行使用メモリ | 7,848 KB |
最終ジャッジ日時 | 2025-04-19 02:51:15 |
合計ジャッジ時間 | 5,415 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 22 |
ソースコード
#include <bits/stdc++.h> //#include "atcoder/all" #pragma GCC target("avx2") #pragma GCC optimize("Ofast") using namespace std; using i64 = long long; template<typename T> class CHT { private: struct node { node *left, *right; static const T inf = numeric_limits<T>::max(); T a, b, l, r; node(const T _a, const T _b, const T _l, const T _r) : left(nullptr), right(nullptr), a(_a), b(_b), l(_l), r(_r){} T f(const T x) const { return (l <= x && x < r) ? (a * x + b) : inf; } T tf(const T x) const { return a * x + b; } bool isLine(const T _l, const T _r) const { return (l == _l) && (r == _r); } }; static void swap(node *x, node *y){ std::swap(x->a, y->a), std::swap(x->b, y->b), std::swap(x->l, y->l), std::swap(x->r, y->r); } void _add_segment(node *cur, node *nw, T l, T r){ while(true){ const T mid = (l + r) / 2; if(cur->f(nw->l) <= nw->tf(nw->l) && cur->f(nw->r - 1) <= nw->tf(nw->r - 1)){ break; } if(cur->tf(cur->l) >= nw->f(cur->l) && cur->tf(cur->r - 1) >= nw->f(cur->r - 1)){ swap(cur, nw); break; } if(cur->isLine(l, r) && nw->isLine(l, r)){ if(nw->tf(l) < cur->tf(l)) swap(cur, nw); if(cur->tf(mid) <= nw->tf(mid)){ if(!cur->right){ cur->right = new node(nw->a, nw->b, mid, r); break; }else{ cur = cur->right, l = mid, nw->l = mid; } }else{ swap(cur, nw); if(!cur->left){ cur->left = new node(nw->a, nw->b, l, mid); break; }else{ cur = cur->left, r = mid, nw->r = mid; } } continue; } if(nw->isLine(l, r)) swap(cur, nw); if(nw->r <= mid){ if(!cur->left){ cur->left = new node(*nw); break; }else{ cur = cur->left, r = mid; } }else if(mid <= nw->l){ if(!cur->right){ cur->right = new node(*nw); break; }else{ cur = cur->right, l = mid; } }else{ node _nw(nw->a, nw->b, mid, nw->r); nw->r = mid; if(!cur->left){ cur->left = new node(*nw); }else{ _add_segment(cur->left, nw, l, mid); } if(!cur->right){ cur->right = new node(_nw); }else{ _add_segment(cur->right, &_nw, mid, r); } break; } } } T query(node *cur, const T k, T l, T r) const { T ans = numeric_limits<T>::max(); while(cur){ ans = min(ans, cur->f(k)); const T mid = (l + r) / 2; if(k < mid){ cur = cur->left, r = mid; }else{ cur = cur->right, l = mid; } } return ans; } void clear(node *cur){ if(cur->left) clear(cur->left); if(cur->right) clear(cur->right); delete cur; } const T lpos, rpos; node *root; public: CHT(const T _lpos, const T _rpos) : lpos(_lpos), rpos(_rpos), root(new node(0, numeric_limits<T>::max(), lpos, rpos)){ assert(lpos < rpos); } // ~CHT(){ clear(root); } // f(x) = a * x + b を挿入 void add_line(const T a, const T b){ node nw(a, b, lpos, rpos); return _add_line(root, &nw, lpos, rpos); } // f(x) = a * x + b (x ∈ [l, r)) を挿入 void add_segment(const T a, const T b, const T l, const T r){ assert(l < r); node nw(a, b, l, r); return _add_segment(root, &nw, lpos, rpos); } // x = k での最小値 T query(const T k) const { return query(root, k, lpos, rpos); } }; void solve(){ int n, q; cin >> n >> q; CHT<i64> cht(0, n + 1); cht.add_segment(1, -1, 1, n + 1); for(int i = 0; i < q; ++i){ int ty; cin >> ty; if(ty == 1){ int x; cin >> x; cout << cht.query(x) << endl; } else{ int x; i64 c; cin >> x >> c; cht.add_segment(1, c - x, x, n + 1); cht.add_segment(-1, c + x, 0, x); } } } signed main(){ solve(); }