#line 2 "/Users/noya2/Desktop/Noya2_library/template/template.hpp" using namespace std; #include #line 1 "/Users/noya2/Desktop/Noya2_library/template/inout_old.hpp" namespace noya2 { template ostream &operator<<(ostream &os, const pair &p){ os << p.first << " " << p.second; return os; } template istream &operator>>(istream &is, pair &p){ is >> p.first >> p.second; return is; } template ostream &operator<<(ostream &os, const vector &v){ int s = (int)v.size(); for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i]; return os; } template istream &operator>>(istream &is, vector &v){ for (auto &x : v) is >> x; return is; } void in() {} template void in(T &t, U &...u){ cin >> t; in(u...); } void out() { cout << "\n"; } template void out(const T &t, const U &...u){ cout << t; if (sizeof...(u)) cout << sep; out(u...); } template void out(const vector> &vv){ int s = (int)vv.size(); for (int i = 0; i < s; i++) out(vv[i]); } struct IoSetup { IoSetup(){ cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); cerr << fixed << setprecision(7); } } iosetup_noya2; } // namespace noya2 #line 1 "/Users/noya2/Desktop/Noya2_library/template/const.hpp" namespace noya2{ const int iinf = 1'000'000'007; const long long linf = 2'000'000'000'000'000'000LL; const long long mod998 = 998244353; const long long mod107 = 1000000007; const long double pi = 3.14159265358979323; const vector dx = {0,1,0,-1,1,1,-1,-1}; const vector dy = {1,0,-1,0,1,-1,-1,1}; const string ALP = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"; const string alp = "abcdefghijklmnopqrstuvwxyz"; const string NUM = "0123456789"; void yes(){ cout << "Yes\n"; } void no(){ cout << "No\n"; } void YES(){ cout << "YES\n"; } void NO(){ cout << "NO\n"; } void yn(bool t){ t ? yes() : no(); } void YN(bool t){ t ? YES() : NO(); } } // namespace noya2 #line 2 "/Users/noya2/Desktop/Noya2_library/template/utils.hpp" #line 6 "/Users/noya2/Desktop/Noya2_library/template/utils.hpp" namespace noya2{ unsigned long long inner_binary_gcd(unsigned long long a, unsigned long long b){ if (a == 0 || b == 0) return a + b; int n = __builtin_ctzll(a); a >>= n; int m = __builtin_ctzll(b); b >>= m; while (a != b) { int mm = __builtin_ctzll(a - b); bool f = a > b; unsigned long long c = f ? a : b; b = f ? b : a; a = (c - b) >> mm; } return a << std::min(n, m); } template T gcd_fast(T a, T b){ return static_cast(inner_binary_gcd(std::abs(a),std::abs(b))); } long long sqrt_fast(long long n) { if (n <= 0) return 0; long long x = sqrt(n); while ((x + 1) * (x + 1) <= n) x++; while (x * x > n) x--; return x; } template T floor_div(const T n, const T d) { assert(d != 0); return n / d - static_cast((n ^ d) < 0 && n % d != 0); } template T ceil_div(const T n, const T d) { assert(d != 0); return n / d + static_cast((n ^ d) >= 0 && n % d != 0); } template void uniq(std::vector &v){ std::sort(v.begin(),v.end()); v.erase(unique(v.begin(),v.end()),v.end()); } template inline bool chmin(T &x, U y) { return (y < x) ? (x = y, true) : false; } template inline bool chmax(T &x, U y) { return (x < y) ? (x = y, true) : false; } template inline bool range(T l, T x, T r){ return l <= x && x < r; } } // namespace noya2 #line 8 "/Users/noya2/Desktop/Noya2_library/template/template.hpp" #define rep(i,n) for (int i = 0; i < (int)(n); i++) #define repp(i,m,n) for (int i = (m); i < (int)(n); i++) #define reb(i,n) for (int i = (int)(n-1); i >= 0; i--) #define all(v) (v).begin(),(v).end() using ll = long long; using ld = long double; using uint = unsigned int; using ull = unsigned long long; using pii = pair; using pll = pair; using pil = pair; using pli = pair; namespace noya2{ /* ~ (. _________ . /) */ } using namespace noya2; #line 2 "c.cpp" #line 2 "/Users/noya2/Desktop/Noya2_library/data_structure/range_tree.hpp" #line 2 "/Users/noya2/Desktop/Noya2_library/data_structure/compress.hpp" #line 5 "/Users/noya2/Desktop/Noya2_library/data_structure/compress.hpp" namespace noya2{ template struct compress { std::vector raws; compress () {} compress (const vector &_raws) : raws(_raws){ build(); } void build(){ std::sort(raws.begin(), raws.end()); raws.erase(std::unique(raws.begin(), raws.end()), raws.end()); } int id(const T &raw){ return lb(raw); } T raw(const int &id){ return raws[id]; } void add(const T &raw){ raws.emplace_back(raw); } void reserve(size_t sz){ raws.reserve(sz); } size_t size(){ return raws.size(); } int lb(const T &raw){ return lower_bound(all(raws),raw) - raws.begin(); } int ub(const T &raw){ return upper_bound(all(raws),raw) - raws.begin(); } bool contains(const T &raw){ int jd = lb(raw); if (jd < (int)size()) return raws[jd] == raw; return false; } int contains_id(const T &raw){ int jd = lb(raw); if (jd < (int)size() && raws[jd] == raw) return jd; return -1; } }; } // namespace noya2 #line 5 "/Users/noya2/Desktop/Noya2_library/data_structure/range_tree.hpp" namespace noya2 { template class data_structure, typename Idx = ll> struct range_tree { using DS = data_structure; using T = typename M::value_type; void join(Idx x, Idx y){ ps.emplace_back(x,y); } void build(){ for (auto &[x, y] : ps) xs.add(x); xs.build(); //siz = bit_ceil(xs.size()); siz = 1; while (siz < (int)(xs.size())) siz <<= 1; ys.resize(siz*2); for (auto &[x, y] : ps){ int xid = xs.id(x) + siz; ys[xid].add(y); while (xid > 1){ xid >>= 1; ys[xid].add(y); } } for (int i = 0; i < 2*siz; i++){ ys[i].build(); ds.emplace_back(ys[i].size()); } } void set(Idx p, Idx q, T val){ int i = xs.id(p) + siz; ds[i].set(ys[i].id(q),val); while (i > 1){ i >>= 1; T lr = M::e(); int i0 = ys[2*i+0].contains_id(q), i1 = ys[2*i+1].contains_id(q); if (i0 != -1) lr = M::op(lr, ds[2*i+0].get(i0)); if (i1 != -1) lr = M::op(lr, ds[2*i+1].get(i1)); ds[i].set(ys[i].id(q),lr); } } T get(Idx p, Idx q){ int ip = xs.contains_id(p); if (ip == -1) return M::e(); int i = ip + siz; int iq = ys[i].contains_id(q); if (iq == -1) return M::e(); return ds[i].get(iq); } T prod(Idx lp, Idx rp, Idx lq, Idx rq){ T res = M::e(); int li = xs.lb(lp), ri = xs.lb(rp); for (li += siz, ri += siz; li < ri; li >>= 1, ri >>= 1){ if (li & 1) res = M::op(res,ds[li].prod(ys[li].lb(lq),ys[li].lb(rq))), ++li; if (ri & 1) --ri, res = M::op(res,ds[ri].prod(ys[ri].lb(lq),ys[ri].lb(rq))); } return res; } int siz; vector> ps; compress xs; vector> ys; vector ds; }; } // namespace noya2 #line 2 "/Users/noya2/Desktop/Noya2_library/data_structure/binary_indexed_tree.hpp" #line 2 "/Users/noya2/Desktop/Noya2_library/misc/concepts.hpp" #include namespace noya2 { template concept Monoid = requires { typename monoid::value_type; {monoid::op(declval(),declval())} -> std::same_as; {monoid::e()} -> std::same_as; }; template concept Group = requires { requires Monoid; {group::inv(declval())} -> std::same_as; }; } // namespace noya2 #line 5 "/Users/noya2/Desktop/Noya2_library/data_structure/binary_indexed_tree.hpp" namespace noya2{ template struct binary_indexed_tree { using T = typename G::value_type; binary_indexed_tree (int _n = 0) : n(_n), d(_n + 1, G::e()) {} void add(int i, T val) { for (int x = i+1; x <= n; x += x & -x) { d[x] = G::op(d[x],val); } } T prod(int r){ return prefix_prod(r); } T prod(int l, int r) { return G::op(G::inv(prefix_prod(l)),prefix_prod(r)); } T get(int i){ return prod(i,i+1); } void set(int i, T val){ add(i,G::op(G::inv(get(i)),val)); } private: int n; std::vector d; T prefix_prod(int i) { assert(0 <= i && i <= n); T ret = G::e(); for (int x = i; x > 0; x -= x & -x) { ret = G::op(ret,d[x]); } return ret; } }; } // namespace noya2 #line 2 "/Users/noya2/Desktop/Noya2_library/misc/monoids.hpp" #line 4 "/Users/noya2/Desktop/Noya2_library/misc/monoids.hpp" namespace noya2{ template struct max_monoid { using value_type = T; static constexpr T op(const T &a, const T &b){ return max(a,b); } static constexpr T e(){ return std::numeric_limits::min(); } static constexpr T inv(const T &a){ return e(); } }; template struct min_monoid { using value_type = T; static constexpr T op(const T &a, const T &b){ return min(a,b); } static constexpr T e(){ return std::numeric_limits::max(); } static constexpr T inv(const T &a){ return e(); } }; template struct plus_group { using value_type = T; static constexpr T op(const T &a, const T &b){ return a + b; } static constexpr T e(){ return T(0); } static constexpr T inv(const T &a){ return -a; } }; template struct xor_group { using value_type = T; static constexpr T op(const T &a, const T &b){ return a ^ b; } static constexpr T e(){ return T(0); } static constexpr T inv(const T &a){ return a; } }; template struct monoid { using value_type = S; static constexpr value_type op(value_type a, value_type b){ return _op(a, b); } static constexpr value_type e(){ return _e(); } }; } // namespace noya2 #line 6 "c.cpp" void solve(){ int n, q; in(n,q); vector a(n); in(a); range_tree,binary_indexed_tree,int> rtsum; range_tree,binary_indexed_tree,int> rtcnt; vector> qs(q); rep(i,n){ rtsum.join(i,a[i]); rtcnt.join(i,a[i]); } rep(i,q){ int t; in(t); if (t == 1){ int x, y; in(x,y); x--; qs[i] = {x,y,-1,-1}; rtsum.join(x,y); rtcnt.join(x,y); } else { int l, r; in(l,r); l--; int ly, ry; in(ly,ry); qs[i] = {l,r,ly,ry}; } } rtsum.build(); rtcnt.build(); auto add = [&](int x, int y){ rtsum.set(x,y,y); rtcnt.set(x,y,1); }; auto del = [&](int x, int y){ rtsum.set(x,y,0); rtcnt.set(x,y,0); }; rep(i,n){ add(i,a[i]); } for (auto [l,r,ly,ry] : qs){ if (ly == -1){ int i = l, x = r; del(i,a[i]); a[i] = r; add(i,a[i]); } else { ll sum = rtsum.prod(l,r,ly,ry); ll ucnt = rtcnt.prod(l,r,0,ly); ll dcnt = rtcnt.prod(l,r,ry,1000000); ll ans = sum + ucnt * ly + dcnt * ry; out(ans); } } } int main(){ int t = 1; //in(t); while (t--) { solve(); } }