3522. [POI2014]Hotel

https://hydro.ac/d/bzoj/record/6322aadfbdf9bc31d1605fa0

BZOJ/3522/3522.cpp

@@ -0,0 +1,64 @@
+#include <iostream>
+#include <algorithm>
+#include <array>
+#include <vector>
+
+using std::cin;
+using std::cout;
+const char endl = '\n';
+
+const int N = 5005;
+
+int n, max_dep;
+long long ans;
+std::vector<int> g[N];
+std::array<int, N> f1, f2, cnt;
+
+void dfs(int u, int f, int d) {
+    cnt[d]++;
+    max_dep = std::max(max_dep, d);
+
+    for (int v : g[u]) {
+        if (v == f) continue;
+
+        dfs(v, u, d + 1);
+    }
+}
+
+int main() {
+    std::ios::sync_with_stdio(false);
+    cin.tie(nullptr);
+
+    cin >> n;
+
+    for (int i = 1, u, v; i < n; i++) {
+        cin >> u >> v;
+
+        g[u].emplace_back(v);
+        g[v].emplace_back(u);
+    }
+
+    dfs(1, 0, 1);
+
+    for (int i = 1; i <= n; i++) {
+        std::fill(f1.begin(), f1.end(), 0);
+        std::fill(f2.begin(), f2.end(), 0);
+
+        for (int u : g[i]) {
+            max_dep = 0;
+            std::fill(cnt.begin(), cnt.end(), 0);
+
+            dfs(u, i, 1);
+
+            for (int j = 1; j <= max_dep; j++) {
+                ans += f1[j] * cnt[j];
+                f1[j] += f2[j] * cnt[j];
+                f2[j] += cnt[j];
+            }
+        }
+    }
+
+    cout << ans << endl;
+
+    return 0;
+}