結果

問題 No.1769 Don't Stop the Game
ユーザー polylogKpolylogK
提出日時 2021-11-03 08:31:16
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,370 ms / 3,000 ms
コード長 2,633 bytes
コンパイル時間 936 ms
コンパイル使用メモリ 80,960 KB
実行使用メモリ 112,344 KB
最終ジャッジ日時 2024-06-29 17:57:17
合計ジャッジ時間 18,573 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 3 ms
6,944 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 3 ms
6,940 KB
testcase_07 AC 3 ms
6,940 KB
testcase_08 AC 326 ms
32,060 KB
testcase_09 AC 386 ms
26,816 KB
testcase_10 AC 753 ms
49,556 KB
testcase_11 AC 231 ms
24,124 KB
testcase_12 AC 255 ms
38,056 KB
testcase_13 AC 258 ms
38,652 KB
testcase_14 AC 263 ms
38,692 KB
testcase_15 AC 306 ms
38,840 KB
testcase_16 AC 432 ms
39,028 KB
testcase_17 AC 770 ms
45,544 KB
testcase_18 AC 1,239 ms
68,380 KB
testcase_19 AC 1,335 ms
79,544 KB
testcase_20 AC 1,370 ms
81,320 KB
testcase_21 AC 1,336 ms
81,536 KB
testcase_22 AC 1,350 ms
81,548 KB
testcase_23 AC 218 ms
38,620 KB
testcase_24 AC 307 ms
38,632 KB
testcase_25 AC 103 ms
35,636 KB
testcase_26 AC 1,049 ms
77,236 KB
testcase_27 AC 209 ms
79,768 KB
testcase_28 AC 1,302 ms
112,344 KB
testcase_29 AC 1,299 ms
93,556 KB
testcase_30 AC 1,208 ms
93,640 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/*
	# Algorithm



	## Time Comlexity

	O(Nlog(N)^2)
*/

#include <stdio.h>
#include <algorithm>
#include <vector>
#include <map>

int main() {
	int n; scanf("%d", &n);
	std::vector<std::vector<int>> g(n);
	std::vector<int> a(n - 1), b(n - 1), c(n - 1);
	for(int i = 0; i < n - 1; i++) {
		scanf("%d%d%d", &a[i], &b[i], &c[i]); a[i]--; b[i]--;

		g[a[i]].push_back(i);
		g[b[i]].push_back(i);
	}

	std::vector<int> x(n), y(n), size(n), in(n), out(n), tour;
	std::vector<std::vector<int>> child(n), child_order(n); // LA
	std::map<int, std::vector<std::vector<int>>> map_xor; {
		std::map<int, int> map;
		auto dfs = [&](auto&& dfs, int v, int par, int xor_val) -> void {
			in[v] = tour.size();
			tour.push_back(v);

			x[v] = xor_val;
			y[v] = map[xor_val]++;
			auto& vec = map_xor[xor_val];
			if(vec.size() <= y[v]) vec.resize(y[v] + 2);
			vec[y[v]].push_back(in[v]);

			size[v] = 1;
			for(int id: g[v]) {
				int to = v ^ a[id] ^ b[id];
				if(to == par) continue;

				dfs(dfs, to, v, xor_val ^ c[id]);
				child[v].push_back(to);
				child_order[v].push_back(in[to]);
				size[v] += size[to];
			}

			tour.push_back(v);
			out[v] = tour.size();

			map[xor_val]--;
		}; dfs(dfs, 0, -1, 0);
	}

	auto get_child = [&](int v, int id) -> int { // 頂点 v の子であって id が部分木に含まれるようなものを返す
		int k = std::prev(upper_bound(begin(child_order[v]), end(child_order[v]), id))
			  	- begin(child_order[v]);
		return child[v][k];
	};

	using i64 = long long;
	i64 ans = 0;
	{ // 解説のパターン 1, 2 を数える
		std::map<int, int> map;
		auto dfs = [&](auto&& dfs, int v, int par, int depth) -> void {
			if(y[v]) ans += map[x[v]] + size[v];

			int tmp = map[x[v]];
			for(int id: g[v]) {
				int to = a[id] ^ b[id] ^ v;
				if(to == par) continue;

				map[x[v]] = n - size[to];
				dfs(dfs, to, v, depth + 1);
			}
			map[x[v]] = tmp;
		}; dfs(dfs, 0, -1, 0);
	}
	{ // 解説のパターン 3 を数える
		std::map<int, int> map;
		for(int i = 0; i < n; i++) map[x[i]] = -1;
		auto dfs = [&](auto&& dfs, int v, int par) -> void {
			if(v and map[x[v]] != par) {
				int root = (y[v] ? get_child(map[x[v]], in[v]) : 0);
				auto& vec = map_xor[x[v]][y[v]];

				i64 cnt = std::distance(lower_bound(begin(vec), end(vec), in[root]),
						 				lower_bound(begin(vec), end(vec), out[root])) - 1;
				ans += cnt * size[v];
			}

			int tmp = map[x[v]]; map[x[v]] = v;
			for(int id: g[v]) {
				int to = a[id] ^ b[id] ^ v;
				if(to == par) continue;

				dfs(dfs, to, v);
			}
			map[x[v]] = tmp;
		}; dfs(dfs, 0, -1);
	}

	ans = (i64)n * (n - 1) - ans;
	printf("%lld\n", ans);
}
0