#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; } }; 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: vec segs; 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; } }; } using i128 = __int128_t; std::istream& operator >>(std::istream &is, i128 &n) { std::string s; std::cin >> s; n = 0; bool minus = false; for (char c : s) { if (c == '-') { minus = true; continue; } n = n * 10 + (c - '0'); } if (minus) n *= -1; return is; } std::ostream& operator <<(std::ostream &os, i128 n) { char buf[50] = { }; int idx = 0; bool minus = false; if (n < 0) { minus = true; buf[idx++] = '-'; n *= -1; } for (; n; n /= 10, idx++) buf[idx] = static_cast('0' + (n % 10)); std::reverse(buf + minus, buf + idx); std::cout << buf; return os; } const i128 inf = 5e15; struct Op { i128 v; Op(i128 v) : v(v) { } Op() : Op(0) { } static Op merge(Op a, Op b) { return Op(a.v + b.v); } }; struct M { i128 sum, msum, l, r; M(i128 sum, i128 msum, i128 l, i128 r) : sum(sum), msum(msum), l(l), r(r) { } M() : M(0, 0, inf, -inf) { } static M merge(M a, M b) { i128 sum = a.sum + b.sum; i128 msum = a.sum * b.sum * 2 + a.msum + b.msum; i128 l = std::min(a.l, b.l); i128 r = std::max(a.r, b.r); return M(sum, msum, l, r); } void apply(Op o) { if (l == inf) return; i128 cnt = r - l; msum += o.v * o.v * cnt * (cnt - 1); msum += sum * (cnt - 1) * 2 * o.v; sum += cnt * o.v; } }; int main() { i128 n; std::cin >> n; vec av(n); for (i128 &e : av) std::cin >> e; i128 q; std::cin >> q; segtree::LazySegmentTree segs([&](i128 i) { return M(av[i], 0, i, i + 1); }, n); while (q--) { i128 t, l, r; std::cin >> t >> l >> r; l--; if (t == 1) { i128 x; std::cin >> x; segs.update_query(l, r, Op(x)); } else { auto m = segs.get_query(l, r); i128 sum = m.sum, msum = m.msum; std::cout << sum * sum - msum << "\n"; } #if 0 { for(i128 i=0;i