假设询问两点中d[x]<d[y]。
若x是y的祖先，那么就是求起点不在x到y方向的第一个点的子树中，且终点在y子树中的跳跃点个数。
若x不是y的祖先，那么就是求起点在x子树中，且终点在y子树中的跳跃点个数。
利用DFS序可以将转化成矩形，于是利用CDQ分治+扫描线+树状数组即可做到$O(nlog^2n)$。
#include<cstdio>
#include<algorithm>
#define N 200010
using namespace std;
int n,m,q,i,op,x,y,T,pos[N],bit[N],ans[N],cb,cc;
int g[N],v[N],nxt[N],ed,f[N],size[N],son[N],top[N],d[N],st[N],en[N],dfn;
struct P{int x,y,t;P(){}P(int _x,int _y,int _t){x=_x,y=_y,t=_t;}}a[N],b[N];
inline bool cmpb(const P&a,const P&b){return a.x<b.x;}
struct C{int x,l,r,t,p;C(){}C(int _x,int _l,int _r,int _t,int _p){x=_x,l=_l,r=_r,t=_t,p=_p;}}c[N];
inline bool cmpc(const C&a,const C&b){return a.x<b.x;}
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
void dfs(int x){
size[x]=1;
for(int i=g[x];i;i=nxt[i])if(v[i]!=f[x]){
d[v[i]]=d[f[v[i]]=x]+1;dfs(v[i]);size[x]+=size[v[i]];
if(size[v[i]]>size[son[x]])son[x]=v[i];
}
}
void dfs2(int x,int y){
st[x]=++dfn;top[x]=y;
if(son[x])dfs2(son[x],y);
for(int i=g[x];i;i=nxt[i])if(v[i]!=son[x]&&v[i]!=f[x])dfs2(v[i],v[i]);
en[x]=dfn;
}
inline int lca2(int x,int y){
int t;
while(top[x]!=top[y])t=top[y],y=f[top[y]];
return x==y?t:son[x];
}
inline void ins(int x,int y){for(;x<=n;x+=x&-x)if(pos[x]<T)pos[x]=T,bit[x]=y;else bit[x]+=y;}
inline int ask(int x){int t=0;for(;x;x-=x&-x)if(pos[x]==T)t+=bit[x];return t;}
void solve(int l,int r){
if(l==r)return;
int mid=(l+r)>>1;
solve(l,mid),solve(mid+1,r);
int i,j;
cb=cc=0;
for(i=l;i<=mid;i++){
if(a[i].t<0){
b[cb++]=P(a[i].x,a[i].y,-1);
b[cb++]=P(a[i].y,a[i].x,-1);
}
if(!a[i].t){
b[cb++]=P(a[i].x,a[i].y,1);
b[cb++]=P(a[i].y,a[i].x,1);
}
}
if(!cb)return;
for(i=r;i>mid;i--){
if(a[i].t==1){
c[cc++]=C(n,st[a[i].y],en[a[i].y],1,i);
c[cc++]=C(st[a[i].x]-1,st[a[i].y],en[a[i].y],1,i);
c[cc++]=C(en[a[i].x],st[a[i].y],en[a[i].y],-1,i);
}
if(a[i].t==2){
c[cc++]=C(st[a[i].x]-1,st[a[i].y],en[a[i].y],-1,i);
c[cc++]=C(en[a[i].x],st[a[i].y],en[a[i].y],1,i);
}
}
if(!cc)return;
if(cb>1)sort(b,b+cb,cmpb);
if(cc>1)sort(c,c+cc,cmpc);
for(T++,i=j=0;i<cc;i++){
while(j<cb&&b[j].x<=c[i].x)ins(b[j].y,b[j].t),j++;
ans[c[i].p]+=c[i].t*(ask(c[i].r)-ask(c[i].l-1));
}
}
int main(){
read(n);
for(i=1;i<n;i++)read(x),read(y),add(x,y),add(y,x);
dfs(1),dfs2(1,1);
read(q);
while(q--)read(x),read(y),a[++m]=P(st[x],st[y],0);
read(q);
while(q--){
read(op),read(x),read(y);
if(op==1)a[++m]=P(st[x],st[y],0);
if(op==2)a[++m]=P(st[x],st[y],-1);
if(op==3){
if(d[x]>d[y])swap(x,y);
if(st[x]<=st[y]&&en[y]<=en[x])a[++m]=P(lca2(x,y),y,1);else a[++m]=P(x,y,2);
}
}
solve(1,m);
for(i=1;i<=m;i++)if(a[i].t>0)printf("%d
",ans[i]);
return 0;
}