結果

問題 No.1038 TreeAddQuery
ユーザー autumn-eelautumn-eel
提出日時 2020-10-09 22:01:17
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 4,270 bytes
コンパイル時間 2,462 ms
コンパイル使用メモリ 202,320 KB
実行使用メモリ 259,748 KB
最終ジャッジ日時 2023-09-27 17:01:51
合計ジャッジ時間 40,578 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 21 ms
52,088 KB
testcase_01 AC 20 ms
52,088 KB
testcase_02 AC 20 ms
52,120 KB
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 AC 525 ms
118,432 KB
testcase_19 AC 607 ms
126,596 KB
testcase_20 AC 631 ms
127,860 KB
testcase_21 AC 673 ms
132,816 KB
testcase_22 AC 955 ms
153,920 KB
testcase_23 AC 1,047 ms
158,260 KB
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
権限があれば一括ダウンロードができます

ソースコード

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));

	for(int u:G[s]){
		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: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);
	}
}
0