poj 1330(LCA)tarjam运用

Poj1330

题意:求一棵树中两点的最近祖先。

思路:直接用tarjan算法。(还可以求出强连通分量)

下面是我对这个算法的理解,假如我们是求u,v的最近公共祖先。

Tarjan主要利用了dfs的性质。我们在进行dfs时,如果我们已经遍历了u,当我们遍历v时,这时他们最近的公共祖先也在栈里边。这时间我们可以在设一个f[i]数组,这个数组是用来确定i的祖先,我们在遍历的时间可以对f进行控制,利用并查集找的答案。



 //最近公共祖先tarjan算法
 #include<stdio.h>
 #include<string.h>
 #include<vector>
 #include<iostream>
 using namespace std;
 const int MAX = 10010;
 vector<int> edge[MAX];
 vector<int> query[MAX];//查询记录与u节点连接的的v
 int f[MAX],n,indegree[MAX];
 int lca[MAX],vis[MAX];


int findx(int x)
{
      int r = x;int i,j;
      while (f[r] != r) //循环结束,则找到根节点
          r = f[r];
        i = x;
      while (i != r) //本循环修改查找路径中所有节点
      {
          j = f[i];
         f[i] = r;
          i = j;
      }
     return r;
}


void Union(int x,int y)
{
    int fx,fy;
    fx = findx(x);
    fy = findx(y);
    if(fx <= fy)
        f[fy] = fx;
        else f[fx]=fy;
}
 void dfs(int u)
 {
      int  i;
      lca[u]=u;
      for(i=0;i<edge[u].size();i++)
      {
           //cout<<u<<endl;
           dfs(edge[u][i]);//以下为回溯部分
   //        cout<<edge[u][i]<<endl;
           Union(u,edge[u][i]);
           lca[findx (u)]=u;
      }
      vis[u]=1;
      for(i=0;i<query[u].size();i++)//遍历过则说明已找到答案
      {
          if(vis[query[u][i]])
          {
              printf("%d\n",lca[findx (query[u][i])]);
              return ;
          }
      }
 }


 void init()
 {
     int i;
     for(i=0;i<=n;i++)
     {
         edge[i].clear();
         query[i].clear();
         indegree[i]=0;vis[i]=0;
         f[i]=i;lca[i]=0;
     }
 }
 int main()
 {
     int t,i,a,b;
     scanf("%d",&t);
     while(t--)
     {
         scanf("%d",&n);
         init();
         for(i=1;i<n;i++)
         {
             scanf("%d%d",&a,&b);
             indegree[b]++;
             edge[a].push_back(b);
         }
         scanf("%d%d",&a,&b);
         query[a].push_back(b);
         query[b].push_back(a);
         for(i=1;i<=n;i++)
         {
              if(!indegree[i])
              {
                  dfs(i);
                  break;
              }
         }
     }
     return 0;
 }

### POJ 1330 题解 POJ 1330 是一道经典的最近公共祖先 (Least Common Ancestor, LCA) 问题。该题的核心在于如何高效地求解两节点之间的最近公共祖先。 #### 倍增算法简介 倍增方法是一种高效的在线求解 LCA 的方式,其时间复杂度为 \(O((n+q)\log n)\),其中 \(n\) 表示树中的节点数量,\(q\) 表示查询次数。此方法通过预处理的方式记录每个节点向上跳跃若干步后的父节点位置,从而加速查询过程。 具体来说,定义数组 `fa[i][j]` 表示从节点 \(i\) 向上跳 \(2^j\) 步所到达的节点编号。为了支持这一操作,我们需要满足如下递推关系: \[ \text{fa}[i][j] = \begin{cases} \text{parent}(i), & j = 0 \\ \text{fa}[\text{fa}[i][j-1]][j-1], & j > 0 \end{cases} \] 这种预处理可以通过深度优先搜索 (DFS) 完成,在 DFS 过程中同时计算每个节点的深度以及它们的倍增父节点表。 #### 查询阶段 在实际查询过程中,假设我们要查找节点 \(a\) 和 \(b\) 的最近公共祖先,则分为以下几个部分: 1. **调整深度一致** 如果当前两个节点的深度不相同,则将较深的节点不断向上提升至两者深度相等的位置。这一步利用了倍增数组 `fa[i][j]` 来快速跳跃多个层次。 2. **同步爬升** 当两个节点处于同一深度时,让它们逐步向上移动直至相遇于某个共同祖先处。为了避免逐层遍历带来的效率低下问题,同样采用倍增策略按指数级跳跃。 最终返回的结果即为最后一次未分叉前的状态作为答案。 以下是基于以上逻辑的一个简单实现代码片段: ```cpp #include <bits/stdc++.h> using namespace std; const int MAXN = 5e4 + 5; vector<int> adj[MAXN]; int depth_[MAXN]; // 存储每个结点的深度 int fa[MAXN][20]; // 记录每个结点向上跳2^k步的父亲 void dfs(int u, int parent){ depth_[u]=depth_[parent]+1; fa[u][0]=parent; for(auto v : adj[u]){ if(v != parent){ dfs(v,u); } } } // 初始化fa[][]表格 void preprocess(int N){ memset(fa,-1,sizeof(fa)); for(int k=1;k<20;k++){ for(int i=1;i<=N;i++){ if(fa[i][k-1]!=-1){ fa[i][k]=fa[fa[i][k-1]][k-1]; } } } } int get_lca(int a,int b){ if(depth_[a]<depth_[b]) swap(a,b); // 将a提到和b相同的高度 for(int k=19;k>=0;k--){ if(fa[a][k]!=-1 && depth_[fa[a][k]] >= depth_[b]){ a=fa[a][k]; } } if(a==b)return a; // 同步爬升 for(int k=19;k>=0;k--){ if(fa[a][k]!=-1 && fa[b][k]!=-1 && fa[a][k]!=fa[b][k]){ a=fa[a][k]; b=fa[b][k]; } } return fa[a][0]; } ``` 上述程序实现了完整的倍增法用于解决LCA问题的功能模块化设计[^2]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值