ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

图论专题-学习笔记:树的直径

2022-04-17 15:33:30  阅读:197  来源: 互联网

标签:f1 图论 专题 int 笔记 Edge 直径 now DP


目录

1. 前言

树的直径是树的一个小板块,但是有着重要的应用。

前置知识:树的基础知识。

2. 详解

例题:SP1437 PT07Z - Longest path in a tree

2.1 定义

树的直径:一棵树上最长的路径叫做树的直径。

比如下面这棵树,带有边权 1 的路径就是树的直径。

在这里插入图片描述

需要注意的是,这棵树的直径不止一条,但是一般情况下我们只取其中一条叫做这棵树的直径。

2.2 求法

那么怎么求树的直径呢?

这里有两种求法:DFS 与树形 DP。

2.2.1 DFS 求解

该算法的大致步骤如下:

  • 首先随便取一个点,做一遍 DFS,求出这个点能够到达的最远点。
  • 然后从这个点再次 DFS,求出这个点能够到达的最远点。
  • 求出来的两个点之间的路径就是树的直径。

步骤简明易懂,那么这个算法为什么是正确的呢?

下面证明假设边权大于 0。

采用反证法:

假设图中还存在一条比我们求出来的路径更长的路径,那么这条路径就是树的直径。

设我们求出的路径为 \(AB\),真正的直径为 \(CD\)。

分为两种情况:

  1. 直径与我们求出的路径相交或者部分重合。

在这里插入图片描述

那么根据上述算法,我们在第一次找到点 \(A\) 的时候必有 \(AE>CE\)。

同理,有 \(AE+EF+FB>AE+EF+FD\)。

考虑将 \(AE>CE\) 带入上述不等式,有 \(AE+EF+FB>CE+EF+FD\),这与 \(CD\) 是直径不符,所以上述情况不成立。

  1. 直径与我们求出的路径不重合。

在这里插入图片描述

则 \(CF+FD>CF+EF+EB\),因此 \(FD>EF+EB\),则有 \(FD>EB\)。

然而根据我们的算法步骤,\(EB>EF+FD\),则有 \(EB>FD\)。

出现了矛盾,因此原假设错误。

综上所述,\(AB\) 必须是直径。

那么为什么说这个证明必须有边权大于 0 呢?

这是因为如果边权小于 0,那么上述所有证明的 \(a>b+c\rightarrow a>c\) 就都不一定成立了。

同时这也揭示了 DFS 求树的直径必须有所有边权大于 0。

DFS 的优点:可以记录树的直径都有哪些点,也就是可以完整记录路径。

DFS 的缺点:不能处理有负边权的树。

2.2.2 树形 DP 求解

这个需要有一定树形 DP 基础,当然没有基础也没有问题。

考虑设 \(f1_i,f2_i\) 分别表示第 \(i\) 个节点到叶子节点路径的最大值与次大值。

那么设计状态转移方程:

设当前的点为 \(u\),当前枚举的儿子节点为 \(v\),边权为 \(val\),那么有如下方程:

  • 如果 \(f1_v+val>f1_u\),那么 \(f2_u\leftarrow f1_u,f1_u\leftarrow f1_v+val\)。
  • 否则,\(f2_u=\max\{f2_u,f1_v+val\}\)。

转移方程还是简明易懂的吧qwq

最后的答案就是 \(\max\{f1_i+f2_i|i \in [1,n]\}\)。

那么这个做法很明显,因为是 DP 思想,可以处理具有负边权的树。

树形 DP 的优点:可以处理具有负边权的树。

树形 DP 的缺点:不能记录路径。

2.3 代码

两种做法的代码如下(Solve_dfs 是 DFS 做法,Solve_DP 是树形 DP 做法):

/*
========= Plozia =========
    Author:Plozia
    Problem:SP1437 PT07Z - Longest path in a tree
    Date:2021/4/27
    Another:树的直径模板题
========= Plozia =========
*/

#include <bits/stdc++.h>

typedef long long LL;
const int MAXN = 1e4 + 10;
int n, Head[MAXN], cnt_Edge = 1, f[MAXN], ans, root, f1[MAXN], f2[MAXN];
struct node { int Next, to, val; } Edge[MAXN << 1];

int read()
{
    int sum = 0, fh = 1; char ch = getchar();
    for (; ch < '0' || ch > '9'; ch = getchar()) fh -= (ch == '-') << 1;
    for (; ch >= '0' && ch <= '9'; ch = getchar()) sum = (sum << 3) + (sum << 1) + (ch ^ 48);
    return sum * fh;
}
int Max(int fir, int sec) { return (fir > sec) ? fir : sec; }
void add_Edge(int x, int y, int z) { ++cnt_Edge; Edge[cnt_Edge] = (node){ Head[x], y, z }; Head[x] = cnt_Edge; }

void dfs(int now, int father)
{
    if (f[now] > ans) { ans = f[now]; root = now; }
    for (int i = Head[now]; i; i = Edge[i].Next)
    {
        int u = Edge[i].to;
        if (u == father) continue ;
        f[u] = f[now] + Edge[i].val;
        dfs(u, now);
    }
}

void Solve_DFS()
{
    root = 1;
    f[root] = 0; dfs(root, 0); ans = 0;
    f[root] = 0; dfs(root, 0);
    printf("%d\n", ans);
}

void DP(int now, int father)
{
    for (int i = Head[now]; i; i = Edge[i].Next)
    {
        int u = Edge[i].to;
        if (u == father) continue ;
        DP(u, now);
        if (f1[u] + Edge[i].val > f1[now]) { f2[now] = f1[now]; f1[now] = f1[u] + Edge[i].val; }
        else { f2[now] = Max(f2[now], f1[u] + Edge[i].val); }
    }
}

void Solve_DP()
{
    DP(1, 0);
    for (int i = 1; i <= n; ++i) ans = Max(ans, f1[i] + f2[i]);
    printf("%d\n", ans);
}

int main()
{
    n = read();
    for (int i = 1; i < n; ++i)
    {
        int u = read(), v = read();
        add_Edge(u, v, 1); add_Edge(v, u, 1);
    }
    // Solve_DFS(); return 0;
    Solve_DP(); return 0;
}

3. 总结

树的直径:树上最长路径。

求法:

  • 两遍 DFS 求解。
    • 优点:可以记录路径。
    • 缺点:不能处理负边权。
  • 树形 DP 求解。
    • 优点:可以处理负边权。
    • 缺点:不能处理路径。

标签:f1,图论,专题,int,笔记,Edge,直径,now,DP
来源: https://www.cnblogs.com/Plozia/p/16155839.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有