【20160815】noip模拟（未完）

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 using namespace std;
 
 const int N=;
 int n,m,len;
 struct node{
     int x,y,d,next;
 }a[*N];
 int first[N],f0[N],f1[N][],g0[N],g1[N][],tot[N];
 /*
 若m=0,
 f[x]：x的子树的节点到x的距离的总和
 g[x]：不是x的子树上的点到x的距离的总和
 f[x]=f[y]+d(x,y)*tot[y] y=all son[x]
 g[x]=f[fa]+g[fa] - (f[x]+tot[x]*d(x,fa)) + (n-tot[x])*d(x,fa)
 
 因为m<=15，异或只会影响最后四位数
 则将f[x]拆分为f0[x]：f[x]除去最后四位数的部分：xx0000
               f1[x][0~15]：f1[x][j]表示f[x]中有多少个最后四位为j的距离（不能像原来一样直接加进去，因为最后要异或）
 g[x]同理。
 */
 
 void ins(int x,int y,int d)
 {
     len++;
     a[len].x=x;a[len].y=y;a[len].d=d;
     a[len].next=first[x];first[x]=len;
 }
 
 void dfs1(int x,int fa)
 {
     f1[x][]=;//debug
     tot[x]=;
     for(int i=first[x];i;i=a[i].next)
     {
         int y=a[i].y;
         if(y!=fa)
         {
             dfs1(y,x);
             tot[x]+=tot[y];
             // f[x]+=f[y]+a[i].d*tot[y];
             int now=a[i].d;
             f0[x]+=f0[y]+tot[y]*(now-(now%));
             for(int j=;j<=;j++)
             {
                 int t=j+now%;
                 f1[x][t%]+=f1[y][j];
                 f0[x]+=(t-(t%))*f1[y][j];
             }
         }
     }
 }
 
 void dfs2(int x,int fa)
 {
     for(int i=first[x];i;i=a[i].next)
     {
         int y=a[i].y;
         if(y!=fa)
         {
             // g[y]=f[x]+g[x]-f[y]+(n-2*tot[y])*a[i].d;
             // g[y]=f[x]+g[x] - (f[y]+tot[y]*a[i].d) + (n-tot[y])*a[i].d
 
             /*
             减法不会出现不够减：
             因为减的东西之前一定加过（所以要考虑它的实际意义）
             */
 
             //z c  f[y]+tot[y]*a[i].d
             int z=,c[];
             memset(c,,sizeof(c));
             int now=a[i].d;
             z+=f0[y]+tot[y]*(now-(now%));
             for(int j=;j<=;j++)
             {
                 int t=j+now%;
                 c[t%]+=f1[y][j];//以y为根的子树中每个点都要加d(x,y),j变为t%16
                 z+=(t-(t%))*f1[y][j];
             }
 
             //g[y]=f[x]+g[x]-(f[y]+tot[y]*a[i].d)
             g0[y]+=f0[x]+g0[x]-z;
             for(int j=;j<=;j++) g1[y][j]+=f1[x][j]+g1[x][j]-c[j];
 
             //g[y]+=(n-tot[y])*a[i].d
             for(int j=;j<=;j++) c[j]=g1[y][j];
             for(int j=;j<=;j++)
             {
                 int t=j+now;
                 g1[y][t%]=c[j];
                 g0[y]+=(t-(t%))*c[j];
             }
             dfs2(y,x);
         }
     }
 }
 
 int main()
 {
     // freopen("a.in","r",stdin);
     freopen("warehouse.in","r",stdin);
     freopen("warehouse.out","w",stdout);
     len=;
     memset(first,,sizeof(first));
     scanf("%d%d",&n,&m);
     for(int i=;i<n;i++)
     {
         int x,y,d;
         scanf("%d%d%d",&x,&y,&d);
         ins(x,y,d),ins(y,x,d);
     }
     memset(f0,,sizeof(f0));
     memset(f1,,sizeof(f1));
     memset(g0,,sizeof(g0));
     memset(g1,,sizeof(g1));
     dfs1(,);dfs2(,);
     int ans=;
     for(int i=;i<=n;i++)
     {
         ans=f0[i]+g0[i];
         for(int j=;j<=;j++)
         {
             int now=j^m;
             ans+=(f1[i][j]+g1[i][j])*now;
         }
         ans-=m;
         printf("%d\n",ans);
     }
     return ;
 }