zoukankan      html  css  js  c++  java
  • tarjan算法求LCA

    tarjan算法求LCA

    LCA(Least Common Ancestors)的意思是最近公共祖先,即在一棵树中,找出两节点最近的公共祖先。

    这里我们使用tarjan算法离线算法解决这个问题。

    离线算法,是指首先读入所有的询问(求一次LCA叫做一次询问),然后重新组织查询处理顺序以便得到更高效的处理方法。Tarjan算法是一个常见的用于解决LCA问题的离线算法,它结合了深度优先遍历和并查集,整个算法为线性处理时间。

    总思路就是每进入一个节点u的深搜,就把整个树的一部分看作以节点u为根节点的小树,再搜索其他的节点。每搜索完一个点后,如果该点和另一个已搜索完点为需要查询LCA的点,则这两点的LCA为另一个点的现在的祖先。

    1.先建立两个链表,一个为树的各条边,另一个是需要查询最近公共祖先的两节点。

    2.建好后,从根节点开始进行一遍深搜。

    3.先把该节点u的father设为他自己(也就是只看大树的一部分,把那一部分看作是一棵树),搜索与此节点相连的所有点v,如果点v没被搜索过,则进入点v的深搜,深搜完后把点v的father设为点u。

    4.深搜完一点u后,开始判断节点u与另一节点v是否满足求LCA的条件,满足则将结果存入数组中。

    5.搜索完所有点,自动退出初始的第一个深搜,输出结果。

    如上图,根据实现算法可以看出,只有当某一棵子树全部遍历处理完成后,才将该子树的根节点标记为有颜色(初始化是白色),假设程序按上面的树形结构进行遍历,首先从节点1开始,然后递归处理根为2的子树,当子树2处理完毕后,节点2, 5, 6均为红色;接着要回溯处理3子树,首先被染色的是节点7(因为节点7作为叶子不用深搜,直接处理),接着节点7就会查看所有询问(7, x)的节点对,假如存在(7, 5),因为节点5已经被染黑,所以就可以断定(7, 5)的最近公共祖先就是find(5),即节点1(因为2子树处理完毕后,子树2和节点1返回了合并后的树的根1,此时树根的祖先的值就是1)。

    代码实现:

     1 #include<cstdio>
     2 #define N 420000
     3 struct hehe{
     4     int next;
     5     int to;
     6     int lca;
     7 };
     8 hehe edge[N];//树的链表
     9 hehe qedge[N];//需要查询LCA的两节点的链表
    10 int n,m,p,x,y;
    11 int num_edge,num_qedge,head[N],qhead[N];
    12 int father[N];
    13 int visit[N];//判断是否被找过 
    14 void add_edge(int from,int to){//建立树的链表 
    15     edge[++num_edge].next=head[from];
    16     edge[num_edge].to=to;
    17     head[from]=num_edge;
    18 }
    19 void add_qedge(int from,int to){//建立需要查询LCA的两节点的链表 
    20     qedge[++num_qedge].next=qhead[from];
    21     qedge[num_qedge].to=to;
    22     qhead[from]=num_qedge;
    23 }
    24 int find(int z){//找爹函数 
    25     if(father[z]!=z)
    26         father[z]=find(father[z]);
    27     return father[z];
    28 }
    29 int dfs(int x){//把整棵树的一部分看作以节点x为根节点的小树 
    30     father[x]=x;//由于节点x被看作是根节点,所以把x的father设为它自己 
    31     visit[x]=1;//标记为已被搜索过 
    32     for(int k=head[x];k;k=edge[k].next)//遍历所有与x相连的节点 
    33         if(!visit[edge[k].to]){//若未被搜索 
    34             dfs(edge[k].to);//以该节点为根节点搞小树 
    35             father[edge[k].to]=x;//把x的孩子节点的father重新设为x 
    36         }
    37     for(int k=qhead[x];k;k=qedge[k].next)//搜索包含节点x的所有询问 
    38         if(visit[qedge[k].to]){//如果另一节点已被搜索过 
    39             qedge[k].lca=find(qedge[k].to);//把另一节点的祖先设为这两个节点的最近公共祖先 
    40             if(k%2)//由于将每一组查询变为两组,所以2n-1和2n的结果是一样的 
    41                 qedge[k+1].lca=qedge[k].lca;
    42             else
    43                 qedge[k-1].lca=qedge[k].lca;
    44         }
    45 }
    46 int main(){
    47     scanf("%d%d%d",&n,&m,&p);//输入节点数,查询数和根节点 
    48     for(int i=1;i<n;++i){
    49         scanf("%d%d",&x,&y);//输入每条边 
    50         add_edge(x,y);
    51         add_edge(y,x);
    52     }
    53     for(int i=1;i<=m;++i){
    54         scanf("%d%d",&x,&y);//输入每次查询,考虑(u,v)时若查找到u但v未被查找,所以将(u,v)(v,u)全部记录 
    55         add_qedge(x,y);
    56         add_qedge(y,x);
    57     }
    58     dfs(p);//进入以p为根节点的树的深搜 
    59     for(int i=1;i<=m;i++)
    60         printf("%d ",qedge[i*2].lca);//两者结果一样,只输出一组即可 
    61     return 0;
    62 }
    View Code

     时间复杂度:O(n+m)。

    离线算法。

  • 相关阅读:
    SSH连接linux的centos报The host '192.168.*.*' is unreachable
    页面格式化数值
    学习webservice之cxf(8):Spring整合CXF
    学习webservice之cxf(7):cxf自定义拦截器
    学习webservice之cxf(6):cxf内置拦截器
    学习webservice之cxf(5):cxf处理map等复杂类型
    学习webservice之cxf(4):cxf处理javabean以及复合类型
    学习webservice之cxf(3):使用webservice实现客户端
    学习webservice之cxf(1):使用cxf实现webservice(使用jdk1.8)
    学习webservice之cxf(2):使用java代码使用webservice(jdk1.8无法使用)
  • 原文地址:https://www.cnblogs.com/jsawz/p/6723221.html
Copyright © 2011-2022 走看看