結果
問題 | No.1705 Mode of long array |
ユーザー |
|
提出日時 | 2021-10-21 21:46:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 289 ms / 3,000 ms |
コード長 | 5,758 bytes |
コンパイル時間 | 5,029 ms |
コンパイル使用メモリ | 255,464 KB |
最終ジャッジ日時 | 2025-01-25 02:36:40 |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 51 |
ソースコード
#pragma region Macros #include <bits/stdc++.h> #if defined(LOCAL) || defined(ONLINE_JUDGE) || defined(_DEBUG) #include <atcoder/all> #endif using namespace std; #define REP(i, n) for(int i=0, i##_len=(n); i<i##_len; ++i) #define REPR(i, n) for(int i=(n); i>=0; --i) #define FOR(i, n, m) for(int i=(m), i##_len=(n); i<i##_len; ++i) #define EACH(i, v) for(const auto& i : v) #define ALL(x) (x).begin(),(x).end() #define ALLR(x) (x).rbegin(),(x).rend() template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } template<class T>using vec = vector<T>; template<class T, class U>using umap = unordered_map<T, U>; template<class T>using uset = unordered_set<T>; using ll = long long; using ld = long double; using P = pair<ll, ll>; using T = tuple<ll, ll, ll>; using vl = vec<ll>; #define fi first #define se second #define el endl constexpr ll INF = numeric_limits<ll>::max()/2-1; #pragma endregion #pragma region IOMacros template<class T> istream &operator>>(istream &stream, vec<T>& o){REP(i, o.size())stream >> o[i];return stream;} template<class T> ostream &operator<<(ostream &stream, vec<T>& objs){REP(i, objs.size())stream << objs[i] << " ";stream << el;return stream;} #define I(T, ...) ;T __VA_ARGS__;__i(__VA_ARGS__); void __i() {} template<class T, class... Ts> void __i(T&& o, Ts&&... args){cin >> o;__i(forward<Ts>(args)...);} void O() {cout << el;} template<class T, class... Ts> void O(T&& o, Ts&&... args){cerr << o << " ";O(forward<Ts>(args)...);} #pragma endregion void Main(); int main(){ std::cin.tie(nullptr); std::cout << std::fixed << std::setprecision(15); Main(); return 0; } // 遅延セグメント木 #pragma region LazySegmentTree // https://ei1333.github.io/algorithm/segment-tree.html // [l, r) namespace LSTOpt{ using Monoid = P; using OpMonoid = ll; // Command基底クラスの宣言 struct CmdM{ virtual Monoid M() = 0; virtual Monoid op(Monoid a, Monoid b) = 0; virtual OpMonoid p(OpMonoid a, size_t n) = 0; }; struct CmdOp{ virtual OpMonoid ID() = 0; virtual Monoid mapping(Monoid a, OpMonoid b) = 0; virtual OpMonoid composition(OpMonoid a, OpMonoid b) = 0; }; } template< typename Monoid = LSTOpt::Monoid, typename OperatorMonoid = LSTOpt::OpMonoid > struct LST { int sz; vector< Monoid > data; vector< OperatorMonoid > lazy; LSTOpt::CmdM* m; LSTOpt::CmdOp* op; LST(int n, LSTOpt::CmdM* cmdm, LSTOpt::CmdOp* cmdop) //: f(cmdm.op), g(cmdop.mapping), h(cmdop.composition), p(cmdm.p), M1(cmdm.M()), OM0(cmdop.ID()) { this->m = cmdm; this->op = cmdop; sz = 1; while(sz < n) sz <<= 1; data.assign(2 * sz, m->M()); lazy.assign(2 * sz, op->ID()); } void set(int k, const Monoid &x) { data[k + sz] = x; } void build(){ for(int k = sz - 1; k > 0; k--) { data[k] = m->op(data[2 * k + 0], data[2 * k + 1]); } } void propagate(int k, int len) { if(lazy[k] != op->ID()) { if(k < sz) { lazy[2 * k + 0] = op->composition(lazy[2 * k + 0], lazy[k]); lazy[2 * k + 1] = op->composition(lazy[2 * k + 1], lazy[k]); } data[k] = op->mapping(data[k], m->p(lazy[k], len)); lazy[k] = op->ID(); } } Monoid update(int a, int b, const OperatorMonoid &x, int k, int l, int r) { propagate(k, r - l); if(r <= a || b <= l) { return data[k]; }else if(a <= l && r <= b) { lazy[k] = op->composition(lazy[k], x); propagate(k, r - l); return data[k]; }else{ return data[k] = m->op(update(a, b, x, 2 * k + 0, l, (l + r) >> 1), update(a, b, x, 2 * k + 1, (l + r) >> 1, r)); } } Monoid update(int a, int b, const OperatorMonoid &x) { return update(a, b, x, 1, 0, sz); } Monoid query(int a, int b, int k, int l, int r){ propagate(k, r - l); if(r <= a || b <= l) { return m->M(); }else if(a <= l && r <= b) { return data[k]; }else{ return m->op(query(a, b, 2 * k + 0, l, (l + r) >> 1), query(a, b, 2 * k + 1, (l + r) >> 1, r)); } } Monoid query(int a, int b) { return query(a, b, 1, 0, sz); } Monoid operator[](const int &k) { return query(k, k + 1); } }; namespace LSTOpt{ struct secadd : CmdOp{ // 単位元 OpMonoid ID(){ return 0; } // 要素と作用素のマージ Monoid mapping(Monoid a, OpMonoid b) { return {a.fi + b, a.se}; }; // 作用素と作用素のマージ OpMonoid composition(OpMonoid a, OpMonoid b) { return a + b; }; } ; // max struct secmax : CmdM{ // 単位元 Monoid M(){ return {-INF, -INF}; } // 演算 Monoid op(Monoid a, Monoid b){ return std::max(a, b); } // 作用素マージの前にかませる関数 OpMonoid p(OpMonoid a, size_t n){ return a; } }; }; #pragma endregion void Main(){ I(ll, N, M); LST lst(M, new LSTOpt::secmax(), new LSTOpt::secadd()); REP(i, M){ I(ll, a); lst.set(i, {a, i+1}); } lst.build(); I(ll, Q); REP(i, Q){ I(ll, T, X, Y); switch (T) { case 1: lst.update(X-1, X, Y);break; case 2: lst.update(X-1, X, -Y);break; case 3: cout << lst.query(0, M).se << el;break; } } }