結果
問題 | No.1750 ラムドスウイルスの感染拡大-hard |
ユーザー |
|
提出日時 | 2021-11-19 22:44:47 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 246 ms / 2,000 ms |
コード長 | 2,008 bytes |
コンパイル時間 | 11,090 ms |
コンパイル使用メモリ | 401,236 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2025-01-01 20:39:16 |
合計ジャッジ時間 | 14,871 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
// https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8macro_rules! input {($($r:tt)*) => {let stdin = std::io::stdin();let mut bytes = std::io::Read::bytes(std::io::BufReader::new(stdin.lock()));let mut next = move || -> String{bytes.by_ref().map(|r|r.unwrap() as char).skip_while(|c|c.is_whitespace()).take_while(|c|!c.is_whitespace()).collect()};input_inner!{next, $($r)*}};}macro_rules! input_inner {($next:expr) => {};($next:expr,) => {};($next:expr, $var:ident : $t:tt $($r:tt)*) => {let $var = read_value!($next, $t);input_inner!{$next $($r)*}};}macro_rules! read_value {($next:expr, ( $($t:tt),* )) => { ($(read_value!($next, $t)),*) };($next:expr, [ $t:tt ; $len:expr ]) => {(0..$len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()};($next:expr, $t:ty) => ($next().parse::<$t>().expect("Parse error"));}const MOD: i64 = 998_244_353;fn squmul(a: &[Vec<i64>], b: &[Vec<i64>], mo: i64) -> Vec<Vec<i64>> {let n = a.len();let mut ret = vec![vec![0; n]; n];for i in 0..n {for j in 0..n {for k in 0..n {ret[i][k] += a[i][j] * b[j][k];ret[i][k] %= mo;}}}ret}fn squpow(a: &[Vec<i64>], mut e: i64, mo: i64) -> Vec<Vec<i64>> {let n = a.len();let mut sum = vec![vec![0; n]; n];for i in 0..n { sum[i][i] = 1; }let mut cur = a.to_vec();while e > 0 {if e % 2 == 1 {sum = squmul(&sum, &cur, mo);}cur = squmul(&cur, &cur, mo);e /= 2;}sum}fn main() {input! {n: usize, m: usize, t: i64,st: [(usize, usize); m],}let mut mat = vec![vec![0; n]; n];for (s, t) in st {mat[s][t] += 1;mat[t][s] += 1;}let pw = squpow(&mat, t, MOD);println!("{}", pw[0][0]);}