結果

問題 No.1038 TreeAddQuery
ユーザー autumn-eelautumn-eel
提出日時 2020-10-09 22:56:29
言語 C++14
(gcc 12.3.0 + boost 1.83.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 23 ms
52,172 KB
testcase_01 AC 22 ms
52,284 KB
testcase_02 AC 24 ms
52,284 KB
testcase_03 AC 40 ms
55,028 KB
testcase_04 AC 42 ms
55,124 KB
testcase_05 AC 40 ms
55,412 KB
testcase_06 AC 37 ms
54,840 KB
testcase_07 AC 42 ms
55,300 KB
testcase_08 AC 1,852 ms
201,612 KB
testcase_09 AC 2,045 ms
212,024 KB
testcase_10 AC 2,121 ms
211,780 KB
testcase_11 AC 2,088 ms
211,436 KB
testcase_12 AC 2,146 ms
214,556 KB
testcase_13 AC 2,629 ms
260,336 KB
testcase_14 AC 2,494 ms
231,296 KB
testcase_15 AC 2,418 ms
227,136 KB
testcase_16 AC 2,439 ms
225,136 KB
testcase_17 AC 2,330 ms
223,616 KB
testcase_18 AC 510 ms
118,860 KB
testcase_19 AC 621 ms
126,108 KB
testcase_20 AC 634 ms
127,360 KB
testcase_21 AC 688 ms
132,324 KB
testcase_22 AC 930 ms
153,556 KB
testcase_23 AC 1,001 ms
157,808 KB
testcase_24 AC 1,514 ms
203,352 KB
testcase_25 AC 2,466 ms
260,848 KB
testcase_26 AC 1,515 ms
190,824 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];
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
*/
0