【20160815】noip模拟（未完）

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int N=100010;
int n,m,len;
struct node{
    int x,y,d,next;
}a[2*N];
int first[N],f0[N],f1[N][20],g0[N],g1[N][20],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][0]=1;//debug
    tot[x]=1;
    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%16));
            for(int j=0;j<=15;j++)
            {
                int t=j+now%16;
                f1[x][t%16]+=f1[y][j];
                f0[x]+=(t-(t%16))*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=0,c[20];
            memset(c,0,sizeof(c));
            int now=a[i].d;
            z+=f0[y]+tot[y]*(now-(now%16));
            for(int j=0;j<=15;j++)
            {
                int t=j+now%16;
                c[t%16]+=f1[y][j];//以y为根的子树中每个点都要加d(x,y),j变为t%16
                z+=(t-(t%16))*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=0;j<=15;j++) g1[y][j]+=f1[x][j]+g1[x][j]-c[j];
            
            //g[y]+=(n-tot[y])*a[i].d
            for(int j=0;j<=15;j++) c[j]=g1[y][j];
            for(int j=0;j<=15;j++) 
            {
                int t=j+now;
                g1[y][t%16]=c[j];
                g0[y]+=(t-(t%16))*c[j];
            }
            dfs2(y,x);
        }
    }
}

int main()
{
    // freopen("a.in","r",stdin);
    freopen("warehouse.in","r",stdin);
    freopen("warehouse.out","w",stdout);
    len=0;
    memset(first,0,sizeof(first));
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int x,y,d;
        scanf("%d%d%d",&x,&y,&d);
        ins(x,y,d),ins(y,x,d);
    }
    memset(f0,0,sizeof(f0));
    memset(f1,0,sizeof(f1));
    memset(g0,0,sizeof(g0));
    memset(g1,0,sizeof(g1));
    dfs1(1,0);dfs2(1,0);
    int ans=0;
    for(int i=1;i<=n;i++) 
    {
        ans=f0[i]+g0[i];
        for(int j=0;j<=15;j++)
        {
            int now=j^m;
            ans+=(f1[i][j]+g1[i][j])*now;
        }
        ans-=m;
        printf("%d
",ans);
    }
    return 0;
}