結果

問題 No.3265 地元に帰れば天才扱い!
ユーザー edon8618
提出日時 2025-09-06 14:39:25
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,994 ms / 2,500 ms
コード長 4,323 bytes
コンパイル時間 6,657 ms
コンパイル使用メモリ 333,900 KB
実行使用メモリ 21,840 KB
最終ジャッジ日時 2025-09-06 14:40:38
合計ジャッジ時間 51,665 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;

// #include <boost/multiprecision/cpp_int.hpp>
// using namespace boost::multiprecision;

// # pragma GCC target("avx2")
// # pragma GCC optimize("O3")
// # pragma GCC optimize("unroll-loops")

#define ll long long
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
#define vi vector<int>
#define vl vector<ll>
#define vd vector<double>
#define vb vector<bool>
#define vs vector<string>
#define vc vector<char>
#define ull unsigned long long
#define chmax(a,b) a=max(a,b)
#define chmin(a,b) a=min(a,b)



ll inf=(1ll<<60);
// const double PI=3.1415926535897932384626433832795028841971;


// ll rui(ll a,ll b){
//     if(b==0)return 1;
//     if(b%2==1) return a*rui(a*a,b/2);
//     return rui(a*a,b/2);
// }

// ll kai(ll n){
//     if(n==0)return 1;
//     return n*kai(n-1);
// }

// using mint = modint998244353;//static_modint<998244353>
// using mint = modint1000000007;//static_modint<1000000007>
// using mint = modint;//mint::set_mod(mod);


// ll const mod=1000000007ll;
// ll const mod=998244353ll;
// ll modrui(ll a,ll b,ll mod){
//     a%=mod;
//     if(b==0)return 1;
//     if(b%2==1) return a*modrui(a*a%mod,b/2,mod)%mod;
//     return modrui(a*a%mod,b/2,mod)%mod;
// }

// ll inv(ll x){
//     x%=mod;
//     return modrui(x,mod-2);
// }

// ll modkai(ll n){
//     ll ret=1;
//     rep(i,n){
//         ret*=(i+1)%mod;
//         ret%=mod;
//     }
//     return ret;
// }


// void incr(vl &v,ll n){
//     ll k=v.size();
//     v[k-1]++;
//     ll now=k-1;
//     while (v[now]>=n)
//     {
//         v[now]=0;
//         if(now==0)break;
//         v[now-1]++;
//         now--;
//     }
//     return;
// }








ll op(ll a,ll b){
    return a+b;
}

ll e(){
    return 0;
}

//RAQ and RSQ
//range_add(l,r,x) a[l],a[l+1]...a[r-1]にxを加算
//get_sum(l,r) a[l],a[l+1]...a[r-1]の和を取得
struct SegTree {
    ll len;
    vector<ll> data, lazy;

    SegTree(ll sz) {
        len = 1;
        while (len < sz) len <<= 1;
        data.assign(2 * len, 0);
        lazy.assign(2 * len, 0);
    }

    void eval(ll k, ll sl, ll sr) {
        if (lazy[k] != 0) {
            data[k] += lazy[k] * (sr - sl);
            if (sr - sl > 1) {
                lazy[2 * k] += lazy[k];
                lazy[2 * k + 1] += lazy[k];
            }
            lazy[k] = 0;
        }
    }

    void range_add(ll l, ll r, ll x, ll sl = 0, ll sr = -1, ll k = 1) {
        if(sr == -1) sr = len;
        eval(k, sl, sr);
        if (r <= sl || sr <= l) return;
        if (l <= sl && sr <= r) {
            lazy[k] += x;
            eval(k, sl, sr);
        } else {
            ll sm = (sl + sr) / 2;
            range_add(l, r, x, sl, sm, 2 * k);
            range_add(l, r, x, sm, sr, 2 * k + 1);
            data[k] = data[2 * k] + data[2 * k + 1];
        }
    }

    ll get_sum(ll l, ll r, ll sl = 0, ll sr = -1, ll k = 1) {
        if(sr == -1) sr = len;
        eval(k, sl, sr);
        if (r <= sl || sr <= l) return 0;
        if (l <= sl && sr <= r) return data[k];
        ll sm = (sl + sr) / 2;
        return get_sum(l, r, sl, sm, 2 * k) + get_sum(l, r, sm, sr, 2 * k + 1);
    }
};

void solve(){
    ll n,m;
    cin >> n >> m;
    vl a(n),l(n),r(n);
    rep(i,n){
        cin >> a[i] >> l[i] >> r[i];
        l[i]--;
        r[i]--;
    }
    vl h(n);
    rep(i,n)h[i]=i;
    segtree<ll,op,e> seg(m);
    SegTree cnt(m);
    rep(i,n){
        seg.set(i,a[i]);
        cnt.range_add(l[i],r[i]+1,1);
    }
    ll ans=0;
    rep(i,n){
        ans+=a[i]*(r[i]-l[i]+1);
        ans-=seg.prod(l[i],r[i]+1);
    }

    ll q;
    cin >> q;
    while(q--){
        ll x,y,u,v;
        cin >> x >> y >> u >> v;
        x--;
        y--;
        u--;
        v--;
        ans+=seg.prod(l[x],r[x]+1);
        ans-=a[x]*(r[x]-l[x]+1);
        cnt.range_add(l[x],r[x]+1,-1);
        ans+=a[x]*(cnt.get_sum(h[x],h[x]+1));
        seg.set(h[x],0);

        h[x]=y;
        l[x]=u;
        r[x]=v;
        ans-=a[x]*cnt.get_sum(y,y+1);
        seg.set(y,a[x]);
        cnt.range_add(u,v+1,1);
        ans+=a[x]*(v-u+1);
        ans-=seg.prod(u,v+1);
        
        cout << ans << endl;
    }
}

int main(){
    // ios::sync_with_stdio(false);
    // std::cin.tie(nullptr);
    ll t=1;
    // cin >> t;
    while(t--)solve();
}
0