Nearest Common Ancestors

poj1330:http://poj.org/problem?id=1330

题意：求一棵树上的两点的最近的公共祖先。

题解：第一次接触LCA，第一道模板题。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <vector>
 using namespace std;
 #define MAXN 10001
 int n,fa[MAXN];
 int indegree[MAXN];
 int vis[MAXN];
 vector<int> hash[MAXN];
 struct Node{
    int v;
    int next;
    int id;
 }edge[MAXN*];
 int head[MAXN],cnt;
 int ances[MAXN],ans[MAXN];//祖先
 void init(){
   memset(head,-,sizeof(head));
   memset(ans,-,sizeof(ans));
   cnt=;
 }
 void add(int u,int v,int w){
     edge[cnt].v=v;
     edge[cnt].id=w;
     edge[cnt].next=head[u];
     head[u]=cnt++;
      edge[cnt].v=u;
     edge[cnt].id=w;
     edge[cnt].next=head[v];
     head[v]=cnt++;
 }

 void init(int n){
     for(int i=;i<=n;i++){
         fa[i]=i;
         indegree[i]=;
         vis[i]=;
         ances[i]=;
         hash[i].clear();
        // Qes[i].clear();
     }
 }
 int find(int x){
    int s;
    for(s=x;s!=fa[s];s=fa[s]);
    while(s!=x){
      int temp=fa[x];
      fa[x]=s;
      x=temp;
    }
    return s;
 }
 void unio(int x,int y){
     int fx=find(x),fy=find(y);
     if(fx==fy) return ;
      fa[fy]=fx;
 }
 void Tarjan(int u){
     ances[u]=u;
     int i,size = hash[u].size();
     for(i=;i<size;i++){
         Tarjan(hash[u][i]);
         unio(u,hash[u][i]);
         ances[find(u)]=u;
     }
     vis[u]=;
     for(int i=head[u];i!=-;i=edge[i].next){
         int v=edge[i].v;
         if(vis[v]==){
             ans[edge[i].id]=ances[find(v)];
             return;
         }
     }

     /*size = Qes[u].size();
     for(i=0;i<size;i++){
         if(vis[Qes[u][i]]==1){
             printf("%d\n",ances[find(Qes[u][i])]);
             return;
         }
     }*/
 }
 int main(){
     int t;
     int i,j;
     scanf("%d",&t);
     while(t--){
         scanf("%d",&n);
         init(n);
         int s,d;
         for(i=;i<=n-;i++){
             scanf("%d%d",&s,&d);
             hash[s].push_back(d);
             indegree[d]++;
         }
             scanf("%d%d",&s,&d);
             init();
             add(s,d,);
             add(d,s,);
         for(j=;j<=n;j++){
             if(indegree[j]==){
                 Tarjan(j);
                 break;
             }
         }
         printf("%d\n",ans[]);
     }
     return ;
 }