#include <bits/stdc++.h>
using namespace std;
long long dp[200005], num[200005][5];
int n, k;
vector<int> g[200005];
void dfs1(int x, int fa) {
  if (g[x].size() == 1 && g[x][0] == fa) {
    return;
  }
  for (int i = 0; i < g[x].size(); i++) {
    int t = g[x][i];
    if (t == fa) continue;
    dfs1(t, x);
    for (int j = 0; j < k; j++) {
      num[x][j] += num[t][(j - 1 + k) % k];
    }
    num[x][1 % k]++;
    dp[x] += dp[t] + 1 + num[t][0];
  }
}
void dfs2(int x, int fa) {
  for (int i = 0; i < g[x].size(); i++) {
    int t = g[x][i];
    if (t == fa) continue;
    for (int j = 0; j < k; j++) {
      num[x][j] -= num[t][(j - 1 + k) % k];
    }
    num[x][1 % k]--;
    dp[t] += dp[x] - dp[t] - 1 - num[t][0];
    dp[t] += 1 + num[x][0];
    for (int j = 0; j < k; j++) {
      num[t][j] += num[x][(j - 1 + k) % k];
    }
    num[t][1 % k]++;
    dfs2(t, x);
    num[t][1 % k]--;
    for (int j = 0; j < k; j++) {
      num[t][j] -= num[x][(j - 1 + k) % k];
    }
    num[x][1 % k]++;
    for (int j = 0; j < k; j++) {
      num[x][j] += num[t][(j - 1 + k) % k];
    }
  }
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  cin >> n >> k;
  for (int i = 1; i < n; i++) {
    int x, y;
    cin >> x >> y;
    g[x].push_back(y);
    g[y].push_back(x);
  }
  dfs1(1, 0);
  dfs2(1, 0);
  long long ans = 0;
  for (int i = 1; i <= n; i++) {
    ans += dp[i];
  }
  cout << ans / 2 << '\n';
  return 0;
}
