树的直径
树上任意两节点之间最长的简单路径即为树的「直径」。
可以用两次深度优先搜索或者动态规划求解
两次DFS
首先从任意节点 y 开始进行第一次 DFS,到达距离其最远的节点,记为 z,然后再从 z 开始做第二次 DFS,到达距离 z 最远的节点,记为 z’,则 \delta(z,z’) 即为树的直径。
显然,如果第一次 DFS 到达的节点 z 是直径的一端,那么第二次 DFS 到达的节点 z’ 一定是直径的一端。我们只需证明在任意情况下,z 必为直径的一端。
定理:在一棵树上,从任意节点 y 开始进行一次 DFS,到达的距离其最远的节点 z 必为直径的一端。
const int N = 10000 + 10;
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;
}
树形DP
我们记录当 1 为树的根时,每个节点作为子树的根向下,所能延伸的最长路径长度 d_1 与次长路径(与最长路径无公共边)长度 d_2,那么直径就是对于每一个点,该点 d_1 + d_2 能取到的值中的最大值。
树形 DP 可以在存在负权边的情况下求解出树的直径。
const int N = 10000 + 10;
int n, d = 0;
int d1[N], d2[N];
vector<int> E[N];
void dfs(int u, int fa) {
d1[u] = d2[u] = 0;
for (int v : E[u]) {
if (v == fa) continue;
dfs(v, u);
int t = d1[v] + 1;
if (t > d1[u])
d2[u] = d1[u], d1[u] = t;
else if (t > d2[u])
d2[u] = t;
}
d = max(d, d1[u] + d2[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);
printf("%d\n", d);
return 0;
}
例题
给你 n 个城市,编号为从 1 到 n 。同时给你一个大小为 n-1 的数组 edges ,其中 edges[i] = [ui, vi] 表示城市 ui 和 vi 之间有一条双向边。题目保证任意城市之间只有唯一的一条路径。换句话说,所有城市形成了一棵 树 。
一棵 子树 是城市的一个子集,且子集中任意城市之间可以通过子集中的其他城市和边到达。两个子树被认为不一样的条件是至少有一个城市在其中一棵子树中存在,但在另一棵子树中不存在。
对于 d 从 1 到 n-1 ,请你找到城市间 最大距离 恰好为 d 的所有子树数目。
请你返回一个大小为 n-1 的数组,其中第 d 个元素(下标从 1 开始)是城市间 最大距离 恰好等于 d 的子树数目。
请注意,两个城市间距离定义为它们之间需要经过的边的数目。
示例 1:
输入:n = 4, edges = [[1,2],[2,3],[2,4]]
输出:[3,4,0]
解释:
子树 {1,2}, {2,3} 和 {2,4} 最大距离都是 1 。
子树 {1,2,3}, {1,2,4}, {2,3,4} 和 {1,2,3,4} 最大距离都为 2 。
不存在城市间最大距离为 3 的子树。
示例 2:
输入:n = 2, edges = [[1,2]]
输出:[1]
示例 3:
输入:n = 3, edges = [[1,2],[2,3]]
输出:[2,1]
提示:
2 <= n <= 15
edges.length == n-1
edges[i].length == 2
1 <= ui, vi <= n
题目保证 (ui, vi) 所表示的边互不相同。
class Solution {
public:
vector<int> countSubgraphsForEachDiameter(int n, vector<vector<int>>& edges) {
vector<vector<int>> adj(n);
for (auto &edge : edges) {
int x = edge[0] - 1;
int y = edge[1] - 1;
adj[x].emplace_back(y);
adj[y].emplace_back(x);
}
function<int(int, int, int)> dfs = [&](int parent, int u, int mask)->int {
int depth = 0;
for (int v : adj[u]) {
if (v != parent && mask & (1 << v)) {
depth = max(depth, 1 + dfs(u, v, mask));
}
}
return depth;
};
vector<int> ans(n - 1);
for (int i = 1; i < (1 << n); i++) {
int x = 32 - __builtin_clz(i) - 1;
int mask = i;
int y = -1;
queue<int> qu;
qu.emplace(x);
mask &= ~(1 << x);
while (!qu.empty()) {
y = qu.front();
qu.pop();
for (int vertex : adj[y]) {
if (mask & (1 << vertex)) {
mask &= ~(1 << vertex);
qu.emplace(vertex);
}
}
}
if (mask == 0) {
int diameter = dfs(-1, y, i);
if (diameter > 0) {
ans[diameter - 1]++;
}
}
}
return ans;
}
};