scutcode123. bx和妹子的游戏

xiaoxiao2021-02-27  790

123. bx和妹子的游戏

题目描述

这天,bxbx和妹子在一个有根树上玩游戏。

这个有根树有nn个节点,编号从11nn,其中11号节点是根节点。除了11号节点每个节点都有一个父亲节点。第ii个节点有权值a_iai

定义s(u, v)s(u,v)uuvv的简单路径上所有点的权值之和(包括uuvv)。

bxbx会随机选一个节点uu,妹子随机选一个节点vv,假设uuvvLCALCAppbxbx能获得\max(s(u, p), s(v, p))max(s(u,p),s(v,p))的愉悦值。

现在你需要回答,考虑所有情况之下(n^2n2种情况),bxbx所获得的愉悦值之和。

一个节点uu的祖先定义为从uu一直往父亲节点走若干次能到达的节点,其中uu是自身的祖先。两个节点的LCALCA(最近公共祖先)定义为既是uu的祖先,也是vv的祖先,并且和uu的距离最近的一个节点。

输入格式

输入第一行包括一个整数TT,表示数据组数。

对于每组数据,第一行一个整数nn

第二行n-1n1个整数,第ii个整数表示节点i+1i+1的父亲p_{i+1}pi+1

第三行nn个整数,第ii个整数表示节点i的权值a_iai

1 \le T \le 201T20

2 \le n \le 1000002n100000

1 \le p_{i + 1} \le i, 1 \leq i \leq n - 11pi+1i,1in1

1 \leq a_i \leq 100001ai10000

输出格式

对每组数据,输出 bx 获得的愉悦值之和。

样例数据

输入

2 3 1 1 1 2 3 2 1 1 2

输出

28 9

max( S[u,p] , S[v,p] ) = max( S[1,u] , S[1,v] ) - ( S[1,p] - a[p] )

显然可以在O(n)内求出S[1,k]

S[1,i]对答案的贡献=

S[1,p]-a[p] 对答案的贡献为 = - ( p 作为 lca 的次数 * (S[1,p]-a[p]) )

#include<iostream> #include<cstdlib> #include<cstdio> #include<string> #include<vector> #include<deque> #include<queue> #include<algorithm> #include<set> #include<map> #include<stack> //#include<cti1me> #include<string.h> #include<math.h> #include<list> using namespace std; #define ll long long #define pii pair<int,int> const int inf = 1e9+7; const int N = 100000 + 5; struct Edge{ int to,next; }edge[N]; int head[N]; int a[N]; void addEdge(int k,int u,int v){ edge[k].to=v; edge[k].next=head[u]; head[u]=k; } ll s[N]; ll childs[N]; ll tS[N]; void dfs(int u,int sum){ s[u]=sum+a[u]; childs[u]=1; for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; dfs(v,s[u]); childs[u]+=childs[v]; } } ll f(int u){//u作为lca有多少种情况 ll ans=1; for(int i=head[u];i!=-1;i=edge[i].next){ int x=edge[i].to; ans+=childs[x]; for(int j=edge[i].next;j!=-1;j=edge[j].next){ int y=edge[j].to; ans+=childs[x]*childs[y]; } } return ans; } ll slove(int n){ dfs(1,0); ll ans=0; for(int i=1;i<=n;++i){ tS[i]=s[i]; } sort(s+1,s+n+1); for(int i=1;i<=n;++i){ ans+=s[i]*i; ans-=f(i)*(tS[i]-a[i]); } return ans; } int main() { //freopen("/home/lu/Documents/r.txt","r",stdin); //freopen("/home/lu/Documents/w.txt","w",stdout); int T; scanf("%d",&T); while(T--){ int n; scanf("%d",&n); fill(head,head+n+1,-1); int nume=0; for(int i=2;i<=n;++i){ int p; scanf("%d",&p); addEdge(nume++,p,i); } ll sub=0; for(int i=1;i<=n;++i){ scanf("%d",&a[i]); sub+=a[i]; } printf("%lld\n",2*slove(n)-sub); } return 0; }

转载请注明原文地址: https://www.6miu.com/read-165.html

最新回复(0)