小奇的仓库：换根dp

一道很好的换根dp题。考场上现场yy十分愉快

给定树，求每个点的到其它所有点的距离异或上m之后的值，n=100000,m<=16

只能线性复杂度求解，m又小得奇怪。或者带一个log像kx一样打一个线段树

我们可以发现，m小的话对距离很大的路径的影响也不会超过16。

那么变化的其实就是最后4个二进制位啊。

所以我们像普通的换根dp一样求出所有距离，在额外处理一下以p为端点的全部路径里路径长度%16之后的值为k的有多少个

设为bits2[k][p]

因为换根dp的主要思路是两遍dfs，第一次处理每个点的子树，第二遍处理以这个点为端点的所有点对的状态

设子树里所有的点到当前根节点的距离是dt[p]，子树所有点到当前根节点的路径%16之后值为k的路径有k个设为bits[k][p]

那么第一轮dfs转移还是很显然的，我稍微说一下第二轮。

dt[p]=dt[p]+dt[f[p]]-(dt[p]+siz[p]*vtf[p])+vtf[p]*(n-siz[p]);其中vtf[p]表示p节点到其父亲的距离

=p的子树里所有点到p的距离+所有点到p的父亲的距离-（p子树内的点到p的距离+这些点从p再到p父亲的距离）+子树p外的所有点从p父亲到p的距离

两个公式的每个部分都是对应的

那么bits数组的转移也是大同小异了（虽说稍复杂一些）

for(int j=;j<=;++j) bits2[j+vtf[p]&][p]=bits[j+vtf[p]&][p]+bits2[j][f[p]]-bits[j-vtf[p]+&][p];
首先是子树里本身就有bits[j+vtf[p]&]个点满足条件
到父亲的距离为j的所有点里去除子树p里的贡献就是bits2[j][f[p]]-bits[j-vtf[p]+&][p]

自己想？

容易重复考虑自己到自己的距离为0而异或后不是0的情况，要减去m

具体实现细节自己想？有什么不会的评论区吧，我扔个代码就跑～

 //这个异或是个狗屎啊！！！但是其实异或的二进制位也不超过4,特殊处理一下？
 //对于每一个节点计算出子树所有点到它的距离和，再计算出有多少点到它的距离&15后是几
 #include<cstdio>
 #define int long long
 int n,m,cnt,fir[],l[],to[],v[],dt[],siz[];
 int bits[][],f[],vtf[],bits2[][];
 void dfs(int p,int fa){
     siz[p]=; bits[][p]=; f[p]=fa;
     for(int i=fir[p];i;i=l[i])if(to[i]!=fa){
         dfs(to[i],p); siz[p]+=siz[to[i]]; dt[p]+=dt[to[i]]+siz[to[i]]*v[i]; vtf[to[i]]=v[i];
         for(int j=;j<=;++j) bits[j+v[i]&][p]+=bits[j][to[i]];
     }
 }
 void DFS(int p){
     if(p!=){
         dt[p]+=dt[f[p]]-(dt[p]+siz[p]*vtf[p])+vtf[p]*(n-siz[p]);
         for(int j=;j<=;++j) bits2[j+vtf[p]&][p]=bits[j+vtf[p]&][p]+bits2[j][f[p]]-bits[j-vtf[p]+&][p];
     }else for(int j=;j<=;++j)bits2[j][]=bits[j][];
     for(int i=fir[p];i;i=l[i]) if(to[i]!=f[p]) DFS(to[i]);
 }
 void link(int a,int b,int vv){l[++cnt]=fir[a];fir[a]=cnt;to[cnt]=b;v[cnt]=vv;}
 signed main(){//freopen("newt3.in","r",stdin);freopen("newcot3.out","w",stdout);
     scanf("%lld%lld",&n,&m);
     for(int i=,x,y,ww;i<n;++i)scanf("%lld%lld%lld",&x,&y,&ww),link(x,y,ww),link(y,x,ww);
     dfs(,);DFS();//for(int i=1;i<=n;++i)printf("%lld\n",dt[i]);
     for(int i=;i<=n;++i){int ans=dt[i];for(int j=;j<=;++j)ans+=bits2[j][i]*((j^m)-j);printf("%lld\n",ans-m);}
 }