【CF860E】Arkady and a Nobody-men
题意：给你一棵n个点的有根树。如果b是a的祖先，定义$r(a,b)$为b的子树中深度小于等于a的深度的点的个数（包括a）。定义$z(a)=\sum\limits r(a,b)$（b是a的祖先）。要你求出每个点的z值。
$n\le 5\times 10^5$
题解：一开始naive的思路：将所有点按深度排序，将深度相同的点统一处理，统计答案时相当于链加，链求和，用树剖+树状数组搞一搞，时间复杂度$O(n\log^2n)$。
后来看题解发现我这个想法简直菜爆了。我们先从树形DP的角度去想，先给出转移方程：
$ans(x)=ans(fa(x))+dep(x)+ans‘(x)$，ans‘(x)表示与a深度相同的点 对a的贡献。
现在问题变成了求ans‘，我们考虑在每个点对的lca处统计贡献。具体地，我们对于每个点x，维护若干个三元组(d,a,cnt)表示x的子树中有cnt个d级子孙，其中一个子孙为a。DP的过程就相当于在父亲节点处将所有儿子节点的三元组合并，在合并时顺便统计贡献。
具体地，合并方式如下：假如x有两个儿子，它们有三元组$(d,a,cnt_a)$和$(d,b,cnt_b)$，则：
1.ans‘(a)+=dep(x)\times cnt_b
2.ans‘(b)+=dep(x)\times cnt_a
3.得到新三元组(d,a,cnt_a+cnt_b)
但是后面的点 对b的贡献呢？我们发现后面的点 对a和b的贡献就是相同的了，所以我们建一个新图，在新图中从a到b连一条长度为$ans‘(b)-ans‘(a)$的边，最后在新图上DFS一下，最最后统计一下ans数组即可。
以上过程采用长链剖分优化，由于一开始的三元组个数为n，则每次合并都会减少一个三元组，所以时间复杂度O(n)。
#include <cstdio> #include <cstring> #include <iostream> #include <vector> using namespace std; const int maxn=500010; typedef long long ll; ll ans[maxn]; struct node { int v,x; node() {} node(int a,int b) {v=a,x=b;} }mem[maxn<<1],*f[maxn],*now=mem; int n,cnt,Cnt,rt; int to[maxn],nxt[maxn],head[maxn],dep[maxn],md[maxn],son[maxn],fa[maxn]; bool vis[maxn]; int To[maxn],Nxt[maxn],Head[maxn]; ll Val[maxn]; inline void add(int a,int b) { to[cnt]=b,nxt[cnt]=head[a],head[a]=cnt++; } inline void Add(int a,int b,ll c) { To[Cnt]=b,Val[Cnt]=c,Nxt[Cnt]=Head[a],Head[a]=Cnt++; } void dfs1(int x) { md[x]=0; for(int i=head[x];i!=-1;i=nxt[i]) { dep[to[i]]=dep[x]+1,dfs1(to[i]); if(md[to[i]]+1>md[x]) md[x]=md[to[i]]+1,son[x]=to[i]; } } void dfs2(int x) { if(f[x]==NULL) f[x]=now,now+=md[x]+2; if(son[x]) f[son[x]]=f[x]+1,dfs2(son[x]); f[x][0]=node(1,x); for(int i=head[x];i!=-1;i=nxt[i]) if(to[i]!=son[x]) { dfs2(to[i]); for(int j=0;j<=md[to[i]];j++) { node a=f[x][j+1],b=f[to[i]][j]; ans[b.x]+=1ll*dep[x]*a.v; ans[a.x]+=1ll*dep[x]*b.v; Add(a.x,b.x,ans[b.x]-ans[a.x]); f[x][j+1]=node(a.v+b.v,a.x); } } } void dfs3(int x) { for(int i=Head[x];i!=-1;i=Nxt[i]) ans[To[i]]=ans[x]+Val[i],dfs3(To[i]); } void dfs4(int x) { for(int i=head[x];i!=-1;i=nxt[i]) ans[to[i]]+=ans[x]+dep[x],dfs4(to[i]); } inline int rd() { int ret=0,f=1; char gc=getchar(); while(gc<‘0‘||gc>‘9‘) {if(gc==‘-‘) f=-f; gc=getchar();} while(gc>=‘0‘&&gc<=‘9‘) ret=ret*10+(gc^‘0‘),gc=getchar(); return ret*f; } int main() { n=rd(); int i; memset(head,-1,sizeof(head)),memset(Head,-1,sizeof(Head)); for(i=1;i<=n;i++) { fa[i]=rd(); if(!fa[i]) rt=i; else add(fa[i],i); } dep[rt]=1,dfs1(rt); dfs2(rt); for(i=0;i<=md[rt];i++) dfs3(f[rt][i].x); dfs4(rt); for(i=1;i<=n;i++) printf("%lld ",ans[i]); return 0; }