CF771C Bear and Tree Jumps（换根dp）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=4e5+10;
int h[N],ne[N],e[N],idx;
ll f[N][6];
ll g[N][6];
int n,k;
ll ans=0;
void add(int a,int b){
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void dfs(int u,int fa){
    int i;
    f[u][k]=0,g[u][k]=1;
    for(i=h[u];i!=-1;i=ne[i]){
        int j=e[i];
        if(j==fa)
            continue;
        dfs(j,u);
        for(int x=1;x<=k;x++){
            if(x==1){
                f[u][x]+=f[j][k]+g[j][k];
                g[u][x]+=g[j][k];
            }
            else{
                f[u][x]+=f[j][x-1];
                g[u][x]+=g[j][x-1];
            }
        }
    }
}
void dfs1(int u,int fa){
    for (int i=1;i<=k;i++)
        ans=ans+1ll*f[u][i];
    for (int i=h[u];i!=-1;i=ne[i]){
        int x=e[i];
        if(x==fa)
            continue;
        ll tmpa[10],tmpb[10],tmpc[10],tmpd[10];
        for(int j=1;j<=k;j++) tmpa[j]=f[u][j];
        for(int j=1;j<=k;j++) tmpc[j]=g[u][j];
        for(int j=1;j<=k;j++)
            if(j==1)
                f[u][j]=f[u][j]-f[x][k]-g[x][k],g[u][j]=g[u][j]-g[x][k];
            else
                f[u][j]-=f[x][j-1],g[u][j]-=g[x][j-1];
        for(int j=1;j<=k;j++)
            if(j==1) f[x][j]+=f[u][k]+g[u][k],g[x][j]+=g[u][k];
            else f[x][j]+=f[u][j-1],g[x][j]+=g[u][j-1];
        dfs1(x,u);
        for (int j=1;j<=k;j++) f[u][j]=tmpa[j];
        for (int j=1;j<=k;j++) g[u][j]=tmpc[j];
    }
}
int main(){
    ios::sync_with_stdio(false);
    memset(h,-1,sizeof h);
    int i;
    cin>>n>>k;
    for(i=1;i<n;i++){
        int a,b;
        cin>>a>>b;
        add(a,b);
        add(b,a);
    }
    dfs(1,-1);
    dfs1(1,-1);
    cout<<ans/2<<endl;
    return 0;
}