結果

問題 No.1038 TreeAddQuery
ユーザー autumn-eelautumn-eel
提出日時 2020-10-10 14:40:02
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 2,198 ms / 4,000 ms
コード長 4,149 bytes
コンパイル時間 2,414 ms
コンパイル使用メモリ 199,572 KB
実行使用メモリ 229,724 KB
最終ジャッジ日時 2023-09-27 21:25:52
合計ジャッジ時間 33,952 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 14 ms
40,432 KB
testcase_01 AC 13 ms
40,436 KB
testcase_02 AC 12 ms
40,428 KB
testcase_03 AC 24 ms
42,884 KB
testcase_04 AC 27 ms
43,048 KB
testcase_05 AC 27 ms
43,368 KB
testcase_06 AC 23 ms
42,744 KB
testcase_07 AC 28 ms
43,200 KB
testcase_08 AC 1,456 ms
176,280 KB
testcase_09 AC 1,669 ms
186,700 KB
testcase_10 AC 1,735 ms
187,044 KB
testcase_11 AC 1,678 ms
186,688 KB
testcase_12 AC 1,740 ms
189,928 KB
testcase_13 AC 2,198 ms
229,508 KB
testcase_14 AC 2,124 ms
205,668 KB
testcase_15 AC 2,066 ms
201,856 KB
testcase_16 AC 2,012 ms
199,952 KB
testcase_17 AC 2,004 ms
197,780 KB
testcase_18 AC 377 ms
99,620 KB
testcase_19 AC 446 ms
107,740 KB
testcase_20 AC 470 ms
108,748 KB
testcase_21 AC 518 ms
115,568 KB
testcase_22 AC 693 ms
138,276 KB
testcase_23 AC 767 ms
139,464 KB
testcase_24 AC 1,234 ms
182,208 KB
testcase_25 AC 2,140 ms
229,724 KB
testcase_26 AC 1,290 ms
166,508 KB
権限があれば一括ダウンロードができます

ソースコード

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];
vector<vector<P>>vs2[200000];

BIT<ll>bit[200000];
vector<BIT<ll>>bit2[200000];

int id[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;
		vs2[s].push_back({});
		build(u,-1,s,g,1);
		auto&vs=vs2[s][g];
		id[G[s][g]]=g;
		sort(vs.begin(),vs.end());
		bit2[s].push_back(BIT<ll>(vs.size()));
		rep(i,vs.size()){
			Q[vs[i].second].push_back({s,g,i});
		}
		g++;
	}
	
	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=id[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);
	}
	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,q){
		int x,y,z;scanf("%d%d%d",&x,&y,&z);x--;
		printf("%lld\n",query(x));
		update(x,y,z);
	}
}
0