結果

問題 No.526 フィボナッチ数列の第N項をMで割った余りを求める
コンテスト
ユーザー かりあげクン
提出日時 2020-09-08 19:44:08
言語 Haskell
(9.10.1)
結果
AC  
実行時間 40 ms / 2,000 ms
コード長 1,318 bytes
コンパイル時間 10,892 ms
コンパイル使用メモリ 197,632 KB
実行使用メモリ 7,668 KB
最終ジャッジ日時 2024-11-30 05:10:03
合計ジャッジ時間 11,741 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます
コンパイルメッセージ
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.Semigroup
import Data.Monoid
import qualified Control.Arrow         as Arrow
import qualified Data.Vector.Unboxed   as VU
import qualified Data.ByteString.Char8 as BSC8

-- https://blog.miz-ar.info/2019/01/fast-fibonacci/

data FibPair a = FibPair !a !a
  deriving (Eq, Show)

instance (Num a) => Semigroup (FibPair a) where
  FibPair a b <> FibPair a' b'
    = FibPair (a * b' + (b - a) * a') (a * a' + b * b')
  stimes = stimesMonoid

instance (Num a) => Monoid (FibPair a) where
  mempty = FibPair 0 1

fibOne :: (Num a) => FibPair a
fibOne = FibPair 1 1

fibPair :: Int -> FibPair Integer
fibPair i = stimesMonoid i fibOne

fastDoubling :: Int -> (Integer, Integer)
fastDoubling 0 = (0, 1)
fastDoubling 1 = (1, 1)
fastDoubling i
  = let (a, b) = fastDoubling (i `quot` 2)
    in if even i
       then (a * (2 * b - a), a * a + b * b)
       else (a * a + b * b, b * (2 * a + b))

fib :: Int -> Integer
fib i = case fastDoubling (i - 1) of
          (a, _) -> a

main :: IO ()
main = do
  (n, mo) <- parse2
  print $ flip mod (fromIntegral mo) $ fib n

type Parser a = BSC8.ByteString -> Maybe (a, BSC8.ByteString)
parseInt :: Parser Int
parseInt = fmap (Arrow.second BSC8.tail) . BSC8.readInt
parse2 :: IO (Int, Int)
parse2 = (\vec -> (vec VU.! 0, vec VU.! 1)) . VU.unfoldrN 2 parseInt <$> BSC8.getLine
0