結果

問題 No.922 東北きりきざむたん
ユーザー arknavearknave
提出日時 2019-11-08 22:17:26
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 112 ms / 2,000 ms
コード長 3,022 bytes
コンパイル時間 1,921 ms
コンパイル使用メモリ 172,944 KB
実行使用メモリ 25,856 KB
最終ジャッジ日時 2024-09-15 01:34:47
合計ジャッジ時間 4,721 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 8 ms
12,928 KB
testcase_01 AC 9 ms
12,928 KB
testcase_02 AC 8 ms
12,800 KB
testcase_03 AC 9 ms
12,928 KB
testcase_04 AC 9 ms
12,928 KB
testcase_05 AC 9 ms
12,800 KB
testcase_06 AC 9 ms
12,672 KB
testcase_07 AC 9 ms
12,800 KB
testcase_08 AC 9 ms
12,800 KB
testcase_09 AC 39 ms
16,384 KB
testcase_10 AC 29 ms
13,824 KB
testcase_11 AC 38 ms
15,700 KB
testcase_12 AC 23 ms
16,076 KB
testcase_13 AC 16 ms
13,824 KB
testcase_14 AC 56 ms
18,304 KB
testcase_15 AC 19 ms
16,100 KB
testcase_16 AC 96 ms
19,200 KB
testcase_17 AC 95 ms
19,072 KB
testcase_18 AC 91 ms
19,072 KB
testcase_19 AC 93 ms
19,072 KB
testcase_20 AC 94 ms
19,072 KB
testcase_21 AC 105 ms
19,328 KB
testcase_22 AC 103 ms
19,328 KB
testcase_23 AC 111 ms
18,944 KB
testcase_24 AC 112 ms
18,944 KB
testcase_25 AC 94 ms
19,712 KB
testcase_26 AC 95 ms
19,584 KB
testcase_27 AC 94 ms
19,584 KB
testcase_28 AC 39 ms
16,660 KB
testcase_29 AC 100 ms
25,856 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

#define all(x) begin(x), end(x)

using ll = long long;
using ld = long double;
using pii = pair<int, int>;
using vi = vector<int>;

constexpr int LOGN  = 18;
constexpr int MAXN = 1e5 + 5;
int n, m, q;
vi tree[MAXN];
int comp[MAXN];
int depth[MAXN];
int st[LOGN][MAXN];

ll tokens[MAXN];
ll sub_cost[MAXN];
ll dp[MAXN];

int lca(int u, int v) {
    assert(comp[u] == comp[v]);
    if (depth[u] > depth[v])
        swap(u, v);
    for (int j = LOGN - 1; j >= 0; --j) {
        if ((depth[v] - depth[u]) & (1 << j)) {
            v = st[j][v];
        }
    }

    if (u == v)
        return u;

    for (int j = LOGN - 1; j >= 0; --j) {
        if (st[j][u] != st[j][v]) {
            u = st[j][u];
            v = st[j][v];
        }
    }

    return st[0][u];
}

int dist(int u, int v) {
    return depth[u] + depth[v] - 2 * depth[lca(u, v)];
}

void dfs1(int u, int p, int r) {
    comp[u] = r;
    for (int v : tree[u]) {
        if (v == p) continue;
        depth[v] = depth[u] + 1;
        st[0][v] = u;
        dfs1(v, u, r);
    }
}

void dfs2(int u, int p) {
    sub_cost[u] = 0;
    for (int v : tree[u]) {
        if (v == p) continue;

        dfs2(v, u);
        tokens[u] += tokens[v];
        sub_cost[u] += sub_cost[v] + tokens[v];
    }
}

void dfs3(int u, int p, ll cost) {
    dp[u] = cost;
    for (int v : tree[u]) {
        if (v == p) continue;

        ll new_cost = cost + tokens[comp[u]] - 2 * tokens[v];
        dfs3(v, u, new_cost);
    }
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);

    // Given a forest and a bunch of travel plans
    // If travel within a CC, fixed contribution
    // Otherwise, do some kind of DFS to get APSP? Tree re-rooting trick?
    // Take min on each component

    cin >> n >> m >> q;
    for (int i = 0; i < m; ++i) {
        int u, v;
        cin >> u >> v;
        --u; --v;
        tree[u].push_back(v);
        tree[v].push_back(u);
    }

    memset(st, -1, sizeof(st));
    vi roots;
    for (int i = 0; i < n; ++i) {
        if (depth[i] == 0) {
            roots.push_back(i);
            dfs1(i, -1, i);
        }
    }

    for (int j = 0; j + 1 < LOGN; ++j) {
        for (int u = 0; u < n; ++u) {
            if (st[j][u] != -1) {
                st[j + 1][u] = st[j][st[j][u]];
            }
        }
    }

    ll ans = 0;
    while (q-- > 0) {
        int a, b;
        cin >> a >> b;
        --a; --b;
        if (comp[a] == comp[b]) {
            ans += dist(a, b);
        } else {
            ++tokens[a];
            ++tokens[b];
        }
    }

    for (int root : roots) {
        dfs2(root, -1);
        dfs3(root, -1, sub_cost[root]);
    }

    constexpr ll INF = 1e18;
    vector<ll> res(roots.size(), INF);
    for (int i = 0; i < n; ++i) {
        int idx = lower_bound(all(roots), comp[i]) - begin(roots);
        res[idx] = min(res[idx], dp[i]);
    }

    for (ll v : res) {
        ans += v;
    }

    cout << ans << '\n';

    return 0;
}
0