求树的直径算法

/*树的直径是指树的最长简单路。求法: 两遍BFS :先任选一个起点BFS找到最长路的终点,再从终点进行BFS,则第二次BFS找到的最长路即为树的直径;
              原理: 设起点为u,第一次BFS找到的终点v一定是树的直径的一个端点
              证明: 1) 如果u 是直径上的点,则v显然是直径的终点(因为如果v不是的话,则必定存在另一个点w使得u到w的距离更长,则于BFS找到了v矛盾)
                      2) 如果u不是直径上的点,则u到v必然于树的直径相交(反证),那么交点到v 必然就是直径的后半段了
                       所以v一定是直径的一个端点,所以从v进行BFS得到的一定是直径长度
*/

#include <bits/stdc++.h>
using namespace std;


#define INF 10000000000

vector <int > G[1000005];
vector<int > E[1000005];
bool vis[1000005];
int d[1000005];

void init() {
	memset(vis, 0, sizeof(vis));
} 

void dfs(int u) {
	vis[u] = 1;
	int size = G[u].size();		//与顶点u相连的点数 
	
	for (int i = 0; i<size; i++) {			//对与顶点u相连的点数进行扫描 
		int v = G[u][i];
		if (!vis[v]) {
			d[v] = d[u] + E[u][i];
			dfs(v);
		}
	}
}

int main() {
	int n;
	cin >> n;
	int u, v, w;
	for (int i = 0; i<n-1; i++) {		//建立树过程 
		scanf("%d%d%d",&u,&v,&w);
		G[u-1].push_back(v-1);				//顶点两边都要记录 
		E[u-1].push_back(w);
		G[v-1].push_back(u-1);
		E[v-1].push_back(w);
	}
	
	init();
	for (int i = 0; i<n; i++)	d[i] = (i == 0?0:INF);
	dfs(0);
	int start = 0;
	int max = -1;
	for (int i = 0; i<n; i++) {
		if (d[i] > max && d[i] != INF) {
			max = d[i];
			start = i;
		}
	}
	
	init();
	for (int i = 0; i<n; i++)	d[i] = (i == start?0:INF);
	dfs(start);
	int ans = -1;
	for (int i = 0; i<n; i++) {
		if (d[i] > ans && d[i] != INF) {
			ans = d[i];
		}
	}
	
	//ans = 10*ans + ans*(ans+1)/2;
	cout << ans << endl;					//ans  即为直径 
	return 0;
}

转载于:https://www.cnblogs.com/Tovi/p/6194764.html

  • 1
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是Python求直径最小生成树的步骤: 1. 首先,需要使用Kruskal算法来生成最小生成树。 2. 然后,需要对最小生成树进行遍历,找到最长的边,即为直径。 3. 最后,输出直径和最小生成树。 下面是Python代码示例: ```python # 定义边的类 class Edge: def __init__(self, u, v, w): self.u = u self.v = v self.w = w # 定义并查集类 class UnionFind: def __init__(self, n): self.parent = list(range(n)) self.rank = [0] * n def find(self, x): if self.parent[x] != x: self.parent[x] = self.find(self.parent[x]) return self.parent[x] def union(self, x, y): px, py = self.find(x), self.find(y) if px == py: return False if self.rank[px] < self.rank[py]: px, py = py, px self.parent[py] = px if self.rank[px] == self.rank[py]: self.rank[px] += 1 return True # 定义Kruskal算法函数 def kruskal(n, edges): uf = UnionFind(n) edges.sort(key=lambda e: e.w) res = [] for e in edges: if uf.union(e.u, e.v): res.append(e) if len(res) == n - 1: break return res # 定义求直径函数 def diameter(n, edges): mst = kruskal(n, edges) dist = [-1] * n dist[mst[0].u] = 0 q = [mst[0].u] while q: u = q.pop(0) for e in mst: if e.u == u and dist[e.v] == -1: dist[e.v] = dist[u] + e.w q.append(e.v) elif e.v == u and dist[e.u] == -1: dist[e.u] = dist[u] + e.w q.append(e.u) u = max(range(n), key=lambda i: dist[i]) dist = [-1] * n dist[u] = 0 q = [u] while q: u = q.pop(0) for e in mst: if e.u == u and dist[e.v] == -1: dist[e.v] = dist[u] + e.w q.append(e.v) elif e.v == u and dist[e.u] == -1: dist[e.u] = dist[u] + e.w q.append(e.u) return max(dist) # 测试 n = 5 edges = [Edge(0, 1, 1), Edge(0, 2, 2), Edge(0, 3, 3), Edge(0, 4, 4), Edge(1, 2, 5), Edge(1, 3, 6), Edge(1, 4, 7), Edge(2, 3, 8), Edge(2, 4, 9), Edge(3, 4, 10)] print("最小生成树:", kruskal(n, edges)) print("直径:", diameter(n, edges)) ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值