結果

問題 No.1769 Don't Stop the Game
ユーザー polylogKpolylogK
提出日時 2021-11-03 13:40:51
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 381 ms / 3,000 ms
コード長 3,028 bytes
コンパイル時間 3,295 ms
コンパイル使用メモリ 239,356 KB
実行使用メモリ 57,180 KB
最終ジャッジ日時 2024-06-29 18:02:18
合計ジャッジ時間 10,914 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 AC 1 ms
6,944 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 1 ms
6,940 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 3 ms
6,940 KB
testcase_08 AC 226 ms
16,268 KB
testcase_09 AC 160 ms
12,556 KB
testcase_10 AC 235 ms
16,224 KB
testcase_11 AC 153 ms
12,800 KB
testcase_12 AC 256 ms
19,100 KB
testcase_13 AC 264 ms
19,392 KB
testcase_14 AC 269 ms
19,292 KB
testcase_15 AC 274 ms
19,380 KB
testcase_16 AC 284 ms
19,276 KB
testcase_17 AC 318 ms
19,532 KB
testcase_18 AC 334 ms
20,460 KB
testcase_19 AC 347 ms
20,880 KB
testcase_20 AC 367 ms
20,864 KB
testcase_21 AC 373 ms
21,012 KB
testcase_22 AC 381 ms
20,804 KB
testcase_23 AC 261 ms
19,376 KB
testcase_24 AC 286 ms
19,288 KB
testcase_25 AC 145 ms
19,744 KB
testcase_26 AC 216 ms
21,156 KB
testcase_27 AC 197 ms
55,672 KB
testcase_28 AC 272 ms
57,180 KB
testcase_29 AC 263 ms
39,964 KB
testcase_30 AC 263 ms
40,012 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<n;++i)

int main() {
    int n; cin >> n;
    vector<vector<pair<int, int>>> g(n);
    rep(i, n - 1) {
        int a, b, c; cin >> a >> b >> c;
        --a, --b;
        g[a].emplace_back(b, c);
        g[b].emplace_back(a, c);
    }

    vector<int> siz(n, 1), x(n);
    auto dfs0 = [&](auto&& f, int cur, int par)->void {
        int par_idx = -1;
        rep(i, g[cur].size()) {
            auto [to, weight] = g[cur][i];
            if (to == par) {
                par_idx = i;
                continue;
            }
            x[to] = x[cur] ^ weight;
            f(f, to, cur);
            siz[cur] += siz[to];
        }
        if (par_idx != -1) {
            swap(g[cur][par_idx], g[cur].back());
            g[cur].pop_back();
        }
    };
    dfs0(dfs0, 0, -1);

    auto sorted_x = x;
    sort(sorted_x.begin(), sorted_x.end());
    sorted_x.erase(unique(sorted_x.begin(), sorted_x.end()), sorted_x.end());
    rep(i, n)x[i] = lower_bound(sorted_x.begin(), sorted_x.end(), x[i]) - sorted_x.begin();//compress

    vector<int> y_sum(sorted_x.size()), cnt_sum(sorted_x.size());
    vector<int> y(n, -1), cnt(n, -1);
    ll ans = ll(n) * (n - 1);
    int parity[2] = { 1,-1 };
    auto dfs1 = [&](auto&& f, int cur, bool type)->void {
        y_sum[x[cur]] += parity[type] * y[cur];
        cnt_sum[x[cur]] += parity[type] * cnt[cur];
        for (int i = 0; i < g[cur].size(); ++i) {
            f(f, g[cur][i].first, type);
        }
    };
    auto dfs2 = [&](auto&& f, int cur, int par, bool type)->void {
        sort(g[cur].begin(), g[cur].end(), [&](const pair<int, int>& i, const pair<int, int>& j) {return siz[i.first] > siz[j.first]; });
        for (int i = 1; i < g[cur].size(); ++i) {
            f(f, g[cur][i].first, cur, true);
        }
        if (!g[cur].empty())f(f, g[cur][0].first, cur, false);
        for (int i = 1; i < g[cur].size(); ++i) {
            dfs1(dfs1, g[cur][i].first, false);
        }

        if (y[cur] == -1) {
            y[cur] = siz[cur] - y_sum[x[cur]];
            cnt[cur] = 1 - cnt_sum[x[cur]];
        }
        ans -= y_sum[x[cur]];//case-1
        y_sum[x[cur]] += y[cur];
        cnt_sum[x[cur]] += cnt[cur];

        if (cur == 0) {
            rep(i, sorted_x.size())if (i != x[cur]) {
                ans -= ll(y_sum[i]) * cnt_sum[i];//case-3
            }
        }
        else {
            ans -= ll(cnt_sum[x[par]]) * (n - siz[cur]);//case-2
            ans -= ll(y_sum[x[par]]) * cnt_sum[x[par]];//case-3
            ans += siz[cur];//case-3
        }

        if (type) {
            for (int i = 0; i < g[cur].size(); ++i) {
                dfs1(dfs1, g[cur][i].first, true);
            }
            y_sum[x[cur]] -= y[cur];
            cnt_sum[x[cur]] -= cnt[cur];
        }
    };
    dfs2(dfs2, 0, -1, false);
    cout << ans << "\n";

    return 0;
}
0