結果

問題 No.922 東北きりきざむたん
ユーザー polylogKpolylogK
提出日時 2019-11-01 18:21:46
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 3,491 bytes
コンパイル時間 2,419 ms
コンパイル使用メモリ 188,688 KB
実行使用メモリ 36,932 KB
最終ジャッジ日時 2024-09-14 22:48:07
合計ジャッジ時間 15,422 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 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 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 3 ms
5,376 KB
testcase_09 AC 459 ms
20,780 KB
testcase_10 AC 38 ms
8,704 KB
testcase_11 AC 181 ms
16,968 KB
testcase_12 AC 1,401 ms
21,868 KB
testcase_13 AC 75 ms
8,192 KB
testcase_14 AC 1,456 ms
30,476 KB
testcase_15 AC 1,805 ms
23,180 KB
testcase_16 AC 332 ms
33,736 KB
testcase_17 AC 367 ms
33,860 KB
testcase_18 AC 363 ms
33,728 KB
testcase_19 AC 346 ms
33,736 KB
testcase_20 AC 334 ms
33,732 KB
testcase_21 AC 206 ms
33,688 KB
testcase_22 AC 208 ms
33,604 KB
testcase_23 AC 227 ms
34,368 KB
testcase_24 AC 239 ms
34,376 KB
testcase_25 AC 190 ms
35,396 KB
testcase_26 AC 191 ms
35,396 KB
testcase_27 AC 192 ms
35,396 KB
testcase_28 TLE -
testcase_29 AC 172 ms
36,932 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std::literals::string_literals;
using i64 = long long;
using std::cout;
using std::endl;
using std::cin;

template<typename T>
std::vector<T> make_v(size_t a){return std::vector<T>(a);}

template<typename T,typename... Ts>
auto make_v(size_t a,Ts... ts){
  return std::vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...));
}

class lowest_common_ancestor {
	std::vector<std::vector<size_t>> g, parent;
	std::vector<size_t> depth;
	
	const size_t logN = 20;

	private:
	void dfs(size_t v, int par = -1, size_t d = 0) {
		parent[0][v] = par;
		depth[v] = d;
		
		for(auto e:g[v]) {
			if(e == par) continue;
			
			dfs(e, v, d + 1);
		}
	}
	void build() {
		for(int i = 0; i < g.size(); i++) if(depth[i] == -1) dfs(i);
		
		for(int k = 0; k < logN - 1; k++) {
			for(int i = 0; i < g.size(); i++) {
				if(parent[k][i] == -1) parent[k + 1][i] = -1;
				else parent[k + 1][i] = parent[k][parent[k][i]];
			}
		}
	}

	public:
	lowest_common_ancestor() {}
	lowest_common_ancestor(int n) : g(n), depth(n, -1) {
		parent.assign(logN, std::vector<size_t>(n, -1));
	}
	void add_edge(const size_t & u, const size_t & v) {
		g[u].push_back(v);
		g[v].push_back(u);
	}
	const size_t query(size_t u, size_t v) {
		static bool built = false;
		if(!built) {
			build();
			built = true;
		}

		if(depth[u] > depth[v]) std::swap(u, v);
		for(int k = 0; k < logN; k++) {
			if(((depth[v] - depth[u]) >> k) & 1) v = parent[k][v];
		}
		if(u == v) return u;

		for(int k = logN - 1; k >= 0; k--) {
			if(parent[k][u] == parent[k][v]) continue;

			u = parent[k][u];
			v = parent[k][v];			
		}		
		return parent[0][u];
	}
};

int main() {
	int n, m, q; scanf("%d%d%d", &n, &m, &q);
	lowest_common_ancestor lca(n);
	std::vector<std::vector<int>> g(n);
	for(int i = 0; i < m; i++) {
		int a, b; scanf("%d%d", &a, &b);

		g[a - 1].push_back(b - 1);
		g[b - 1].push_back(a - 1);
		lca.add_edge(a - 1, b - 1);
	}
	
	int cnt = 0;
	std::vector<int> comp(n, -1), dist(n);
	for(int i = 0; i < n; i++) {
		if(comp[i] != -1) continue;
		
		auto kiri = [&](auto && kiri, int v, int par) -> void {
			comp[v] = cnt;
			
			for(auto e: g[v]) {
				if(e == par) continue;

				dist[e] = dist[v] + 1;
				kiri(kiri, e, v);
			}
		};
		kiri(kiri, i, -1);
		cnt++;
	}
	auto calc = [dist, &lca](int a, int b) -> i64 {
		return dist[a] + dist[b] - 2 * dist[lca.query(a, b)];
	};

	i64 ans = 0;
	std::vector<i64> p(n);
	while(q--) {
		int a, b; scanf("%d%d", &a, &b); a--; b--;

		if(comp[a] != comp[b]) {
			p[a]++;
			p[b]++;
		} else {
			ans += calc(a, b);
		}
	}
	
	std::vector<bool> used(cnt, false);
	for(int i = 0; i < n; i++) {
		if(used[comp[i]]) continue;
		used[comp[i]] = true;

		i64 now = 0;
		auto dfs_A = [&](auto && dfs_A, int v, int par, i64 d) -> void {
			now += d * p[v];

			for(auto e: g[v]) {
				if(e == par) continue;

				dfs_A(dfs_A, e, v, d + 1);
			}
		};
		dfs_A(dfs_A, i, -1, 0);

		std::vector<i64> sz(n, 0);
		auto dfs_B = [&](auto && dfs_B, int v, int par) -> void {
			sz[v] += p[v];
			
			for(auto e: g[v]) {
				if(e == par) continue;

				dfs_B(dfs_B, e, v);
				sz[v] += sz[e];
			}
		};
		dfs_B(dfs_B, i, -1);

		i64 tmp = now;
		auto dfs_C = [&](auto && dfs_C, int v, int par) -> void {
			tmp = std::min(tmp, now);
			
			for(auto e: g[v]) {
				if(e == par) continue;

				now += sz[i] - 2 * sz[e];
				dfs_C(dfs_C, e, v);
				now -= sz[i] - 2 * sz[e];
			}
		};
		dfs_C(dfs_C, i, -1);

		ans += tmp;
	}

	printf("%lld\n", ans);
	return 0;
}
0