結果

問題 No.399 動的な領主
コンテスト
ユーザー ryofjt
提出日時 2016-09-12 11:07:43
言語 C++11(old_compat)
(gcc 12.4.0 + boost 1.89.0)
コンパイル:
g++-12 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -include bits/stdc++.h -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 247 ms / 2,000 ms
コード長 1,326 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,280 ms
コンパイル使用メモリ 170,568 KB
実行使用メモリ 23,140 KB
最終ジャッジ日時 2026-03-08 16:07:24
合計ジャッジ時間 4,404 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>

#define rep(i, n) for(int i = 0; i < (n); ++i)

using namespace std;

typedef long long ll;
typedef vector<int> V;

int n;
V t[100000];
int d[100000];
int k[100000][17];
ll v[100000];

void dfs(int x){
	for(int i: t[x]){
		if(d[i] == -1){
			d[i] = d[x] + 1;
			k[i][0] = x;
			dfs(i);			
		}
	}
}
void init(){
	fill_n(d, n, -1);
	d[0] = 0;
	k[0][0] = -1;
	dfs(0);
	for(int j = 1; j < 17; ++j){
		rep(i, n){
			k[i][j] = k[i][j - 1] >= 0 ? k[k[i][j - 1]][j - 1] : -1;
		}
	}
}

int up(int v, int x){
	rep(i, 17){
		if(x >> i & 1){
			v = k[v][i];
		}
	}
	return v;
}
int lca(int u, int v){
	if(d[u] > d[v]){
		swap(u, v);
	}
	v = up(v, d[v] - d[u]);
	if(u == v){
		return u;
	}
	for(int i = 16; i >= 0; --i){
		if(k[u][i] != k[v][i]){
			u = k[u][i];
			v = k[v][i];
		}
	}
	return k[u][0];
}

int cs(int x = 0){
	for(int i: t[x]){
		if(d[x] < d[i]){
			v[x] += cs(i);
		}
	}
	return v[x];
}

int main(){
	cin >> n;
	rep(i, n - 1){
		int u, v;
		cin >> u >> v;
		t[--u].push_back(--v);
		t[v].push_back(u);
	}

	init();

	int q;
	cin >> q;
	rep(i, q){
		int a, b;
		cin >> a >> b;
		int c = lca(--a, --b);
		--v[c];
		if(k[c][0] != -1){
			--v[k[c][0]];
		}
		++v[a];
		++v[b];
	}

	cs();

	ll ans = 0;
	rep(i, n){
		ans += v[i] * (v[i] + 1) / 2;
	}
	cout << ans << endl;
	return 0;
}
0