傳送門

題解：
好題。又是一種經典模型。

先假設沒有終止節點。
記fi$f_i$表示i$i$到fai$f{a}_i$的期望步數，gi$g_i$表示fai$f{a}_i$到i$i$的期望步數。

我們發現：每次詢問兩點，u$u$到lca$lca$的fi$f_i$,v$v$到lca$lca$的gi$g_i$不會失效（可以想想為什麼）。 然後處理出這兩個陣列就直接鏈查詢了。

處理就非常Naive了，按照以前的套路，可以推得：
fi=degi+∑j≠faifj$f_i=de{g}_i+\sum _{j\ne f{a}_i}{f}_j$
gi=degfai+∑j≠i,j≠fafaifj+gfai$g_i=de{g}_{f{a}_i}+\sum _{j\ne i,j\ne f{a}_{f{a}_i}}{f}_j+g_{f{a}_i}$。

O(n)$O(n)$DP即可。

#include <bits/stdc++.h>
 

using namespace std;
typedef long long LL;
typedef long double LD;

const int RLEN=1<<18|1;
inline char nc() {
    static char ibuf[RLEN],*ib,*ob;
    (ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
    char ch=nc(); int i=0,f=1
 
;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
    return i*f;
}
inline void W(LL x) {
    static int buf[50];
    if(!x) {putchar('0'); return;}
    if(x<0) {putchar('-'); x=-x;}
    while(x) {buf[++buf[0]]=x%10; x/=10;}
    while
 
(buf[0]) putchar(buf[buf[0]--]+'0');
}

const int N=1e5+50;

int n,src,des;
vector <int> edge[N];
int deg[N], sze[N], son[N], top[N], fa[N], dep[N];
LL f[N], g[N], sf[N], sg[N];
inline void dfs(int x,int ff) {
    sze[x]=1; fa[x]=ff; dep[x]=dep[ff]+1;
    if(x!=1 && deg[x]==1) return (void)(f[x]=1);
    LL sum=0;
    for(auto v:edge[x]) {
        if(v==ff) continue;
        dfs(v,x); sum+=f[v];
        sze[x]+=sze[v];
        son[x]=(sze[son[x]] < sze[v]) ? v : son[x];
    }
    f[x]=deg[x]+sum;
}
inline void dfs2(int x,int ff) {
    LL sum=0; sum+=g[x];
    for(auto v:edge[x]) if(v!=ff) sum+=f[v];
    for(auto v:edge[x]) if(v!=ff) g[v]=deg[x]+sum-f[v];
    for(auto v:edge[x]) if(v!=ff) dfs2(v,x);
}
inline void dfs3(int x,int ff) {
    sf[x]=sf[ff]+f[x]; sg[x]=sg[ff]+g[x];
    for(auto v:edge[x]) {
        if(v==ff) continue;
        top[v]=(v==son[x]) ? top[x] : v;
        dfs3(v,x);
    }
}
inline int lca(int x,int y) {
    while(top[x]!=top[y]) {
        (dep[top[x]] > dep[top[y]]) ? (x=fa[top[x]]) : (y=fa[top[y]]);
    } return (dep[x]>dep[y]) ? y : x;
}
int main() {
    n=rd(); 
    for(int i=1;i<n;i++) {
        int x=rd(), y=rd();
        edge[x].push_back(y);
        edge[y].push_back(x);
        ++deg[x]; ++deg[y];
    }
    dfs(1,0); dfs2(1,0); 
    top[1]=1; dfs3(1,0);
    for(int q=rd();q;q--) {
        int u=rd(), v=rd(), l=lca(u,v);
        W((sf[u]-sf[l]+sg[v]-sg[l]));
        putchar('\n');
    }
}