#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #ifdef CPP17 #include #endif #define endl codeforces #define ALL(v) std::begin(v), std::end(v) #define ALLR(v) std::rbegin(v), std::rend(v) using ll = std::int64_t; using ull = std::uint64_t; using pii = std::pair; using tii = std::tuple; using pll = std::pair; using tll = std::tuple; using size_type = ssize_t; template using vec = std::vector; template using vvec = vec>; template const T& var_min(const T &t) { return t; } template const T& var_max(const T &t) { return t; } template const T& var_min(const T &t, const Tail&... tail) { return std::min(t, var_min(tail...)); } template const T& var_max(const T &t, const Tail&... tail) { return std::max(t, var_max(tail...)); } template void chmin(T &t, const Tail&... tail) { t = var_min(t, tail...); } template void chmax(T &t, const Tail&... tail) { t = var_max(t, tail...); } template struct multi_dim_array { using type = std::array::type, Head>; }; template struct multi_dim_array { using type = std::array; }; template using mdarray = typename multi_dim_array::type; #ifdef CPP17 template void fill_seq(T &t, F f, Args... args) { if constexpr (std::is_invocable::value) { t = f(args...); } else { for (size_type i = 0; i < t.size(); i++) fill_seq(t[i], f, args..., i); } } #endif template vec make_v(size_type sz) { return vec(sz); } template auto make_v(size_type hs, Tail&&... ts) { auto v = std::move(make_v(std::forward(ts)...)); return vec(hs, v); } namespace init__ { struct InitIO { InitIO() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(30); } } init_io; } namespace utility { template using validate_integer = typename std::enable_if::value, ll>::type; template auto popcount(T n) -> validate_integer { return __builtin_popcount(n); } // 0 indexed template auto msb(T n) -> validate_integer { return 64 - __builtin_clzll(n) - 1; } template constexpr auto ceil_pow2(T s) -> validate_integer { ll ret = 1; while (ret < s) ret *= 2; return ret; } } namespace utility { struct has_id_ele { template auto operator ()(T &&t) -> decltype(T(), std::true_type()) { return std::true_type(); } std::false_type operator ()(...) { return std::false_type(); } }; struct has_merge { template auto operator ()(T &&t) -> decltype(T::merge(std::declval(), std::declval()), std::true_type()) { return std::true_type(); } std::false_type operator ()(...) { return std::false_type(); } }; struct has_apply { template auto operator ()(M &&m, Op &&op) -> decltype(m.apply(op), std::true_type()) { return std::true_type(); } std::false_type operator ()(...) { return std::false_type(); } }; template using callable = std::is_same::type, std::true_type>; template using is_monoid = std::conjunction< callable, callable>; template using enable_apply = callable; } namespace segtree { template class LazySegmentTree { static_assert(utility::is_monoid::value, "M must be monoid."); static_assert(utility::is_monoid::value, "Op must be monoid."); static_assert(utility::enable_apply::value, "Op is not operator of M."); struct segment { M m; Op op; bool has_lazy; segment(M m = M()) : m(m), op(Op()), has_lazy(false) { } void update_op(Op o) { m.apply(o); op = Op::merge(op, o); has_lazy = true; } void init_op() { op = Op(); has_lazy = false; } }; vec segs; size_type height; void push(size_type idx) { auto &s = segs[idx]; if (!s.has_lazy) return; for (int i = 0; i < 2; i++) { auto cidx = 2 * idx + i; if (segs.size() <= cidx) break; auto &cs = segs[cidx]; cs.update_op(s.op); } s.init_op(); } void propagate_from_top(size_type idx) { for (int i = height; 1 <= i; i--) push(idx >> i); } void update_from_bottom(size_type idx) { while (true) { auto pidx = idx / 2; if (pidx == 0) break; size_type c0 = 2 * pidx + 0, c1 = 2 * pidx + 1; push(c0); push(c1); segs[pidx].m = M::merge(segs[c0].m, segs[c1].m); idx = pidx; } } size_type get_endpoint_seg(size_type i) { i += size(); return i / (i & -i); } public: template LazySegmentTree(F f, size_type sz) { size_type sz2 = utility::ceil_pow2(sz); segs.resize(sz2 * 2); height = utility::msb(sz2); for (size_type i = 0; i < sz; i++) segs[i + sz2] = f(i); for (size_type i = sz2 - 1; 1 <= i; i--) segs[i] = M::merge(segs[2 * i].m, segs[2 * i + 1].m); } template LazySegmentTree(const vec &v) : LazySegmentTree([&](size_type i) { return v[i]; }, v.size()) { } size_type size() const { return segs.size() / 2; } template void update_query(size_type ql, size_type qr, const T &t) { Op op(t); auto l0 = get_endpoint_seg(ql); auto r0 = get_endpoint_seg(qr); propagate_from_top(l0); propagate_from_top(r0); size_type lnode = ql + size(), rnode = qr + size(); while (lnode < rnode) { if (lnode & 1) { segs[lnode].update_op(op); push(lnode); lnode++; } if (rnode & 1) { rnode--; segs[rnode].update_op(op); push(rnode); } lnode /= 2; rnode /= 2; } update_from_bottom(l0); update_from_bottom(r0); } M get_query(ll ql, ll qr) { auto ret = M(); auto l0 = get_endpoint_seg(ql); auto r0 = get_endpoint_seg(qr); propagate_from_top(l0); propagate_from_top(r0); size_type lnode = ql + size(), rnode = qr + size(); while (lnode < rnode) { if (lnode & 1) { push(lnode); ret = M::merge(segs[lnode].m, ret); lnode++; } if (rnode & 1) { rnode--; push(rnode); ret = M::merge(ret, segs[rnode].m); } lnode /= 2; rnode /= 2; } return ret; } }; } const ll inf = 5e15; struct Op { ll v; Op(ll v) : v(v) { } Op() : Op(0) { } static Op merge(Op a, Op b) { return Op(a.v + b.v); } }; struct M { ll sum, msum, l, r; M(ll sum, ll msum, ll l, ll r) : sum(sum), msum(msum), l(l), r(r) { } M() : M(0, 0, inf, -inf) { } static M merge(M a, M b) { ll sum = a.sum + b.sum; ll msum = a.sum * b.sum * 2 + a.msum + b.msum; ll l = std::min(a.l, b.l); ll r = std::max(a.r, b.r); return M(sum, msum, l, r); } void apply(Op o) { if (l == inf) return; ll cnt = r - l; sum += cnt * o.v; msum += o.v * o.v * (cnt - 1) * (cnt - 1); msum += sum * (cnt - 1) * o.v; } }; int main() { ll n; std::cin >> n; vec av(n); for (ll &e : av) std::cin >> e; ll q; std::cin >> q; segtree::LazySegmentTree segs([&](ll i) { return M(av[i], 0, i, i + 1); }, n); while (q--) { ll t, l, r; std::cin >> t >> l >> r; l--; if (t == 1) { ll x; std::cin >> x; segs.update_query(l, r, Op(x)); } else { auto m = segs.get_query(l, r); ll sum = m.sum, msum = m.msum; std::cout << sum * sum - msum << "\n"; } } return 0; }