結果
問題 | No.1094 木登り / Climbing tree |
ユーザー |
|
提出日時 | 2022-05-28 00:15:02 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,420 ms / 2,000 ms |
コード長 | 2,580 bytes |
コンパイル時間 | 4,171 ms |
コンパイル使用メモリ | 240,296 KB |
実行使用メモリ | 71,364 KB |
最終ジャッジ日時 | 2024-09-20 17:41:23 |
合計ジャッジ時間 | 38,535 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 26 |
コンパイルメッセージ
main.cpp: In function 'int main()': main.cpp:51:15: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions] 51 | for(auto[u,c]:E[v])if(chmin(depth[u],depth[v]+1)){ | ^
ソースコード
#include<bits/stdc++.h> #include<atcoder/all> using namespace std; using namespace atcoder; typedef long long int ll; typedef long double ld; #define FOR(i,l,r) for(ll i=l;i<r;i++) #define REP(i,n) FOR(i,0,n) #define RFOR(i,l,r) for(ll i=r-1;i>=l;i--) #define RREP(i,n) RFOR(i,0,n) #define ALL(x) x.begin(),x.end() #define PA pair<ll,ll> #define F first #define S second #define BS(A,x) binary_search(ALL(A),x) #define LB(A,x) (ll)(lower_bound(ALL(A),x)-A.begin()) #define UB(A,x) (ll)(upper_bound(ALL(A),x)-A.begin()) #define COU(A,x) (UB(A,x)-LB(A,x)) #define sz(c) ((ll)(c).size()) //#include<boost/multiprecision/cpp_int.hpp> //namespace mp=boost::multiprecision; //using Bint=mp::cpp_int; template<typename T>using min_priority_queue=priority_queue<T,vector<T>,greater<T>>; template<typename T1,typename T2>ostream&operator<<(ostream&os,pair<T1,T2>&p){os<<p.F<<" "<<p.S;return os;} template<typename T1,typename T2>istream&operator>>(istream&is,pair<T1,T2>&p){is>>p.F>>p.S;return is;} template<typename T>ostream&operator<<(ostream&os,vector<T>&v){REP(i,sz(v))os<<v[i]<<(i+1!=sz(v)?" ":"");return os;} template<typename T>istream&operator>>(istream&is,vector<T>&v){for(T&in:v)is>>in;return is;} template<class T>bool chmax(T&a,T b){if(a<b){a=b;return 1;}return 0;} template<class T>bool chmin(T&a,T b){if(b<a){a=b;return 1;}return 0;} ld dist(ld x1,ld y1,ld x2,ld y2){return sqrtl((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));} //using mint=modint1000000007; using mint=modint998244353; //const ll mod=1000000007; const ll mod=998244353; //using mint=modint; ostream&operator<<(ostream&os,mint&a){os<<a.val();return os;} int main(){ ll N;cin>>N; vector<vector<PA>>E(N); REP(i,N-1){ ll a,b,c;cin>>a>>b>>c;a--;b--; E[a].emplace_back(PA(b,c)); E[b].emplace_back(PA(a,c)); } vector<ll>D(N),P(N,-1),depth(N,1e18); depth[0]=P[0]=0; queue<ll>Q; Q.emplace(0); while(sz(Q)){ ll v=Q.front();Q.pop(); for(auto[u,c]:E[v])if(chmin(depth[u],depth[v]+1)){ D[u]=D[v]+c;P[u]=v;Q.emplace(u); } } vector<vector<ll>>A(30,vector<ll>(N));A[0]=P; FOR(i,1,30)REP(v,N)A[i][v]=A[i-1][A[i-1][v]]; ll _;cin>>_; while(_--){ ll u,v;cin>>u>>v;u--;v--; if(depth[u]>depth[v])swap(u,v); ll w=v; RREP(i,30)if(depth[u]<=depth[w]-(1<<i))w=A[i][w]; ll l=0,r=1e9; while(l!=r){ ll m=(l+r)/2; ll _w=w,_u=u; REP(i,30)if((m>>i)%2)_w=A[i][_w],_u=A[i][_u]; if(_w==_u)r=m;else l=m+1; } REP(i,30)if((l>>i)%2)w=A[i][w]; cout<<D[u]+D[v]-2*D[w]<<endl; } return 0; }