51nod 1405 树的距离之和

给定一棵无根树,假设它有n个节点,节点编号从1到n, 求任意两点之间的距离(最短路径)之和。
Input
第一行包含一个正整数n (n <= 100000),表示节点个数。
后面(n - 1)行,每行两个整数表示树的边。
Output
每行一个整数,第i(i = 1,2,...n)行表示所有节点到第i个点的距离之和。
Input示例
4
1 2
3 2
4 2
Output示例
5
3
5
5
官方题解:
数学题,也可以说dp,不太难。
   
   
(1)我们给树规定一个根。假设所有节点编号是0-(n-1),我们可以简单地把0当作根,这样下来父子关系就确定了。
(2)定义数组num[x]表示以节点x为根的子树有多少个节点,dp[x]是我们所求的——所有节点到节点x的距离之和。
(3)在步骤(1)中,其实我们同时可以计算出 num[x],还可以计算出每个节点的深度(每个到根节点0的距离),累加全部节点深度得到的其实就是是dp[0]。
(4) 假设一个非根节点x,它的父亲节点是y, 并且dp[y]已经计算好了,我们如何计算dp[x]?  
以x为根的子树中那些节点,到x的距离比到y的距离少1, 这样的节点有num[x]个。
        其余节点到x的距离比到y的距离多1,这样的节点有(n - num[x])个。
        于是我们有 dp[x] = dp[y] - num[x] + (n - num[x])
                                      = dp[y] + n - num[x] * 2
        因为树的根节点dp[0]在步骤(3)已经计算出来了,根据所有的父子关系和这个上式,我们可以按照顺序计算出整个dp数组。
注意点: 重要的步骤都是简单的dfs,但是一半递归实现可能导致堆栈溢出。
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<map>
#include<vector>
#include<cstring>
#include<cmath>
#pragma comment(linker, "/STACK:10240000,10240000")
using namespace std;
typedef long long ll;
const int inf =0x3f3f3f3f;
const double  pi = acos(-1.0);
const int N = 1e5 + 10;
int root[N], num[N], len[N], vis[N], n;
vector<int>path[N];
ll dp[N];
void dfs1(int x)
{
    num[x]  = 1;
    for(int i = 0, sz = (int)path[x].size(); i<sz; i++)
    {
        int y = path[x][i];
        if(!vis[y])
        {
            vis[y] = 1;
            len[y] = len[x] + 1;
            dp[1] += len[y];
            dfs1(y);
            num[x] += num[y];
        }
    }
}
void dfs2(int x)
{
    for(int i = 0, sz = (int)path[x].size(); i<sz; i++)
    {
        int y = path[x][i];
        if(!vis[y])
        {
            dp[y] = dp[x] - num[y] + n - num[y];
            vis[y] = 1;
            dfs2(y);
        }
    }
}
int main()
{
    scanf("%d", &n);
    int u, v;
    for(int i = 1; i<n; i++)
    {
        scanf("%d%d", &u, &v);
        path[u].push_back(v);
        path[v].push_back(u);
    }
    vis[1] = 1;
    dfs1(1);
    memset(vis, 0, sizeof(vis));
  //  for(int i = 1; i<=n; i++)
  //      printf("%d\n", num[i]);
    vis[1] = 1;
    dfs2(1);
    for(int i = 1; i<=n; i++)
        printf("%I64d\n", dp[i]);
        return 0;
}





曹鹏  (题目提供者)
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值