結果

問題 No.526 フィボナッチ数列の第N項をMで割った余りを求める
ユーザー ken miura
提出日時 2021-10-18 21:55:52
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 1 ms / 2,000 ms
コード長 1,091 bytes
コンパイル時間 14,684 ms
コンパイル使用メモリ 379,248 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-19 18:26:02
合計ジャッジ時間 13,630 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff #

use std::io::*;
use std::str::FromStr;

fn read<T: FromStr>() -> T {
  let s = stdin();
  let s = s.lock();
  let s: String = s.bytes()
    .map(|c| c.expect("failed reading char") as char)
    .skip_while(|c| c.is_whitespace())
    .take_while(|c| !c.is_whitespace())
    .collect();
  s.parse().ok().expect("failed parsing")
}

fn mul(
  a: &Vec<Vec<usize>>,
  b: &Vec<Vec<usize>>,
  m: usize,
) -> Vec<Vec<usize>> {
  let mut ret = vec![
    vec![0, 0],
    vec![0, 0],
  ];
  for i in 0..2 {
    for j in 0..2 {
      for k in 0..2 {
        ret[i][j] += a[i][k] * b[k][j];
        ret[i][j] %= m;
      }
    }
  }
  ret
}

fn pow(
  a: Vec<Vec<usize>>,
  n: usize,
  m: usize,
) -> usize {
  let mut a = a.clone();
  let mut n = n-1;
  let mut ret = vec![
    vec![1, 0],
    vec![0, 1],
  ];
  while n > 0 {
    if n & 1 > 0 {
      ret = mul(&a, &ret, m);
    }
    a = mul(&a, &a, m);
    n >>= 1;
  }
  ret[1][0]
}

fn main() {
  let n: usize = read();
  let m: usize = read();
  let a = vec![
    vec![1, 1],
    vec![1, 0],
  ];
  let ans = pow(a, n, m);
  println!("{}", ans);
}
0