#437247 (C++17) No.1002 Twotone

提出ソース

結果

問題	No.1002 Twotone
ユーザー	ferin
提出日時	2020-03-03 00:40:11
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	WA
実行時間	-
コード長	5,686 bytes
コンパイル時間	3,583 ms
コンパイル使用メモリ	220,688 KB
実行使用メモリ	40,992 KB
最終ジャッジ日時	2024-04-21 23:00:18
合計ジャッジ時間	25,127 ms
ジャッジサーバーID （参考情報）	judge5 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,376 KB
testcase_03	WA	-
testcase_04	WA	-
testcase_05	WA	-
testcase_06	WA	-
testcase_07	WA	-
testcase_08	WA	-
testcase_09	WA	-
testcase_10	WA	-
testcase_11	WA	-
testcase_12	WA	-
testcase_13	WA	-
testcase_14	WA	-
testcase_15	WA	-
testcase_16	WA	-
testcase_17	WA	-
testcase_18	WA	-
testcase_19	WA	-
testcase_20	WA	-
testcase_21	WA	-
testcase_22	WA	-
testcase_23	WA	-
testcase_24	WA	-
testcase_25	WA	-
testcase_26	WA	-
testcase_27	WA	-
testcase_28	WA	-
testcase_29	WA	-
testcase_30	WA	-
testcase_31	WA	-
testcase_32	WA	-
testcase_33	WA	-
testcase_34	WA	-

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using PII = pair<ll, ll>;
#define FOR(i, a, n) for (ll i = (ll)a; i < (ll)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
template<typename T> void chmin(T &a, const T &b) { a = min(a, b); }
template<typename T> void chmax(T &a, const T &b) { a = max(a, b); }
struct FastIO {FastIO() { cin.tie(0); ios::sync_with_stdio(0); }}fastiofastio;
#ifdef DEBUG_ 
#include "../program_contest_library/memo/dump.hpp"
#else
#define dump(...)
#endif
const ll INF = 1LL<<60;

int main(void) {
    ll n, k;
    cin >> n >> k;
    vector<vector<PII>> g(n);
    REP(i, n-1) {
        ll u, v, c;
        cin >> u >> v >> c;
        u--, v--, c--;
        g[u].push_back({v, c});
        g[v].push_back({u, c});
    }

    vector<ll> sz(n), dead(n);
    auto find_centroid = [&](ll root) {
        auto get_size = [&](auto &&self, ll v, ll p) -> void {
            sz[v] = 1;
            for(auto to: g[v]) if(to.first != p && !dead[to.first]) {
                self(self, to.first, v);
                sz[v] += sz[to.first];
            }
        };
        get_size(get_size, root, -1);
        auto dfs = [&](auto &&self, ll v, ll p) -> ll {
            for(auto to: g[v]) if(to.first != p && !dead[to.first]) {
                if(sz[to.first] > sz[root]/2) return self(self, to.first, v);
            }
            return v;
        };
        return dfs(dfs, root, -1);
    };

    auto get_edges = [&](ll v) -> vector<vector<PII>> {
        vector<vector<PII>> ret;

        auto dfs = [&](auto &&self, ll v, ll p, PII col) -> void {
            for(auto to: g[v]) if(to.first != p && !dead[to.first]) {                
                PII now = col;
                if(to.second == now.first || to.second == now.second) {
                    ret.back().push_back(now);
                    self(self, to.first, v, now);
                } else if(now.first == -1) {
                    now.first = to.second;
                    if(now.first > now.second) swap(now.first, now.second);
                    ret.back().push_back(now);
                    self(self, to.first, v, now);
                }
            }
        };

        for(auto to: g[v]) if(!dead[to.first]) {
            ret.emplace_back();
            ret.back().emplace_back(-1, to.second);
            dfs(dfs, to.first, -1, PII(-1, to.second));
            sort(ALL(ret.back()));
        }
        return ret;
    };

    ll ans = 0;
    auto centroid_decomposition = [&](auto &&self, ll root) -> void {
        ll c = find_centroid(root);
        dead[c] = true;
        for(auto to: g[c]) if(!dead[to.first]) self(self, to.first);

        // (色1,色2) を含むパスとマージさせられるパスは (色1,色2) or (-1,色1) or (-1,色2) のどれか
        // (色1,色2) のペアが何個あるのか数えるための用意
        auto colors = get_edges(c);
        vector<PII> flat;
        for(auto v: colors) for(auto p: v) flat.push_back(p);
        sort(ALL(flat));
        dump(c);
        dump(colors);
        dump(flat);

        for(auto v: colors) for(auto p: v) {
            dump(p, ans);
            // 2色存在するペアと合わせられるものを数える
            if(p.first != -1 && p.second != -1) {
                // 頂点cまでのパスを数える
                ans += 2;
                // 同じ色のペア
                {
                    auto [lb, ub] = equal_range(ALL(flat), p);
                    ans += ub - lb;
                    auto [lb2, ub2] = equal_range(ALL(v), p);
                    ans -= ub2 - lb2; // 同じ方向のパスから2回選ぶのはだめ
                }
                // (-1,色1)
                {
                    auto [lb, ub] = equal_range(ALL(flat), PII(-1, p.first));
                    ans += 2 * (ub - lb);
                    auto [lb2, ub2] = equal_range(ALL(v), PII(-1, p.first));
                    ans -= 2 * (ub2 - lb2);
                }
                // (-1,色2)
                {
                    auto [lb, ub] = equal_range(ALL(flat), PII(-1, p.second));
                    ans += 2 * (ub - lb);
                    auto [lb2, ub2] = equal_range(ALL(v), PII(-1, p.second));
                    ans -= 2 * (ub2 - lb2);
                }
            } 
            // (-1, 色1) と合わせられるやつを探す
            else {
                // (-1, 色i) (色i!=色1) 
                {
                    auto ub = upper_bound(ALL(flat), PII(-1, INF));
                    auto lb = lower_bound(ALL(flat), PII(-1, 0));
                    ans += ub - lb;
                    auto ub2 = upper_bound(ALL(v), PII(-1, INF));
                    auto lb2 = lower_bound(ALL(v), PII(-1, 0));
                    ans -= ub2 - lb2;
                    // auto [lb3, ub3] = equal_range(ALL(flat), PII(-1, p.second));
                    // ans -= ub3 - lb3;
                    // auto [lb4, ub4] = equal_range(ALL(v), PII(-1, p.second));
                    // ans += ub4 - lb4;
                }
                // (色1, 色*)
                {
                    auto ub = upper_bound(ALL(flat), PII(p.second, INF));
                    auto lb = lower_bound(ALL(flat), PII(p.second, 0));
                    ans += ub - lb;
                    auto ub2 = upper_bound(ALL(v), PII(p.second, INF));
                    auto lb2 = lower_bound(ALL(v), PII(p.second, 0));
                    ans -= ub2 - lb2;
                }
            }
            dump(ans);
        }

        dead[c] = false;
    };

    centroid_decomposition(centroid_decomposition, 0);
    cout << ans/2 << endl;

    return 0;
}

yukicoder

結果

テストケース

ソースコード