結果

問題 No.1074 増殖
ユーザー 沙耶花
提出日時 2020-06-05 22:23:30
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 415 ms / 2,000 ms
コード長 3,837 bytes
コンパイル時間 2,166 ms
コンパイル使用メモリ 199,504 KB
最終ジャッジ日時 2025-01-10 22:43:00
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:171:10: warning: use of ‘auto’ in parameter declaration only available with ‘-std=c++20’ or ‘-fconcepts’
  171 | int calc(auto &seg,int maxi,int value){
      |          ^~~~

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
#define modulo 1000000007
#define mod(mod_x) ((((long long)mod_x+modulo))%modulo)
#define Inf 1000000000000000010
template <typename T0,typename T1,typename F0,typename F1,typename F2>
struct lazysegtree{
    //元データx[i]はv[n+i]
    //v[i]の親はv[i/2],子はv[i*2]とv[i*2+1]
    F0 func0;
    F1 func1;
    F2 func2;
    vector<T0> v0;
    vector<T1> v1;
    int n;
    int cnt;
    
    T0 init_value0;
    T1 init_value1;
    
    lazysegtree(int sz,F0 f0,F1 f1,F2 f2,T0 iv0,T1 iv1):func0(f0),func1(f1),func2(f2){
        init_value0 = iv0;
        init_value1 = iv1;
        n=1;
        cnt=0;
        while(true){
            if(n>=sz)break;
            n*=2;
            cnt++;
        }
        v0.resize(2*n,init_value0);
        v1.resize(2*n,init_value1);
        
    }
    
    lazysegtree(vector<T0> &x,F0 f0,F1 f1,F2 f2,T0 iv0,T1 iv1):func0(f0),func1(f1),func2(f2){
        init_value0 = iv0;
        init_value1 = iv1;
        n=1;
        cnt=0;
        while(true){
            if(n>=x.size())break;
            n*=2;
            cnt++;
        }
        v0.resize(2*n,init_value0);
        v1.resize(2*n,init_value1);
        for(int i=0;i<x.size();i++){
            v0[n+i]=x[i];
        }
        for(int i=n-1;i>=0;i--){
            v0[i]=func0(v0[i<<1],v0[(i<<1)+1]);
        }   
    }
    //2人の子供に伝える
    void propagate(int ind){
        update(ind<<1,v1[ind]);
        update((ind<<1)+1,v1[ind]);
        v1[ind] = init_value1;
    }
    
    //あるノードに対し先祖から伝播  
    void reflect(int l,int r){
        int j = cnt;
        while(l>>j==r>>j&&j>=1){
            propagate(l>>j);
            j--;
        }
        for(;j>=1;j--){
            propagate(l>>j);
            propagate(r>>j);
        }
    }
    
    //子供の値を親に伝える
    void mergeChildren(int ind){
        v0[ind] = func1(func0(v0[ind<<1],v0[(ind<<1)+1]),v1[ind],n>>(31-__builtin_clz(ind)));
    }
    
    //ある要素について作用させる
    void update(int ind,T1 x){
        v0[ind] = func1(v0[ind],x,n>>(31-__builtin_clz(ind)));
        v1[ind] = func2(v1[ind],x);
    }
    
    //[l,r)に対して作用させる
    void update(int l,int r,T1 x){
        if(l>=r)return;
        int L = l,R = r;
        l+=n;
        r+=n;
        reflect(l,r-1);
        while(true){
            if(l&1){
                update(l,x);
                l++;
            }
            if(r&1){
                update(r-1,x);
                r--;
            }
            if(l>=r)break;
            l>>=1;
            r>>=1;
        }
        
        l = L + n;
        r = R + n-1;
        
        while(true){
            l>>=1;
            r>>=1;
            if(l<=0)break;
            if(l==r){
                while(true){
                    mergeChildren(l);
                    l>>=1;
                    if(!l)return;
                }
            }
            else{
                mergeChildren(l);
                mergeChildren(r);
            }
        }
        
    }
    
    //区間[l,r)におけるクエリ処理
    T0 query(int l,int r){
        T0 res1 = init_value0;
        T0 res2 = init_value0;
        if(l>=r)return res1;
        l+=n;
        r+=n;
        reflect(l,r-1);
        while(true){
            if(l&1){
                res1=func0(res1,v0[l]);
                l++;
            }
            if(r&1){
                res2=func0(v0[r-1],res2);
                r--;
            }
            if(l>=r)break;
            l>>=1;
            r>>=1;
        }
        return func0(res1,res2);
    }
    void show(){
        int n = 1;
        for(int i=1;i<v0.size();i++){
            for(int j=0;j<n;j++){
                if(j!=0)cout<<' ';
                cout<<v0[i+j];
            }
            cout<<endl;
            i+=n-1;
            n*=2;
        }
    }
    
};
int calc(auto &seg,int maxi,int value){
    int ok = maxi+1;
    int ng = -1;
    while(ok-ng>1){
        int mid = (ok+ng)/2;
        int t = seg.query(mid,mid+1);
        if(t>=value)ng = mid;
        else ok = mid;
    }
    //cout<<ok<<endl;
    int ret =  ((maxi+1)-ok)*value-seg.query(ok,maxi+1);
    seg.update(ok,maxi+1,value);
    return ret;
}
int main(){
    
    auto f0 = [](long long a,long long b){
        return a+b;
    };
    
    auto f1 = [](long long a,long long b,long long sz){
        if(b==-1)return a;
        return b*sz;
    };
    
    auto f2 = [](long long a,long long b){
        if(b==-1)return a;
        return b;
    };
    
    lazysegtree<long long,long long,decltype(f0),decltype(f1),decltype(f2)> rp(20000,f0,f1,f2,0,-1);
    auto rm=rp,lp=rp,lm=rp;
    
    int N;
    cin>>N;
    
    for(int i=0;i<N;i++){
        int Xa,Ya,Xb,Yb;
        cin>>Xa>>Ya>>Xb>>Yb;
        
        Xa++,Ya++,Xb--,Yb--;
        
        int ans = 0;
        ans += calc(rp,Yb,Xb+1);
        ans += calc(lp,Yb,-Xa+1);
        ans += calc(rm,-Ya,Xb+1);
        ans += calc(lm,-Ya,-Xa+1);
        cout<<ans<<endl;
    }
    
    
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0