結果
問題 | No.898 tri-βutree |
ユーザー |
|
提出日時 | 2019-10-04 21:43:45 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 329 ms / 4,000 ms |
コード長 | 1,760 bytes |
コンパイル時間 | 1,894 ms |
コンパイル使用メモリ | 180,716 KB |
実行使用メモリ | 22,436 KB |
最終ジャッジ日時 | 2024-11-08 21:57:46 |
合計ジャッジ時間 | 8,847 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 21 |
ソースコード
#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,int>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;}};int main(){int n;cin>>n;LCA lca(n);rep(i,n-1){int a,b,c;scanf("%d%d%d",&a,&b,&c);lca.add_edge(a,b,c);}int q;cin>>q;rep(i,q){int x,y,z;scanf("%d%d%d",&x,&y,&z);int p=lca.lca(x,lca.lca(y,z));ll ans=lca.dist(p,x)+lca.dist(p,y)+lca.dist(p,z);ans-=lca.dist(lca.lca(x,y),p)+lca.dist(lca.lca(y,z),p)+lca.dist(lca.lca(z,x),p);printf("%lld\n",ans);}}