結果

問題 No.2949 Product on Tree
ユーザー NakLon131
提出日時 2024-10-26 10:51:24
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 161 ms / 2,000 ms
コード長 1,733 bytes
コンパイル時間 12,475 ms
コンパイル使用メモリ 401,052 KB
実行使用メモリ 35,068 KB
最終ジャッジ日時 2024-10-26 10:51:52
合計ジャッジ時間 22,006 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

fn main() {
    input!{
		n: usize,
		a: [usize; n],
		uv: [(Usize1, Usize1); n-1],
	}

	let mut g = vec![Vec::new(); n];
	for &(u, v) in &uv {
		g[u].push(v);
		g[v].push(u);
	}

	let mut ans = 0;
	dfs(0, MOD93, &g, &a, &mut ans);
	println!("{}", ans);

	fn dfs(pos: usize, ppos: usize, g: &Vec<Vec<usize>>, a: &Vec<usize>, ans: &mut usize) -> usize {
		let mut sm1 = 0;
		let mut sm2 = 0;
		for &npos in &g[pos] {
			if npos == ppos { continue; }
			let esm = dfs(npos, pos, g, a, ans) % MOD93;
			sm1 += esm; 		sm1 %= MOD93;
			sm2 += esm * esm; 	sm2 %= MOD93;
		}

		// ans+=a[p]*(sm1+(sm1*sm1-sm2)/2);
		let multi = (sm1*sm1+MOD93-sm2) % MOD93;
		let divs = div(multi, 2, MOD93);
		*ans += a[pos]* (sm1+divs); *ans %= MOD93;

		a[pos]*(sm1+1) % MOD93
	}
}

// 繰り返し二乗法
fn mod_pow(mut a: usize, mut b: usize, mods: usize) -> usize {
	let mut res = 1;
	while b > 0 {
		if b % 2 == 1 {
			res *= a;
			res %= mods;
		}
		a *= a;
		a %= mods;
		b >>= 1;
	}
	res
}

// 繰り返し二乗法による除算
fn div(a: usize, b: usize, mods: usize) -> usize {
	let mut ans = a;
	let x = mod_pow(b, mods - 2, mods);
	ans *= x;
	ans %= mods;
	ans
}

// const MOD17: usize = 1000000007;
const MOD93: usize = 998244353;
// const INF: usize = 1 << 60;
// let dx = vec![!0, 0, 1, 0]; // 上左下右
// let dy = vec![0, !0, 0, 1]; // 上左下右
// let d = vec!{(!0, 0), (0, !0), (1, 0), (0, 1)}; // 上左下右

#[allow(unused)]
use proconio::{input, marker::Chars, marker::Usize1};

#[allow(unused)]
use std::{
	mem::swap,
	cmp::min, cmp::max,
	cmp::Reverse,
	collections::HashSet, collections::BTreeSet,
	collections::HashMap, collections::BTreeMap,
	collections::BinaryHeap,
	collections::VecDeque,
	iter::FromIterator,
};
0