zoukankan      html  css  js  c++  java
  • LCA---模板

    LCA:最近公共祖先

    性质:1、LCA( u ) = u。

               2、u 是 v 的祖先,当且仅当 LCA( u , v ) = u。

               3、如果 u 不为 v 的祖先并且 v 不为 u 的祖先,那么这两个点属于不同的子树。

               4、两点的最近公共祖先必定处在两点间的树上最短路上。

               5、d( u , v ) = h( u ) + h( v ) - 2*h( LCA( u , v )) ,其中 d 是树上两点间的最短距离, h 是该点与树根之间的距离。

    用途:求树从 x 到 y 节点最短路径上所有节点的值之和。

    代码:

     1 #include<bits/stdc++.h>
     2 #define ll long long
     3 #define MOD 1e9+7 
     4 #define INF 0x3f3f3f3f
     5 #define mem(a,x) memset(a,x,sizeof(a))  
     6 #define _for(i,a,b) for(int i=a; i< b; i++)
     7 #define _rep(i,a,b) for(int i=a; i<=b; i++)
     8 #define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
     9 
    10 using namespace std;
    11 const int MAXN = 1000005 ;
    12 inline int rd() 
    13 {
    14     int res = 0,flag = 0;
    15     char ch;
    16     if ((ch = getchar()) == '-')flag = 1;
    17     else if(ch >= '0' && ch <= '9')res = ch - '0';
    18     while ((ch = getchar()) >= '0' && ch <= '9')res = (res<<1) + (res<<3) + (ch - '0');
    19     return flag ? -res : res;
    20 }
    21 //
    22 int head[MAXN]; 
    23 int num=0;
    24 struct edg{
    25    int next,to;
    26 }edge[MAXN];
    27 void edge_add(int u,int v)   //链式前向星存图
    28 {
    29    num++;
    30    edge[num].next=head[u];edge[num].to=v;head[u]=num;
    31    edge[++num].next=head[v];edge[num].to=u;head[v]=num;
    32 }
    33 //---------------------
    34 // lca部分
    35 int n,m,s;
    36 int dep[MAXN]={0},f[MAXN][23];
    37 void dfs(int u,int father)//对应深搜预处理f数组 
    38 {
    39     dep[u]=dep[father]+1;
    40     for(int i=1;(1<<i)<=dep[u];i++)
    41     {
    42         f[u][i]=f[f[u][i-1]][i-1];
    43     }
    44     for(int i=head[u];i;i=edge[i].next)
    45     {
    46         int v=edge[i].to;
    47         if(v==father)continue;//双向图需要判断是不是父亲节点 
    48         f[v][0]=u;
    49         dfs(v,u);
    50     }
    51 }
    52 int lca(int x,int y)
    53 {
    54     if(dep[x]<dep[y])swap(x,y);
    55     for(int i=20;i>=0;i--)//从大到小枚举使x和y到了同一层 
    56     {
    57         if(dep[f[x][i]]>=dep[y])x=f[x][i];
    58         if(x==y)return x;
    59     }
    60     for(int i=20;i>=0;i--)//从大到小枚举 
    61     {
    62         if(f[x][i]!=f[y][i])//尽可能接近 
    63         {
    64             x=f[x][i];y=f[y][i];
    65         } 
    66     } 
    67     return f[x][0];//随便找一个**输出 
    68 }
    69 int main(){
    70     scanf("%d%d%d",&n,&m,&s);
    71     for(int i=1;i<n;i++)
    72     {
    73         scanf("%d",&a1);scanf("%d",&a2);
    74         edge_add(a1,a2);//链式存边 
    75     }
    76     dfs(s,0);
    77     for(int i=1;i<=m;i++)
    78     {
    79         scanf("%d %d",&a1,&a2);
    80         printf("%d
    ",lca(a1,a2));//求两个节点的LCA 
    81     }
    82     return 0;
    83 } 
    View Code
    越自律,越自由
  • 相关阅读:
    安装 macbook 双系统( OS X 和 Ubuntu )
    微信小程序初探
    [Javascript] Promise
    [AS/400] 基本概念
    [AS/400] Control Language(CL) 基本概念
    [AS/400] Control Language
    [github] 创建个人网页
    [tools] sublime 使用记录
    [Java] TreeMap
    exe转msi
  • 原文地址:https://www.cnblogs.com/ha-chuochuo/p/13522713.html
Copyright © 2011-2022 走看看