2021年度训练联盟热身训练赛第一场 I. Full Depth Morning Show

链接

https://ac.nowcoder.com/acm/contest/12606/I

题意

记 $d_{u,v}$ 为节点 $u$ 到节点 $v$ 的最短距离，$t_u$ 为节点 $u$ 的权值，对树上任意节点 $u$ 求 ${\sum }_{v=1}^n[d_{u,v}\ast (t_u+t_v)]$。

思路

换根 DP。

记 $re{s}_u={\sum }_{v=1}^n[d_{u,v}\ast (t_u+t_v)]$：

$re{s}_u={\sum }_{v=1}^n[d_{u,v}\ast (t_u+t_v)]=t_u\ast {\sum }_{v=1}^n{d}_{u,v}+{\sum }_{v=1}^n(d_{u,v}\ast t_v)$。

记 $a_u={\sum }_{v=1}^n{d}_{u,v}$，$b_u={\sum }_{v=1}^n(d_{u,v}\ast t_v)$，$s{z}_u$ 为节点 $u$ 的子树节点数，$su{m}_u$ 为节点 $u$ 的子树节点权值和，$u^{\prime }$ 为节点 $u$ 的儿子节点：

$a_{u^{\prime }}=a_u-d_{u,u^{\prime }}\ast s{z}_u+d_{u,u^{\prime }}\ast (n-s{z}_{u^{\prime }})=a_u+d_{u,u^{\prime }}\ast (n-2\ast s{z}_{u^{\prime }})$，

$b_{u^{\prime }}=b_u-d_{u,u^{\prime }}\ast su{m}_u+d_{u,u^{\prime }}\ast ({\sum }_{i=1}^n{t}_i-su{m}_{u^{\prime }})=b_u+d_{u,u^{\prime }}\ast ({\sum }_{i=1}^n{t}_i-2\ast su{m}_{u^{\prime }})$，

$re{s}_{u^{\prime }}=t_{u^{\prime }}\ast a_{u^{\prime }}+b_{u^{\prime }}$。

两次 DFS，以 节点 $1$ 为根，先计算出 $a_1$，$b_1$，再推出其余节点的 $res$，$O(n)$ 求解。

代码

#include <bits/stdc++.h>
#define SZ(x) (int)(x).size()
#define ALL(x) (x).begin(),(x).end()
#define PB push_back
#define EB emplace_back
#define MP make_pair
#define FI first
#define SE second
using namespace std;
typedef double DB;
typedef long double LD;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
typedef vector<int> VI;
typedef vector<PII> VPII;
// head
const int N=1e5+5;
int n;
LL t[N],d[N],sz[N],sum[N],a[N],b[N],res[N];
VPII g[N];
void dfs1(int u,int fa) {
    sz[u]=1;
    sum[u]=t[u];
    for(auto x:g[u]) {
        int v=x.FI,w=x.SE;
        if(v==fa) continue;
        d[v]=d[u]+w;
        a[1]+=d[v];
        b[1]+=d[v]*t[v];
        dfs1(v,u);
        sz[u]+=sz[v];
        sum[u]+=sum[v];
    }
}
void dfs2(int u,int fa) {
    for(auto x:g[u]) {
        int v=x.FI,w=x.SE;
        if(v==fa) continue;
        a[v]=a[u]+w*(n-2*sz[v]);
        b[v]=b[u]+w*(sum[1]-2*sum[v]);
        res[v]=a[v]*t[v]+b[v];
        dfs2(v,u);
    }
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin>>n;
    for(int i=1;i<=n;i++) cin>>t[i];
    for(int i=1;i<n;i++) {
        int u,v,w;
        cin>>u>>v>>w;
        g[u].EB(v,w);
        g[v].EB(u,w);
    }
    dfs1(1,0);
    dfs2(1,0);
    res[1]=a[1]*t[1]+b[1];
    for(int i=1;i<=n;i++) cout<<res[i]<<'\n';
    return 0;
}