結果

問題 No.1038 TreeAddQuery
ユーザー autumn-eel
提出日時 2020-10-09 22:56:29
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2,629 ms / 4,000 ms
コード長 4,394 bytes
コンパイル時間 2,332 ms
コンパイル使用メモリ 205,256 KB
実行使用メモリ 260,848 KB
最終ジャッジ日時 2024-07-20 13:37:06
合計ジャッジ時間 38,653 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

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

#include <bits/stdc++.h>
#define rep(i,n)for(int i=0;i<(n);i++)
using namespace std;
typedef long long ll;
typedef pair<ll,ll>P;
class LCA{
public:
    vector<vector<pair<int,int>>>E;
    vector<vector<int>>par;
    vector<int>d1;
    vector<long long>d2;
    int MAX_LOG=22;
    int n;
    LCA(){}
    LCA(int N){
        n=N;
        E=vector<vector<pair<int,int>>>(n);
        par=vector<vector<int>>(MAX_LOG,vector<int>(n));
        d1=vector<int>(n);
        d2=vector<long long>(n);
    }
    void add_edge(int a,int b,int c){
        E[a].push_back(make_pair(c,b));
        E[b].push_back(make_pair(c,a));
    }
    void add_edge(int a,int b){
        E[a].push_back(make_pair(1,b));
        E[b].push_back(make_pair(1,a));
    }
private:
    void dfs(int u,int p,int a,long long b){
        par[0][u]=p;d1[u]=a;d2[u]=b;
        for(auto&v:E[u]){
            if(v.second!=p)dfs(v.second,u,a+1,b+v.first);
        }
    }
    bool flag=false;
    void init(){
        dfs(0,-1,0,0);
        for(int i=1;i<MAX_LOG;i++)for(int j=0;j<n;j++){
            if(par[i-1][j]==-1)par[i][j]=-1;
            else par[i][j]=par[i-1][par[i-1][j]];
        }
        flag=true;
    }
public:
    int lca(int u,int v){
        if(!flag)init();
        if(d1[u]>d1[v])swap(u,v);
        for(int i=0;i<MAX_LOG;i++){
            if((d1[v]-d1[u])>>i&1)v=par[i][v];
        }
        if(u==v)return u;
        for(int i=MAX_LOG-1;i>=0;i--){
            if(par[i][u]!=par[i][v]){
                u=par[i][u];
                v=par[i][v];
            }
        }
        return par[0][u];
    }
    long long dist(int u,int v){
        int l=lca(u,v);
        return d2[u]+d2[v]-d2[l]*2;
    }
};
template<class T>
class BIT{
public:
    vector<T>bit;
    T cnt;
    BIT(){}
    BIT(int n){
        bit=vector<T>(n+10);
        cnt=0;
    }
    void add(int k,T x){
        k++;
        while(k<(int)bit.size()){
            bit[k]+=x;
            k+=k&-k;
        }
        cnt+=x;
    }
    T sum(int k){
        k++;
        T res=0;
        while(k){
            res+=bit[k];
            k-=k&-k;
        }
        return res;
    }
    T sum2(int k){
        return cnt-sum(k);
    }
};
LCA lca;
vector<int>E[200000],G[200000];
bool used[200000];
int N;
int centroid;
int sz[200000];
int par[200000];
vector<P>vs[200000];
unordered_map<int,vector<P>>vs2[200000];
BIT<ll>bit[200000];
unordered_map<int,BIT<ll>>bit2[200000];
struct st{
    int a,b,c;
};
vector<st>Q[200000];
void szdfs(int v,int p){
    sz[v]=1;
    for(int u:E[v]){
        if(u==p||used[u])continue;
        szdfs(u,v);
        sz[v]+=sz[u];
    }
}
void search_centroid(int v,int p){
    bool ok=true;
    for(int u:E[v]){
        if(u==p||used[u])continue;
        search_centroid(u,v);
        if(sz[u]>N/2)ok=false;
    }
    if(N-sz[v]>N/2)ok=false;
    if(ok)centroid=v;
}
void build(int v,int p,int s,int s2,int d){
    vs[s].push_back(P(d,v));
    vs2[s][s2].push_back(P(d,v));
    
    for(int u:E[v]){
        if(u==p||used[u])continue;
        build(u,v,s,s2,d+1);
    }
}
void solve(int v,int p){
    szdfs(v,-1);
    N=sz[v];
    search_centroid(v,-1);
    int s=centroid;
    used[s]=true;
    par[s]=p;
    if(p!=-1)G[p].push_back(s);
    for(int u:E[s]){
        if(used[u])continue;
        solve(u,s);
    }
}
void solve2(int s){
    used[s]=true;
    vs[s].push_back(P(0,s));
    int g=0;
    rep(i,E[s].size()){
        int u=E[s][i];
        if(used[u])continue;
        int t=G[s][g++];
        build(u,-1,s,t,1);
        auto&vs=vs2[s][t];
        sort(vs.begin(),vs.end());
        bit2[s][t]=BIT<ll>(vs.size());
        rep(i,vs.size()){
            Q[vs[i].second].push_back({s,t,i});
        }
    }
    
    sort(vs[s].begin(),vs[s].end());
    bit[s]=BIT<ll>(vs[s].size());
    rep(i,vs[s].size()){
        Q[vs[s][i].second].push_back({s,-1,i});
    }
    for(int u:G[s]){
        solve2(u);
    }
}
void update(int x,int y,int z){
    int X=x;
    int prv=-1;
    while(x!=-1){
        int d=lca.dist(X,x);
        int r=lower_bound(vs[x].begin(),vs[x].end(),P(y-d,INT_MAX))-vs[x].begin();
        bit[x].add(0,z);bit[x].add(r,-z);
        if(prv!=-1){
            auto&vs=vs2[x][prv];
            auto&bit=bit2[x][prv];
            int r2=lower_bound(vs.begin(),vs.end(),P(y-d,INT_MAX))-vs.begin();
            bit.add(0,z);bit.add(r2,-z);
        }
        prv=x;
        x=par[x];
    }
}
ll query(int x){
    ll ans=0;
    for(auto&s:Q[x]){
        if(s.b==-1)ans+=bit[s.a].sum(s.c);
        else ans-=bit2[s.a][s.b].sum(s.c);
        //~ cerr<<s.a<<' '<<s.b<<' '<<s.c<<' '<<(s.b==-1?bit[s.a].sum(s.c):bit2[s.a][s.b].sum(s.c))<<endl;
    }
    return ans;
}
int main(){
    int n,q;cin>>n>>q;
    lca=LCA(n);
    rep(i,n-1){
        int a,b;scanf("%d%d",&a,&b);a--;b--;
        E[a].push_back(b);
        E[b].push_back(a);
        lca.add_edge(a,b);
    }
    solve(0,-1);
    memset(used,false,sizeof(used));
    solve2(find(par,par+n,-1)-par);
    //~ rep(i,n){
        //~ cout<<par[i]<<' ';
        //~ for(auto j:vs[i])cout<<'('<<j.first<<','<<j.second<<") ";
        //~ cout<<endl;
    //~ }
    rep(i,q){
        int x,y,z;scanf("%d%d%d",&x,&y,&z);x--;
        printf("%lld\n",query(x));
        update(x,y,z);
    }
}
/*
6 3
1 5
2 3
2 4
3 5
5 6
1 5 14
5 1 23
2 2 29
*/
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0