最近需要统计树形结构中任意两点间的最短距离,而且是针对大规模的树形结构,无法调用floyd最短路径算法,所以想要根据树形结构节点间路径的唯一性来求解。在网上搜到了相关的代码,故贴在这里。
转载,转自:
http://www.cnblogs.com/wally/archive/2013/06/03/3116020.html
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2376
思路:
引:如果暴力枚举两点再求距离是显然会超时的。转换一下思路,我们可以对每条边,求所有可能的路径经过此边的次数:设这条边两端的点数分别为A和B,那 么这条边被经过的次数就是A*B,它对总的距离和的贡献就是(A*B*此边长度)。我们把所有边的贡献求总和,再除以总路径数N*(N-1)/2,即为最 后所求。
每条边两端的点数的计算,实际上是可以用一次dfs解决的。任取一点为根,在dfs的过程中,对每个点k记录其子树包含的点数(包括其自身),设点数为a[k],则k的父亲一侧的点数即为N-a[k]。这个统计可以和遍历同时进行。故时间复杂度为O(n)。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
#define MAXN 10000+10
struct Node{
int v,len;
};
vector<Node>vet[MAXN];
int n;
double dp[MAXN];
int sum[MAXN];
void dfs(int root,int father){
sum[root]=1;
for(int i=0;i<vet[root].size();i++){
int son=vet[root][i].v;
int len=vet[root][i].len;
if(son==father)continue;
dfs(son,root);
sum[root]+=sum[son];
dp[root]+=dp[son]+(sum[son]*(n-sum[son]))*(double)len;
}
}
int main(){
// freopen("1.txt","r",stdin);
int _case,u,v,len;
scanf("%d",&_case);
while(_case--){
scanf("%d",&n);
for(int i=0;i<=n;i++)vet[i].clear();
for(int i=1;i<=n-1;i++){
scanf("%d%d%d",&u,&v,&len);
Node p1,p2;
p1.v=v,p2.v=u;
p1.len=p2.len=len;
vet[u].push_back(p1);
vet[v].push_back(p2);
}
memset(sum,0,sizeof(sum));
memset(dp,0,sizeof(dp));
dfs(0,-1);
int s=(n*(n-1)/2);
printf("%1lf\n",(double)dp[0]/s);
}
return 0;
}