#include using namespace std; using ull = unsigned long long; class SpecialSparseTable{ //LCA専用. private: int log = 0; pair op(pair &a,pair &b){return min(a,b);} pair add(pair &a,pair &b){return {a.first+b.first,a.second+b.second};} vector check,belong,Second; vector> Add; vector>>> Group; vector>> table; void maketable(vector> &A){ int n = A.size(); table.resize(n+1); check.resize(n+1); int p2 = 1; for(int i=1; i<=n; i++){ if(i == p2) check.at(i) = i,p2 *= 2; else check.at(i) = check.at(i-1); } table.at(1) = A; for(int i=2; i<=n; i*=2){ table.at(i).resize(n); for(int k=0; k<=n-i; k++) table[i][k] = op(table[i>>1][k],table[i>>1][k+(i>>1)]); } } void makeGroup(vector> &A){ int loop = 1<<(log-1); Group.resize(loop,vector>>(log,vector>(log))); for(int i=0; i now(1); for(int k=0; k now.at(r)) mina = now.at(r),pos = r; Group.at(i).at(l).at(r) = {mina,pos}; } } } } pair tablerange(int l,int r){//[L,R)だよ. int len = r-l,pos = check.at(len); return op(table.at(pos).at(l),table.at(pos).at(r-pos)); } public: void make(vector> &A){ int p2 = 1,n = A.size(); while(p2 < n) p2 *= 2,log++; log = max(1,(log+1)/2); vector> sepa; for(int i=0; i now = {1001001001,-1}; for(int k=i; k A.at(k)) now = {A.at(k).first,k}; } Add.push_back({A.at(i).first,i}); sepa.push_back(now); belong.push_back(g); } maketable(sepa); makeGroup(A); } int rangeans(int l,int r){ //[l,r)! int l2 = (l+log-1)/log,r2 = r/log; pair mind = {1001001001,-1}; if(l2 > r2) mind = add(Group[belong[r2]][l%log][(r-1)%log],Add[r2]); else{ if(l2 < r2) mind = tablerange(l2,r2); if(l%log) mind = min(mind,add(Group[belong[l2-1]][l%log][log-1],Add[l2-1])); if(r%log) mind = min(mind,add(Group[belong[r2]][0][(r-1)%log],Add[r2])); } return Second.at(mind.second); } }; class LCA{ private: vector dist,in,out; SpecialSparseTable ST; public: void make1(vector &P){// p0=-1でp1以降だけ int r = 0; vector> G(P.size()+1); for(int i=0; i &P){//pi=-1 iは固定されていない. int r = -1; vector> G(P.size()); for(int i=0; i> &Graph,int root){ //直接グラフを渡す. int t = 0,dep = 0,n = Graph.size(),pos = root; dist.resize(n),in.resize(n),out.resize(n); vector> depth; vector Itr(n,-1),P(n,-1); while(pos != -1){ depth.push_back({dep,pos}); if(++Itr.at(pos) == 0) in.at(pos) = t; int to = P.at(pos); if(Itr.at(pos) == Graph.at(pos).size()) out.at(pos) = ++t; else{ to = Graph.at(pos).at(Itr.at(pos)); if(to == P.at(pos)){ Itr.at(pos)++; if(Itr.at(pos) == Graph.at(pos).size()) out.at(pos) = ++t; else to = Graph.at(pos).at(Itr.at(pos)); } } if(to != P.at(pos)) dist.at(to) = dist.at(pos)+1,P.at(to) = pos,t++,dep++; else dep--; pos = to; } ST.make(depth); } int lca(int u,int v){ int tu = in.at(u),tv = in.at(v); if(tu > tv) swap(tu,tv); return ST.rangeans(tu,tv+1); } int twodist(int u,int v){return dist.at(u)+dist.at(v)-2*dist.at(lca(u,v));} int onedist(int u){return dist.at(u);} }; int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; vector> Graph(200000),Graph2(200000); for(int i=0; i> u >> v; u--; v--; Graph.at(u).push_back(v); Graph.at(v).push_back(u); } for(int i=0; i> u >> v; u--; v--; Graph2.at(u).push_back(v); Graph2.at(v).push_back(u); } int n = 200000; for(int i=N; i