#include <bits/stdc++.h>
using namespace std;
long long cnt[200100][5], tot[200100], ans = 0;
int k, n;
vector<int> g[200100];
void dfs(int u, int p, int d) {
  int dmodk, res;
  tot[u] = cnt[u][d % k] = 1;
  for (int i = 0; i < g[u].size(); i++) {
    if (g[u][i] != p) {
      dfs(g[u][i], u, d + 1);
      for (int j = 0; j < k; j++) {
        for (int l = 0; l < k; l++) {
          dmodk = ((j + l - 2 * d) % k + k) % k;
          res = (((k - dmodk) % k) + k) % k;
          ans += 1ll * res * cnt[u][j] * cnt[g[u][i]][l];
        }
      }
      for (int j = 0; j < k; j++) cnt[u][j] += cnt[g[u][i]][j];
      tot[u] += tot[g[u][i]];
    }
  }
  ans += tot[u] * (n - tot[u]);
}
int main() {
  int i, u, v;
  scanf("%d%d", &n, &k);
  for (i = 0; i < n - 1; i++) {
    scanf("%d%d", &u, &v);
    g[u].push_back(v);
    g[v].push_back(u);
  }
  dfs(1, 0, 0);
  printf("%I64d\n", ans / k);
}
