結果

問題 No.1002 Twotone
ユーザー keymoonkeymoon
提出日時 2024-02-23 03:56:54
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,224 ms / 5,000 ms
コード長 2,856 bytes
コンパイル時間 3,931 ms
コンパイル使用メモリ 277,316 KB
実行使用メモリ 71,956 KB
最終ジャッジ日時 2024-09-29 04:46:28
合計ジャッジ時間 20,530 ms
ジャッジサーバーID
(参考情報)
judge2 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 2 ms
6,820 KB
testcase_03 AC 459 ms
37,104 KB
testcase_04 AC 633 ms
47,352 KB
testcase_05 AC 633 ms
47,356 KB
testcase_06 AC 2 ms
6,820 KB
testcase_07 AC 316 ms
30,116 KB
testcase_08 AC 552 ms
47,368 KB
testcase_09 AC 576 ms
47,232 KB
testcase_10 AC 2 ms
6,820 KB
testcase_11 AC 632 ms
37,936 KB
testcase_12 AC 886 ms
48,020 KB
testcase_13 AC 880 ms
48,020 KB
testcase_14 AC 2 ms
6,820 KB
testcase_15 AC 387 ms
28,860 KB
testcase_16 AC 851 ms
47,760 KB
testcase_17 AC 833 ms
47,760 KB
testcase_18 AC 2 ms
6,816 KB
testcase_19 AC 661 ms
41,296 KB
testcase_20 AC 774 ms
47,224 KB
testcase_21 AC 800 ms
47,352 KB
testcase_22 AC 3 ms
6,820 KB
testcase_23 AC 652 ms
40,944 KB
testcase_24 AC 1,224 ms
64,252 KB
testcase_25 AC 1,221 ms
64,508 KB
testcase_26 AC 2 ms
6,816 KB
testcase_27 AC 123 ms
33,492 KB
testcase_28 AC 243 ms
47,872 KB
testcase_29 AC 199 ms
47,864 KB
testcase_30 AC 2 ms
6,816 KB
testcase_31 AC 185 ms
47,556 KB
testcase_32 AC 273 ms
47,796 KB
testcase_33 AC 194 ms
47,836 KB
testcase_34 AC 991 ms
71,956 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using ll = long long;

struct CentroidDecomposition {
  int n, r;
  vector<vector<int>> g, cd;
  vector<int> ord, szs;
  CentroidDecomposition(vector<vector<int>> &gg): g(gg), n(gg.size()), cd(n), szs(n,-1) {
    calc(0); dfs(0,-1,n,-1);
  }
  void calc(int v) {
    stack<int> s; s.emplace(v);
    while (s.size()) {
      int i = s.top(); s.pop();
      if (i<0) szs[~i]=1;
      else { szs[i]=0; s.emplace(~i); }
      for(auto e:g[i<0?~i:i]) {
        if (!szs[e]) continue;
        if (i<0) szs[~i] += szs[e];
        else s.emplace(e);
      }
    }
  }
  void dfs(int v, int p, int sz, int root){
    if (szs[v]==0) return;
    begin:
    for(auto e:g[v]) {
      if (p==e or szs[e]<sz-sz/2) continue;
      p=v; v=e; goto begin;
    }
    if(root!=-1)cd[root].emplace_back(v);
    else r=v;
    ord.push_back(v);
    calc(v); szs[v]=0;
    for(auto e:g[v])dfs(e,v,szs[e],v);
    return;
  }
};

int main() {
  cin.tie(0)->sync_with_stdio(false);
  int n, k;
  cin >> n >> k;
  using pii = pair<int, int>;
  vector<vector<int>> g(n);
  vector<vector<pii>> edges(n);
  for (int i = 0; i < n - 1; i++) {
    int u, v, c;
    cin >> u >> v >> c;
    u--, v--, c--;
    g[u].push_back(v);
    g[v].push_back(u);
    edges[u].emplace_back(v, c);
    edges[v].emplace_back(u, c);
  }
  CentroidDecomposition cd(g);
  vector<int> visited(n);
  ll ans = 0;
  for (int ce : cd.ord) {
    int cnt0 = 0; // *, -
    int cnt1 = 0; // *, *
    map<pii, int> cnt2; // x, y
    map<int, int> cnt3; // x, *
    map<int, int> cnt4; // x, -
    auto f1 = [&](auto self, int node, int parent, int c1, int c2) -> void {
      if (c2 == -1) {
        cnt0++;
        cnt4[c1]++;
      } else {
        cnt1++;
        cnt2[make_pair(min(c1, c2), max(c1, c2))]++;
        cnt3[c1]++;
        cnt3[c2]++;
      }
      for (auto [child, c] : edges[node]) {
        if (child == parent || visited[child]) continue;
        if (c1 == c) {
          self(self, child, node, c1, c2);
        } else if (c2 == c || c2 == -1) {
          self(self, child, node, c1, c);
        }
      }
    };
    auto f2 = [&](auto self, int node, int parent, int c1, int c2) -> void {
      if (c2 == -1) {
        ans += cnt0 - cnt4[c1] + cnt3[c1];
      } else {
        ans += cnt2[make_pair(min(c1, c2), max(c1, c2))] + cnt4[c1] + cnt4[c2];
      }
      for (auto [child, c] : edges[node]) {
        if (child == parent || visited[child]) continue;
        if (c1 == c) {
          self(self, child, node, c1, c2);
        } else if (c2 == c || c2 == -1) {
          self(self, child, node, c1, c);
        }
      }
    };
    for (auto [child, c] : edges[ce]) {
      if (visited[child]) continue;
      f2(f2, child, ce, c, -1);
      f1(f1, child, ce, c, -1);
    }
    ans += cnt1;
    visited[ce] = true;
  }
  cout << ans << endl;
}
0