poj 1330 LCA

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 using namespace std;
 const int MAXN=;
 int rmq[*MAXN];//rmq数组，就是欧拉序列对应的深度序列
 struct ST
 {
     int mm[*MAXN];
     int dp[*MAXN][];//最小值对应的下标
     void init(int n)
     {
         mm[]=-;
         for(int i=;i<=n;i++)
         {
             mm[i]=((i&(i-))==)?mm[i-]+:mm[i-];
             dp[i][]=i;
         }
         for(int j=;j<=mm[n];j++)
         for(int i=;i+(<<j)-<=n;i++)
         dp[i][j]=rmq[dp[i][j-]]<rmq[dp[i+(<<(j-))][j-]]?dp[i][j-]:dp[i+(<<(j-))][j-];
     }
     int query(int a,int b)//查询[a,b]之间最小值的下标
     {
         if(a > b)swap(a,b);
         int k=mm[b-a+];
         return rmq[dp[a][k]]<=rmq[dp[b-(<<k)+][k]]?dp[a][k]:dp[b-(<<k)+][k];
     }
 };
 //边的结构体定义
 struct Edge
 {
     int to,next;
 };
 Edge edge[MAXN*];
 int tot,head[MAXN];
 int F[MAXN*];//欧拉序列，就是dfs遍历的顺序，长度为2*n-1,下标从1开始
 int P[MAXN];//P[i]表示点i在F中第一次出现的位置
 int cnt;
 ST st;
 void init()
 {
     tot=;
     memset(head,-,sizeof(head));
 }
 void addedge(int u,int v)//加边，无向边需要加两次
 {
     edge[tot].to=v;
     edge[tot].next=head[u];
     head[u]=tot++;
 }
 void dfs(int u,int pre,int dep)
 {
     F[++cnt]=u;
     rmq[cnt]=dep;
     P[u]=cnt;
     for(int i=head[u];i!=-;i=edge[i].next)
     {
         int v=edge[i].to;
         if(v==pre)continue;
         dfs(v,u,dep+);
         F[++cnt]=u;
         rmq[cnt]=dep;
     }
 }
 void LCA_init(int root,int node_num)//查询LCA前的初始化
 {
     cnt=;
     dfs(root,root,);
     st.init(*node_num-);
 }
 int query_lca(int u,int v)//查询u,v的lca编号
 {
     return F[st.query(P[u],P[v])];
 }
 bool flag[MAXN];
 int main()
 {
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     int T;
     int N;
     int u,v;
     scanf("%d",&T);
     while(T--)
     {
         scanf("%d",&N);
         init();
         memset(flag,false,sizeof(flag));
         for(int i=;i<N;i++)
         {
             scanf("%d%d",&u,&v);
             addedge(u,v);
             addedge(v,u);
             flag[v]=true;
         }
         int root;
         for(int i=;i<=N;i++)
         if(!flag[i])
         {
             root=i;
             break;
         }
         LCA_init(root,N);
         scanf("%d%d",&u,&v);
         printf("%d\n",query_lca(u,v));
     }
     return ;
 }