NC201400 树学题解
一.题解
这道题又是一道换根dp板子题,代码结构与 Accumulation Degree 这道题基本一致,唯一不同的就是转移了【不过转移的时候,因为方程的原因不需要特殊考虑叶节点】
我们先套路的设\(dp[i]\)表示以\(i\)为根的子树中,所有点的深度和,现在,我们来想想转移。
我们发现,如果我们要从\(i\)的一个儿子v转移到i的话,以v为根的子树中的所有节点的深度都加了1,现在我们就不能够直接求值了,怎么办呢?很简单,我们发现,既然以v为根的子树中的所有节点的深度都加了1,那么,一共增加的值,不就正好是以v为根的子树中节点的个数嘛?
所以,我们只需要再维护一个以i为根的子树的大小\(siz[i]\)即可
那么,我们就可以很简单的推出转移:
\(dp[i]=\)\(\sum_{v\in son(i)}(dp[v]+siz[v])\)
\(siz[i]=\)\(1(i本身)+\sum_{v\in son(i)}siz[v]\)
这样,我们一遍dfs就可以求出所有\(dp[i]\)了!
现在,我们考虑换根dp
我们假设将根\(x\)换成了它的一个儿子\(y\),那么,同样的,改变的只是\(x\)和\(y\)的子树,我们将这个"影响"修改一下即可~
首先,\(dp[i]-=(dp[v]+siz[v])\)(减去v的贡献)
然后,\(siz[i]-=siz[v],siz[v]+=siz[i]\)(修改两个点的子树大小)
最后,\(dp[v]+=(dp[i]+siz[i])\)(加上i的贡献)
为什么是这个顺序呢?因为修改\(dp[i]\)的时候,\(dp[v]\)和\(siz[v]\)必须和以前一样(因为我们是从这两个值转移过来的),然后要修改\(dp[v]\)的话,\(dp[i]\)和\(siz[i]\)又必须是修改后,没加v的贡献的值,所以,这样安排修改顺序是很好的。(当然,你非要用其他顺序的话,我也无话可说,只是要注意修改的方法有可能会不同)
最后,我们每次统计下以任意点为根时,\(dp[i]\)的值的最小值即可~
代码:
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e6+1;
struct node{
int v,nex;
}t[N<<1];
int dp[N],siz[N];
int las[N],len,ans=1e17;
inline void add(int u,int v){
t[++len]=(node){v,las[u]},las[u]=len;
}
inline void dfs1(int now,int fa){
siz[now]=1;
for(int i=las[now];i;i=t[i].nex){
int v=t[i].v;
if(v!=fa){
dfs1(v,now);
siz[now]+=siz[v];
dp[now]+=(dp[v]+siz[v]);
}
}
}
inline void dfs2(int now,int fa){
ans=min(ans,dp[now]);
for(int i=las[now];i;i=t[i].nex){
int v=t[i].v;
if(v!=fa){
int pas1=siz[now],pas2=siz[v];
int Pas1=dp[now],Pas2=dp[v];//继续偷懒
dp[now]-=(dp[v]+siz[v]);
siz[now]-=siz[v];siz[v]+=siz[now];
dp[v]+=(dp[now]+siz[now]);
dfs2(v,now);
siz[now]=pas1,siz[v]=pas2;
dp[now]=Pas1,dp[v]=Pas2;
}
}
}
signed main(){
int n;
scanf("%lld",&n);
for(int i=1;i<n;++i){
int u,v;
scanf("%lld%lld",&u,&v);
add(u,v),add(v,u);
}
dfs1(1,1);dfs2(1,1);
printf("%lld",ans);
return 0;
}
二.闲话
其实一开始,我变量都开的是int,然而,交上去后,算了下,每次的贡献最大可以达到\(n^2\)水平【链】,于是发现很可能爆int,于是,我马上改成long long,然后,改完后正准备提交。。。
恭喜通过~
我:???