結果

問題 No.526 フィボナッチ数列の第N項をMで割った余りを求める
ユーザー HaarHaar
提出日時 2018-09-26 09:47:53
言語 Haskell
(9.8.2)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 506 bytes
コンパイル時間 2,107 ms
コンパイル使用メモリ 171,136 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-09 10:48:35
合計ジャッジ時間 2,861 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,820 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 2 ms
6,820 KB
testcase_03 AC 1 ms
6,820 KB
testcase_04 AC 1 ms
6,820 KB
testcase_05 AC 2 ms
6,820 KB
testcase_06 AC 1 ms
6,820 KB
testcase_07 AC 1 ms
6,816 KB
testcase_08 AC 2 ms
6,820 KB
testcase_09 AC 1 ms
6,820 KB
testcase_10 AC 2 ms
6,820 KB
testcase_11 AC 2 ms
6,820 KB
testcase_12 AC 1 ms
6,820 KB
testcase_13 AC 2 ms
6,816 KB
testcase_14 AC 2 ms
6,820 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Loaded package environment from /home/judge/.ghc/x86_64-linux-9.8.2/environments/default
[1 of 2] Compiling Main             ( Main.hs, Main.o )
[2 of 2] Linking a.out

ソースコード

diff # 

import Data.Int
import Data.Array

main = do
  [n,m] <- map read . words <$> getLine :: IO [Int64]
  print $ fibonacci (n-1) m

{-- Fibonacci数の一般項 --}
fibonacci 0 _ = 0
fibonacci 1 _ = 1
fibonacci n p = let [x,_,_,_] = matpow [1,1,1,0] (n-1) in x `mod` p
  where
    matprod [a,b,c,d] [w,x,y,z] = map (`mod` p) [a*w+b*y,a*x+b*z,c*w+d*y,c*x+d*z]
    matpow a 0 = [1,0,0,1]
    matpow a 1 = a
    matpow a p = let t = matpow a (div p 2) in (if even p then id else (`matprod` a)) (t `matprod` t)
0