#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; class range {private: struct I{int x;int operator*(){return x;}bool operator!=(I& lhs){return x ostream& operator<<(ostream& os, const pair& p){ return os << "{" << p.first << ", " << p.second << "}"; } template ostream& operator<<(ostream& os, const vector& obj) { os << "{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const set& obj) { os << "set{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const map& obj) { os << "map{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template void take(vector& vec, int n) { vec.resize(n); for (int i = 0; i < n; ++i) cin >> vec[i]; } #ifdef ONLINE_JUDGE #define dump(expr) ; #else #define dump(expr) { cerr << "\033[33m#L" << __LINE__ << ": " << expr << "\033[39m" << endl; } #endif // range update, range query segment tree // based on https://maspypy.com/segment-tree-%e3%81%ae%e3%81%8a%e5%8b%89%e5%bc%b72 template struct LazySegmentTree { int n; vector xs; vector as; function fun; // (x @ y) fold function function op; // (x ^ a) operator function comp; // (a $ b) composition X unit_x; A unit_a; // Requirements: // (x @ y) @ z = x @ (y @ z) // (x^a) @ (y^a) = (x @ y)^a // (x^a)^b = x^(ab) // (a $ b) $ c = a $ (b $ c) LazySegmentTree( int n_, function fun_, function op_, function comp_, X unit_x_, A unit_a_ ) : n(n_), xs(2 * n, unit_x_), as(2 * n, unit_a_), fun(fun_), op(op_), comp(comp_), unit_x(unit_x_), unit_a(unit_a_) {} void build(vector seq) { assert((int)seq.size() == n); for (int i : range(n)) xs[n + i] = seq[i]; for (int i = n - 1; i >= 1; --i) xs[i] = fun(xs[i << 1], xs[(i << 1) | 1]); } void setval(int i, X x) { // xs[i] := x i += n; propagate_above(i); xs[i] = x; as[i] = unit_a; recalc_above(i); } X at(int i) { // returns xs[i] i += n; propagate_above(i); return eval_at(i); } X fold(int l, int r) { // xs[l] @ xs[l + 1] @ ... @ xs[r - 1] l += n; r += n; propagate_above(l / (l & -l)); propagate_above(r / (r & -r) - 1); X vl = unit_x; X vr = unit_x; while (l < r) { if (l & 1) { vl = fun(vl, eval_at(l)); l++; } if (r & 1) { r--; vr = fun(eval_at(r), vr); } l >>= 1; r >>= 1; } return fun(vl, vr); } void operate_range(int l, int r, A a) { // x[l ] := x[l ] ^ a // x[l + 1] := x[l + 1] ^ a // ... // x[r - 1] := x[l - 1] ^ a l += n; r += n; propagate_above(l / (l & -l)); propagate_above(r / (r & -r) - 1); while (l < r) { if (l & 1) { as[l] = comp(as[l], a); l++; } if (r & 1) { r--; as[r] = comp(as[r], a); } l >>= 1; r >>= 1; } recalc_above(l / (l & -l)); recalc_above(r / (r & -r) - 1); } //--------------------------------------- // auxiliary methods //--------------------------------------- X eval_at(int i) { return op(xs[i], as[i]); } void propagate_at(int i) { xs[i] = eval_at(i); as[i << 1] = comp(as[i << 1], as[i]); as[(i << 1) | 1] = comp(as[(i << 1) | 1], as[i]); as[i] = unit_a; } void propagate_above(int i) { const int max_bitlen = 25; for (int h = max_bitlen; h >= 1; --h) { if (i >> h) propagate_at(i >> h); } } void recalc_above(int i) { while (i > 1) { i >>= 1; xs[i] = fun(eval_at(i << 1), eval_at((i << 1) | 1)); } } }; template LazySegmentTree add_lazy_segment_tree_factory(int n) { return LazySegmentTree( n, [](ll x, ll y) -> T { return x + y; }, [](ll x, ll a) -> T { return x + a; }, [](ll a, ll b) -> T { return a + b; }, (T)0, (T)0 ); } namespace solver { template struct In2 { T1 a; T2 b; friend std::istream& operator>>(std::istream& is, In2& obj) { T1 t1; T2 t2; is >> t1 >> t2; obj = {t1, t2}; return is; } }; template struct In3 { T1 a; T2 b; T3 c; friend std::istream& operator>>(std::istream& is, In3& obj) { T1 t1; T2 t2; T3 t3; is >> t1 >> t2 >> t3; obj = {t1, t2, t3}; return is; } }; ll n, q; vector vs; vector> vq; void read() { cin >> n >> q; take(vs, n); take(vq, q); } using RetType = void; RetType run0() { using Pll = pair; LazySegmentTree seg( n, [](Pll, Pll) -> Pll { return {0, 0}; }, [](Pll p, ll a) -> Pll { return {p.first, a * p.first + p.second}; }, [](ll a, ll b) -> char { return a + b; }, Pll(0, 0), 0 ); vector init(n); for (int i : range(n)) init[i].first = vs[i]; seg.build(init); for (auto query : vq) { if (query.a == 'A') { // add ll i = query.b - 1; Pll p = seg.at(i); p.first += query.c; seg.setval(i, p); } else { // act ll i = query.b - 1; ll j = query.c - 1; seg.operate_range(i, j + 1, 1); } // dump(seg.xs); // dump(seg.as); } for (int i : range(n)) cout << seg.at(i).second << " "; cout << endl; } RetType run() { auto seg = add_lazy_segment_tree_factory(n); vector bs(n); for (int k = q - 1; k >= 0; --k) { auto query = vq[k]; if (query.a == 'A') { int i = query.b - 1; ll count = seg.fold(i, i + 1); bs[i] += count * query.c; } else { ll i = query.b - 1; ll j = query.c - 1; seg.operate_range(i, j + 1, 1); } } for (int i : range(n)) { ll count = seg.fold(i, i + 1); bs[i] += count * vs[i]; } dump(bs); for (ll b : bs) cout << b << " "; cout << endl; } } // namespace template void run(F f) { if constexpr (std::is_same_v) f(); else cout << f() << endl; } int main(int argc, char** argv) { cerr << fixed << setprecision(12); cout << fixed << setprecision(12); int testcase = 1; if (argc > 1) testcase = atoi(argv[1]); while (testcase--) { solver::read(); } run(solver::run); }