結果

問題 No.650 行列木クエリ
ユーザー eve__fuyukieve__fuyuki
提出日時 2024-10-25 14:36:18
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 89 ms / 2,000 ms
コード長 2,680 bytes
コンパイル時間 3,262 ms
コンパイル使用メモリ 232,116 KB
実行使用メモリ 22,916 KB
最終ジャッジ日時 2024-10-25 14:36:23
合計ジャッジ時間 4,722 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 28 ms
6,820 KB
testcase_02 AC 89 ms
17,336 KB
testcase_03 AC 2 ms
6,816 KB
testcase_04 AC 29 ms
6,820 KB
testcase_05 AC 84 ms
17,244 KB
testcase_06 AC 2 ms
6,816 KB
testcase_07 AC 2 ms
6,820 KB
testcase_08 AC 28 ms
7,200 KB
testcase_09 AC 72 ms
22,916 KB
testcase_10 AC 2 ms
6,824 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

void fast_io() {
	ios_base::sync_with_stdio(false);
	cin.tie(nullptr);
}
#include <atcoder/modint>
#include <atcoder/segtree>
using namespace atcoder;
using mint = modint1000000007;
struct S {
	mint x00, x01, x10, x11;
};
S op(S a, S b) {
	return S{a.x00 * b.x00 + a.x01 * b.x10, a.x00 * b.x01 + a.x01 * b.x11,
			 a.x10 * b.x00 + a.x11 * b.x10, a.x10 * b.x01 + a.x11 * b.x11};
}
S e() { return S{1, 0, 0, 1}; }
int main() {
	fast_io();
	int n;
	cin >> n;
	vector<vector<int>> g(n);
	vector<pair<int, int>> idx(n - 1);
	for (int i = 0; i < n - 1; i++) {
		int a, b;
		cin >> a >> b;
		g[a].push_back(b);
		g[b].push_back(a);
		idx[i] = {a, b};
	}
	vector<int> sz(n), par(n);

	{
		auto dfs = [&](auto self, int u, int p) -> void {
			sz[u] = 1;
			for (int v : g[u]) {
				if (v == p) continue;
				par[v] = u;
				self(self, v, u);
				sz[u] += sz[v];
			}
		};
		dfs(dfs, 0, -1);
		vector<vector<int>> g_new(n);
		for (int i = 0; i < n; i++) {
			for (int j : g[i]) {
				if (j != par[i]) {
					g_new[i].push_back(j);
				}
			}
			sort(g_new[i].begin(), g_new[i].end(),
				 [&](int a, int b) { return sz[a] > sz[b]; });
		}
		g = g_new;
	}
	vector<int> dfs_ord(n), root(n);
	{
		iota(root.begin(), root.end(), 0);
		vector<int> ord;
		auto dfs = [&](auto self, int u) -> void {
			dfs_ord[u] = ord.size();
			ord.push_back(u);
			for (int v : g[u]) {
				if (v == g[u][0]) {
					root[v] = root[u];
				}
				self(self, v);
			}
		};
		dfs(dfs, 0);
	}
	for (int i = 0; i < n - 1; i++) {
		if (idx[i].first != par[idx[i].second]) {
			swap(idx[i].first, idx[i].second);
		}
	}
	vector<int> bfs_ord(n);
	{
		vector<int> ord;
		deque<int> dq{0};
		int cnt = 0;
		while (!dq.empty()) {
			int u = dq.front();
			dq.pop_front();
			bfs_ord[u] = cnt++;
			ord.push_back(u);
			for (int v : g[u]) {
				dq.push_back(v);
			}
		}
	}
	vector<S> seg_data(n, e());
	segtree<S, op, e> seg(seg_data);
	int q;
	cin >> q;
	for (; q--;) {
		char com;
		cin >> com;
		if (com == 'x') {
			int i, x00, x01, x10, x11;
			cin >> i >> x00 >> x01 >> x10 >> x11;
			int j = dfs_ord[idx[i].second];
			seg.set(j, S{x00, x01, x10, x11});
		}
		if (com == 'g') {
			int i, j;
			cin >> i >> j;
			S ans = e();
			while (i != j) {
				if (root[j] == j) {
					ans = op(seg.get(dfs_ord[j]), ans);
					j = par[j];
				} else if (root[j] != root[i]) {
					ans = op(seg.prod(dfs_ord[root[j]], dfs_ord[j] + 1), ans);
					j = par[root[j]];
				} else {
					ans = op(seg.prod(dfs_ord[i] + 1, dfs_ord[j] + 1), ans);
					break;
				}
			}
			cout << ans.x00.val() << " " << ans.x01.val() << " "
				 << ans.x10.val() << " " << ans.x11.val() << "\n";
		}
	}
}
0