結果

問題 No.1099 Range Square Sum
ユーザー IKyoproIKyopro
提出日時 2020-06-26 22:00:51
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 153 ms / 2,000 ms
コード長 3,216 bytes
コンパイル時間 1,974 ms
コンパイル使用メモリ 207,424 KB
実行使用メモリ 17,664 KB
最終ジャッジ日時 2024-07-04 21:08:39
合計ジャッジ時間 4,701 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 3 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 3 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 3 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 124 ms
17,524 KB
testcase_22 AC 153 ms
17,664 KB
testcase_23 AC 128 ms
17,536 KB
testcase_24 AC 126 ms
17,536 KB
testcase_25 AC 128 ms
17,408 KB
testcase_26 AC 103 ms
17,536 KB
testcase_27 AC 109 ms
17,664 KB
testcase_28 AC 102 ms
17,536 KB
testcase_29 AC 100 ms
17,536 KB
testcase_30 AC 105 ms
17,536 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
template <class T, class U> using Pa = pair<T, U>;
template <class T> using vec = vector<T>;
template <class T> using vvec = vector<vec<T>>;

template<typename Monoid,typename OperatorMonoid,typename F,typename G,typename H>
class LazySegmentTree {
private:
    int sz,height;
    vec<Monoid> data;
    vec<OperatorMonoid> lazy;
    const F op;
    const G homo;
    const H comp;
    const Monoid e;
    const OperatorMonoid Oe;
public:
    LazySegmentTree(int n,const F op,const G homo,const H comp,
                    const Monoid &e,const OperatorMonoid Oe)
        : op(op),homo(homo),comp(comp),e(e),Oe(Oe) {
        sz = 1;
        height = 0;
        while(sz<=n) sz <<= 1,height++;
        data.assign(2*sz,e);
        lazy.assign(2*sz,Oe);
    }

    void set(int k,const Monoid &x) {
        data[k+sz] = x;
    }

    void build() {
        for(int k=sz-1;k>0;k--) {
            data[k] = op(data[2*k], data[2*k+1]);
        }
    }

    inline void propagate(int k) {
        if(lazy[k]!=Oe) {
            lazy[2*k] = comp(lazy[2*k], lazy[k]);
            lazy[2*k+1] = comp(lazy[2*k+1], lazy[k]);
            data[k] = reflect(k);
            lazy[k] = Oe;
        }
    }

    inline Monoid reflect(int k) {
        return lazy[k] == Oe? data[k]:homo(data[k],lazy[k]);
    }

    inline void recalc(int k) {
        while(k>>=1) data[k] = op(reflect(2*k), reflect(2*k+1));
    }

    inline void thrust(int k) {
        for(int i=height;i>0;i--) propagate(k>>i);
    }

    void update(int a, int b, const OperatorMonoid &x) {
        thrust(a+=sz);
        thrust(b+=sz-1);
        for(int l=a,r=b+1;l<r;l>>=1,r>>=1) {
            if(l&1) lazy[l] = comp(lazy[l],x),++l;
            if(r&1) --r, lazy[r] = comp(lazy[r],x);
        }
        recalc(a);
        recalc(b);
    }

    Monoid query(int a, int b) {
        thrust(a+=sz);
        thrust(b+=sz-1);
        Monoid L = e, R = e;
        for(int l=a, r=b+1;l<r;l>>= 1,r>>=1) {
            if(l&1) L = op(L,reflect(l++));
            if(r&1) R = op(reflect(--r),R);
        }
        return op(L,R);
    }

    Monoid operator[](const int &k) {
        return query(k,k+1);
    }
};

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);
    struct state{
        ll sum,sum2,len;
    };

    auto op = [](state L,state R)->state{
        ll sum = L.sum+R.sum;
        ll sum2 = L.sum2+R.sum2;
        ll len = L.len+R.len;
        return {sum,sum2,len};
    };

    auto func = [](state S,ll x)->state{
        ll sum = S.sum+x*S.len;
        ll sum2 = S.sum2+2*S.sum*x+x*x*S.len;
        return {sum,sum2,S.len};
    };

    auto comp = [](ll x,ll y){return x+y;};

    int N;
    cin >> N;
    state e = {0,0,0};
    LazySegmentTree seg(N,op,func,comp,e,0);
    
    for(int i=0;i<N;i++){
        ll a;
        cin >> a;
        seg.set(i,(state){a,a*a,1});
    }
    seg.build();
    int Q;
    cin >> Q;
    while(Q--){
        int t,l,r;
        cin >> t >> l >> r;
        l--,r--;
        if(t==1){
            ll x;
            cin >> x;
            seg.update(l,r+1,x);
        }else{
            cout << seg.query(l,r+1).sum2 << "\n";
        }
    }
    
}
0