A-Radio-Prize

  1. 1. 题意
  2. 2. 分析
  3. 3. 代码

题意

给定一棵树,树上每条边及每个点都有权值。定义$P_u=\sum_{v}d(u,v)(t_u+t_v)$其中$d(u,v)$为树上u到v的距离,$t_u$和$t_v$分别是两点权值,求每个点的P值。

分析

首先我们可以看到$P_u$可以分成 $t_u\sum_{v}d(u,v)+\sum_vt_vd(u,v)$ ,则我们只需要知道对于每一个节点$\sum_{v}d(u,v)$和$\sum_vt_vd(u,v)$的值即可。直接计算复杂度太高,我们可以考虑通过一个已经计算出值的节点计算当前节点值。

当我们已经计算得到$u$的值,并想要得到与$u$相邻$u’$的值时,可以将除了$u’$的其他节点分为两部分,一部分不经过$u$,另一部分经过$u$。

  • 计算$\sum_{v}d(u,v)$时,对于这两部分节点来说,不经过$u$的节点的$d(u’,v)$相比$d(u,v)$减少了$w_i$,经过$u$的节点的$d(u’,v)$相比$d(u,v)$增加了$w_i$,故只需知道两部分节点的分别数量即可得到$\sum_{v}d(u,v)$。

  • 计算$\sum_vt_vd(u,v)$时,对于这两部分节点来说,不经过$u$的节点的$t_vd(u’,v)$相比$t_vd(u,v)$减少了$w_i\sum_vtv$,经过$u$的节点的$t_vd(u’,v)$相比$t_vd(u,v)$增加了了$w_i\sum_vtv$,所以这部分只需知道两部分节点的点值和即可得到$\sum_vt_vd(u,v)$

有了转移方法后,我们需要获得对于每一个节点两部分分别的值,故可以指定任意节点为根节点,通过一次dfs,得到子节点数量和及子节点点值和后,通过所有节点数量和所有节点点值和得到另一部分节点的对应值,即可进行转移。

代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int N = 5e5 + 10;
long long sum1[N],sum2[N],num[N],sum[N],Sum;
struct node{
int u,v,w,next;
}a[N];
int n;
int f[N],cnt;
long long value[N],dis[N];
void add(int u,int v,int w){
a[++cnt].v =v;
a[cnt].w = w;
a[cnt].next = f[u];
f[u] = cnt;
}
int vis[N];

void dfs1(int u){
num[u] = 1;
sum[u] = value[u];
for(int i = f[u];i;i = a[i].next){
int v = a[i].v,w = a[i].w;
if(!vis[v]){
vis[v] = 1;
dis[v] = dis[u]+w;
sum1[1] += dis[v];
sum2[1] += dis[v]*value[v];
dfs1(v);
num[u] += num[v];
sum[u] += sum[v];
}
}
}
long long ans[N];
void dfs2(int u){
ans[u] = sum1[u]*value[u] + sum2[u];
for(int i = f[u];i;i = a[i].next){
int v = a[i].v,w = a[i].w;
if(!vis[v]){
vis[v] = 1;
sum1[v] = sum1[u] + (n-num[v]*2)*w;
sum2[v] = sum2[u] + (Sum-sum[v]*2)*w;
dfs2(v);
}
}
}
int main() {
cin>>n;
for(int i = 1;i<=n;i++){
scanf("%d",&value[i]);
Sum += value[i];
}
int u,v,w;
for(int i = 1;i<n;i++){
scanf("%d %d %d",&u,&v,&w);
add(u,v,w);
add(v,u,w);
}
memset(vis,0,sizeof vis);
vis[1] = 1;
dfs1(1);
memset(vis,0,sizeof vis);
vis[1] = 1;
dfs2(1);
for(int i = 1;i<=n;i++){
printf("%lld\n",ans[i]);
}
return 0;
}