結果

問題 No.1112 冥界の音楽
ユーザー StrorkisStrorkis
提出日時 2020-07-11 13:41:44
言語 Rust
(1.77.0)
結果
AC  
実行時間 28 ms / 2,000 ms
コード長 2,316 bytes
コンパイル時間 3,599 ms
コンパイル使用メモリ 147,892 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-03 02:23:52
合計ジャッジ時間 2,739 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,384 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 2 ms
4,384 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 1 ms
4,376 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 1 ms
4,380 KB
testcase_11 AC 2 ms
4,384 KB
testcase_12 AC 1 ms
4,380 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 8 ms
4,380 KB
testcase_15 AC 3 ms
4,380 KB
testcase_16 AC 1 ms
4,376 KB
testcase_17 AC 4 ms
4,376 KB
testcase_18 AC 3 ms
4,376 KB
testcase_19 AC 1 ms
4,380 KB
testcase_20 AC 1 ms
4,380 KB
testcase_21 AC 8 ms
4,384 KB
testcase_22 AC 7 ms
4,380 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 1 ms
4,380 KB
testcase_25 AC 1 ms
4,376 KB
testcase_26 AC 1 ms
4,384 KB
testcase_27 AC 1 ms
4,376 KB
testcase_28 AC 2 ms
4,380 KB
testcase_29 AC 1 ms
4,380 KB
testcase_30 AC 8 ms
4,376 KB
testcase_31 AC 1 ms
4,384 KB
testcase_32 AC 10 ms
4,380 KB
testcase_33 AC 1 ms
4,380 KB
testcase_34 AC 1 ms
4,380 KB
testcase_35 AC 1 ms
4,376 KB
testcase_36 AC 28 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

use std::ops::{Mul, MulAssign};

#[derive(Clone)]
struct Matrix(Vec<Vec<usize>>);

impl Matrix {
    const MOD: usize = 1_000_000_007;

    fn pow(self, mut n: usize) -> Matrix {
        let p = self.0.len();

        let mut a = vec![vec![0; p]; p];
        for i in 0..p {
            a[i][i] = 1;
        }

        let mut res = Matrix(a);
        let mut x = self.clone();

        while n > 0 {
            if n & 1 == 1 { res *= x.clone(); }
            x *= x.clone();
            n >>= 1;
        }

        res
    }
}

impl Mul for Matrix {
    type Output = Matrix;

    fn mul(self, other: Matrix) -> Matrix {
        let p = self.0.len();
        let q = other.0[0].len();
        let r = other.0.len();

        let mut a = vec![vec![0; q]; p];

        for i in 0..p {
            for j in 0..q {
                for k in 0..r {
                    a[i][j] += self.0[i][k] * other.0[k][j] % Matrix::MOD;
                    a[i][j] %= Matrix::MOD;
                }
            }
        }

        Matrix(a)
    }
}

impl MulAssign for Matrix {
    fn mul_assign(&mut self, rhs: Matrix) {
        *self = (*self).clone() * rhs;
    }
}

fn main() {
    const MOD: usize = 1_000_000_007;

    let (k, m, n): (usize, usize, usize) = {
        let mut buf = String::new();
        std::io::stdin().read_line(&mut buf).unwrap();
        let mut iter = buf.split_whitespace();
        (
            iter.next().unwrap().parse().unwrap(),
            iter.next().unwrap().parse().unwrap(),
            iter.next().unwrap().parse().unwrap(),
        )
    };

    let mut a = Matrix(vec![vec![0; k * k]; k * k]);
    for _ in 0..m {
        let (p, q, r): (usize, usize, usize) = {
            let mut buf = String::new();
            std::io::stdin().read_line(&mut buf).unwrap();
            let mut iter = buf.split_whitespace();
            (
                iter.next().unwrap().parse::<usize>().unwrap() - 1,
                iter.next().unwrap().parse::<usize>().unwrap() - 1,
                iter.next().unwrap().parse::<usize>().unwrap() - 1,
            )
        };
        a.0[p * k + q][q * k + r] = 1;
    }
    let res = a.pow(n - 2);

    let mut ans = 0;
    for i in 0..k {
        for j in 0..k {
            ans = (ans + res.0[i][j * k]) % MOD;
        }
    }
    println!("{}", ans);
}
0