結果

問題 No.399 動的な領主
コンテスト
ユーザー Pachicobue
提出日時 2017-11-09 14:08:53
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 4,002 bytes
コンパイル時間 2,239 ms
コンパイル使用メモリ 185,820 KB
実行使用メモリ 64,756 KB
最終ジャッジ日時 2024-11-24 07:25:51
合計ジャッジ時間 9,796 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 RE * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define show(x) cerr << #x << " = " << x << endl

using namespace std;
using ll = long long;
using pii = pair<int, int>;
using vi = vector<int>;

template <typename T>
ostream& operator<<(ostream& os, const vector<T>& v)
{
    os << "sz=" << v.size() << "\n[";
    for (const auto& p : v) {
        os << p << ",";
    }
    os << "]\n";
    return os;
}

template <typename S, typename T>
ostream& operator<<(ostream& os, const pair<S, T>& p)
{
    os << "(" << p.first << "," << p.second
       << ")";
    return os;
}


constexpr ll MOD = 1e9 + 7;

template <typename T>
constexpr T INF = numeric_limits<T>::max() / 100;

struct Graph {
    Graph(const int n)
    {
        edge.resize(n);
    }
    void addEdge(const int from, const int to)
    {
        edge[from].push_back(to);
    }
    vector<vector<int>> edge;
};

template <typename Base>
class SparseTable
{
public:
    using T = typename Base::T;
    using SemiLattice = Base;
    SparseTable(const vector<T>& val) : size(val.size()), lg2(size + 1, 0)
    {
        for (int i = 2; i <= size; i++) {
            lg2[i] = lg2[i / 2] + 1;
        }
        table.resize(size, vector<T>(lg2[size] + 1));
        for (int i = 0; i < size; i++) {
            table[i][0] = val[i];
        }
        for (int j = 0; j < lg2[size]; j++) {
            const int w = 1 << j;
            for (int i = 0; i <= size - (w << 1); i++) {
                T tl = table[i][j], tr = table[i + w][j];
                table[i][j + 1] = op(tl, tr);
            }
        }
    }

    T accumulate(const int l, const int r) const
    {
        assert(0 <= l and l < r and r <= size);
        const int j = lg2[r - l];
        return op(table[l][j], table[r - (1 << j)][j]);
    }


private:
    const int size;
    vector<int> lg2;
    vector<vector<T>> table;
    const SemiLattice op{};
};

struct Min {
    using T = pii;
    T operator()(const T& a, const T& b) const
    {
        return min(a, b);
    }
};

void make_dfs(Graph& g, const Graph& g_, vector<bool>& used, vector<int>& prev, const int s)
{
    used[s] = true;
    for (const int to : g_.edge[s]) {
        if (not used[to]) {
            g.addEdge(s, to);
            prev[to] = s;
            make_dfs(g, g_, used, prev, to);
        }
    }
}

void make_dfs(const Graph& g, vector<pii>& nodes, vector<int>& left, const int s, const int depth)
{
    nodes.push_back(make_pair(depth, s));
    left[s] = nodes.size() - 1;
    for (const int to : g.edge[s]) {
        make_dfs(g, nodes, left, to, depth + 1);
        nodes.push_back(make_pair(depth, s));
    }
}

int acc_dfs(const Graph& g, vector<int>& imos, const int s)
{
    int sum = 0;
    for (const int to : g.edge[s]) {
        sum += acc_dfs(g, imos, to);
    }
    imos[s] += sum;
    return imos[s];
}

int LCA(const SparseTable<Min>& stable, const vector<int>& left, const int u, const int v)
{
    const int l = left[u];
    const int r = left[v] + 1;
    return stable.accumulate(l, r).second;
}

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);
    int N;
    cin >> N;
    Graph g_(N);
    for (int i = 0; i < N - 1; i++) {
        int u, v;
        cin >> u >> v;
        u--, v--;
        g_.addEdge(u, v);
        g_.addEdge(v, u);
    }
    Graph g(N);
    vector<bool> used(N, false);
    vector<int> prev(N, -1);
    make_dfs(g, g_, used, prev, 0);
    vector<pii> nodes;
    vector<int> left(N, 0);
    make_dfs(g, nodes, left, 0, 0);
    SparseTable<Min> stable(nodes);

    vector<int> imos(N, 0);
    int Q;
    cin >> Q;
    for (int i = 0; i < Q; i++) {
        int a, b;
        cin >> a >> b;
        a--, b--;
        imos[a]++;
        imos[b]++;
        const int lca = LCA(stable, left, a, b);
        imos[lca]--;
        if (prev[lca] != -1) {
            imos[prev[lca]]--;
        }
    }
    acc_dfs(g, imos, 0);
    ll sum = 0;
    for (int i = 0; i < N; i++) {
        sum += (ll)imos[i] * ((ll)imos[i] + 1) / 2;
    }
    cout << sum << endl;

    return 0;
}
0