首页 > 编程语言 >树的直径算法

树的直径算法

时间:2022-12-26 13:22:46浏览次数:26  
标签:int dfs 算法 delta 直径 Rightarrow d2 d1

树上任意两节点之间最长的简单路径即为树的「直径」。记 \(\delta(s, t)\) 为树的真实直径。

算法1,两遍dfs

首先从任意节点 \(y\) 开始进行第一次 DFS,到达距离其最远的节点,记为 \(z\),然后再从 \(z\) 开始做第二次 DFS,到达距离 \(z\) 最远的节点,记为 \(z'\),则 即为树的直径。

那么为什么此算法可行?

我们只需证明对于任意一点 \(x\) 经过第一次dfs之后到达的点 \(y\) 一定是 \(\delta(s,t)\) 的一端 \(s\) 或 \(t\)。

定理:在一棵树上,从任意节点 \(x\) 开始进行一次 DFS,到达的距离其最远的节点 \(z\) 必为直径的一端。

我们采用反证法。记出发点为 \(y\),真实直径 \(\delta(s,t)\),设从 \(y\) 进行的第一次dfs到达的距离其最远的点 \(z\) 不为 \(t\) 或 \(s\)。
那么有三种情况:

  • 若 \(y\) 在 \(\delta(s,t)\) 上:那么 \(\delta(y,z)>\delta(y,t)\Rightarrow \delta(x,z)>\delta(x,t)\Rightarrow \delta(s,z)>\delta(s,t)\) 矛盾。其中 \(x\) 是 \(\delta(s,t)\) 上一中转点
    image
  • 若 \(y\) 不在直径上,且 \(\delta(y,z)\) 与 \(\delta(s,t)\) 有重合路径
    image
    那么 \(\delta(y,z)>\delta(y,t)\Rightarrow \delta(x,z)>\delta(x,t)\Rightarrow \delta(s,z)>\delta(s,,t)\) 矛盾。
  • 若 \(y\) 不在 \(\delta(s,t)\) 上,且 \(\delta(y,z)\) 与 \(\delta(s,t)\) 不存在重复路径:
    image
    那么 \(\delta(y,z)>\delta(y,t)\Rightarrow \delta(x',z)>\delta(x',t)\Rightarrow\delta(x,z)>\delta(x,t)\Rightarrow\delta(s,z)>\delta(s,t)\)矛盾。
  • 三种情况都矛盾,得证。

注意负权边不适用

那么代码实现就很简单了

const int N = 10005;

int n, c, d[N];
vector<int> e[N];

void dfs(int u, int fa) {
  for (int v :e[u]) {
    if (v == fa) continue;
    d[v] = d[u] + 1;
    if (d[v] > d[c]) c = v;
    dfs(v, u);
  }
}

int main() {
  scanf("%d", &n);
  for (int i = 1; i < n; i++) {
    int u, v;
    scanf("%d %d", &u, &v);
    e[u].push_back(v), e[v].push_back(u);
  }
  dfs(1, 0);
  d[c] = 0, dfs(c, 0);
  printf("%d\n", d[c]);
  return 0;
}

\(O(N)\)?

树形dp

每个节点作为子树的根向下,所能延伸的最远距离 \(d1\),和次远距离 \(d2\),那么直径就是所有 \(d1+d2\) 的最大值。

树形 DP 可以在存在负权边的情况下求解出树的直径。

const int N=10010,M=20010;
int n,a,b,c,ans;
struct edge{int v,w;};
vector<edge> e[N];

int dfs(int x,int fa){
  int d1=0,d2=0; 
  for(auto ed : e[x]){
    int y=ed.v, z=ed.w;
    if(y==fa) continue;
    int d=dfs(y,x)+z;
    if(d>=d1) d2=d1,d1=d;
    else if(d>d2) d2=d;
    // printf("回%d d1=%d d2=%d\n",x,d1,d2);
  }
  ans=max(ans,d1+d2);
  // printf("离%d ans=%d\n",x,ans);
  return d1;
}
int main(){
  cin>>n;
  for(int i=1; i<n; i++){
    cin>>a>>b>>c;
    e[a].push_back({b,c});
    e[b].push_back({a,c});
  }
  dfs(1,-1);
  cout << ans << endl;
  return 0;
}

最后的最后,始终要记住,树是无环的,这一性质很重要,在两个算法中都用到了。

标签:int,dfs,算法,delta,直径,Rightarrow,d2,d1
From: https://www.cnblogs.com/CYLSY/p/17005582.html

相关文章

  • 每日算法之左旋转字符串
    JZ58左旋转字符串题目汇编语言中有一种移位指令叫做循环左移(ROL),现在有个简单的任务,就是用字符串模拟这个指令的运算结果。对于一个给定的字符序列S,请你把其循环左移......
  • C++不知算法系列之解析回溯算法中的人文哲学
    1.前言回溯算法让我想起“退一步海阔天空”的名言。当事情的发展到了绝境或是边缘时,可以试着后退一步,换一个方向、换一种策略,或许会看到新的出路或生机。回溯算法的精髓:......
  • 共识算法——Paxos算法
    故事Lamport描述了一个名为Paxos的希腊城邦(算法得名于此),这个城邦是按照民主的议会制度来进行选举的,所有的居民进行提议和投票来选出决议。但是居民们不想花时间一直在选举上......
  • Paxos算法原理及理解
    Paxos算法产生的背景Paxos算法是基于消息传递且具有高度容错特性的一致性算法,是目前公认的解决分布式一致性问题最有效的算法之一,其解决的问题就是在分布式系统中如何就某个......
  • 代码随想录算法训练营Day24|77. 组合
    代码随想录算法训练营Day24|77.组合回溯基础回溯的本质是穷举,穷举所有可能,然后选出我们想要的答案,常见的问题类型为:组合问题:N个数里面按一定规则找出k个数的集合切割......
  • 代码随想录算法训练营Day25|216. 组合总和 III、17. 电话号码的字母组合
    代码随想录算法训练营Day25|216.组合总和III、17.电话号码的字母组合216.组合总和III216.组合总和III与「77.组合」类似,但区别在于题干要求的变化:只使用数字1......
  • 周而复始,往复循环,递归、尾递归算法与无限极层级结构的探究和使用(Golang1.18)
    所有人都听过这样一个歌谣:从前有座山,山里有座庙,庙里有个和尚在讲故事:从前有座山。。。。,虽然这个歌谣并没有一个递归边界条件跳出循环,但无疑地,这是递归算法最朴素的落地实......
  • 使用键值对数组构造的无重复随机数算法
    --@paramlist_length生成的数组长度--@parammax_random_length随机数的最大范围math.generate=function(list_length,max_random_length) localrandom={}......
  • 精通visual c++指纹模式识别系统算法及实现
    通过学习,掌握以下几个问题:1、核心算法,并且向GVF衍生;2、核心库封装的方法2016年11月16日06:52:51昨日实现了梯度场和频率场的计算。最大的感觉就是建立基础代码库的重要性。......
  • ()找圆算法((HoughCircles)总结与优化
       Opencv内部提供了一个基于Hough变换理论的找圆算法,HoughCircle与一般的拟合圆算法比起来,各有优势:优势:HoughCircle对噪声点不怎么敏感,并且可以在同一个图中......