結果
| 問題 |
No.1038 TreeAddQuery
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-10-09 20:16:59 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,703 bytes |
| コンパイル時間 | 2,320 ms |
| コンパイル使用メモリ | 193,360 KB |
| 実行使用メモリ | 242,444 KB |
| 最終ジャッジ日時 | 2024-07-20 07:54:26 |
| 合計ジャッジ時間 | 37,390 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 6 WA * 18 |
ソースコード
#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;
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);
}
};
vector<int>E[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;
vs[s].push_back(P(0,s));
for(int u:E[s]){
if(used[u])continue;
build(u,s,s,u,1);
auto&vs=vs2[s][u];
sort(vs.begin(),vs.end());
bit2[s][u]=BIT<ll>(vs.size());
rep(i,vs.size()){
Q[vs[i].second].push_back({s,u,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:E[s]){
if(used[u])continue;
solve(u,s);
}
}
void update(int x,int y,int z){
int prv=-1;
int d=0;
while(x!=-1){
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];
d++;
}
}
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);
}
return ans;
}
int main(){
int n,q;cin>>n>>q;
rep(i,n-1){
int a,b;scanf("%d%d",&a,&b);a--;b--;
E[a].push_back(b);
E[b].push_back(a);
}
solve(0,-1);
//~ 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);
}
}