結果

問題 No.399 動的な領主
ユーザー ctyl_0ctyl_0
提出日時 2016-07-16 17:37:11
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 195 ms / 2,000 ms
コード長 3,520 bytes
コンパイル時間 1,267 ms
コンパイル使用メモリ 105,528 KB
実行使用メモリ 31,296 KB
最終ジャッジ日時 2024-04-25 08:58:17
合計ジャッジ時間 3,872 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 12 ms
5,376 KB
testcase_06 AC 195 ms
24,248 KB
testcase_07 AC 180 ms
24,252 KB
testcase_08 AC 168 ms
24,380 KB
testcase_09 AC 170 ms
24,508 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 10 ms
5,376 KB
testcase_12 AC 119 ms
25,400 KB
testcase_13 AC 117 ms
25,276 KB
testcase_14 AC 77 ms
31,292 KB
testcase_15 AC 86 ms
31,296 KB
testcase_16 AC 92 ms
28,216 KB
testcase_17 AC 164 ms
24,380 KB
testcase_18 AC 168 ms
24,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <iomanip>
#include <vector>
#include <algorithm>
#include <numeric>
#include <functional>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <sstream>
#include <string>
#define repd(i,a,b) for (int i=(int)(a);i<(int)(b);i++)
#define rep(i,n) repd(i,0,n)
#define all(x) (x).begin(),(x).end()
#define mod 1000000007
#define inf 2000000007
#define mp make_pair
#define pb push_back
typedef long long ll;
using namespace std;
template <typename T>
inline void output(T a, int p) {
    if(p) cout << fixed << setprecision(p)  << a << "\n";
    else cout << a << "\n";
}
// end of template

// least common ancestor
class LCA{
public:
    int V;
    int logV, id = 0;
    vector<vector<pair<int, int>>> G; // graph: vertex, distance
    vector<int> depth, dist;
    vector<vector<int>> parent;
    vector<int> L, R, I;
    
    LCA(int n){
        V = n;
        G.resize(V);
        depth.resize(V);
        dist.resize(V);
        L.resize(V), R.resize(V), I.resize(V);
        logV = 0;
        while (V >= (1 << logV)) logV++;
        parent.resize(logV, vector<int>(V));
        
    }
    
    // determine depth and distance from root
    void init(int v = 0, int par = -1, int d = 0, int l = 0) { // root: v = 0, init(0, -1, 0, 0)
        depth[v] = d;
        parent[0][v] = par;
        dist[v] = l;
        I[id] = v;
        L[v] = id++;
        rep(i, G[v].size()){
            int w = G[v][i].first;
            int lc = G[v][i].second;
            if (w == par) continue;
            init(w, v, d + 1, lc + l);
        }
        R[v] = id;
    }
    
    void build() {
        for (int k = 0; k + 1 < logV; k++) {
            for (int v = 0; v < V; v++) {
                if (parent[k][v] < 0) parent[k+1][v] = -1;
                else parent[k+1][v] = parent[k][parent[k][v]];
            }
        }
    }
    
    void add(int v1, int v2, int dist = 0){
        G[v1].pb(mp(v2, dist));
        G[v2].pb(mp(v1, dist));
    }
    
    int query(int u, int v) {
        if (depth[u] > depth[v]) swap(u, v);
        for (int k = 0; k < logV; k++) {
            if ((depth[v] - depth[u]) >> k & 1)
                v = parent[k][v];
        }
        if (u == v) return u;
        
        for (int k = logV - 1; k >= 0; k--) {
            if (parent[k][u] != parent[k][v]) {
                u = parent[k][u];
                v = parent[k][v];
            }
        }
        return parent[0][u];
    }
    
};


void dfs(vector<vector<pair<int, int>>> &G, vector<ll> &dist, int pre = -1, int cur = 0){
    for(auto v: G[cur]){
        if(v.first != pre){
            dfs(G, dist, cur, v.first);
            dist[cur] += dist[v.first];
        }
    }
}



int main() {
    cin.tie(0);
    ios::sync_with_stdio(0);
    // source code
    int N;
    cin >> N;
    LCA lca(N);
    rep(i, N - 1){
        int A, B;
        cin >> A >> B;
        A--, B--;
        lca.add(A, B);
    }
    lca.init();
    lca.build();
    auto G = lca.G;
    vector<ll> dist(N, 0);
    
    int M;
    cin >> M;
    ll ret = 0;
    rep(i, M){
        int A, B;
        cin >> A >> B;
        A--, B--;
        dist[A]++;
        dist[B]++;
        int l = lca.query(A, B);
        dist[l]--;
        if(lca.parent[0][l] != -1) dist[lca.parent[0][l]]--;
//        cout << A << ", " << B << ", " << l << ", " << lca.parent[0][l] << endl;
    }
    
    dfs(G, dist);
    
    rep(i, N) ret += dist[i] * (dist[i] + 1) / 2;
    
    output(ret, 0);
    
    return 0;
}
0