結果

問題 No.225 文字列変更(medium)
コンテスト
ユーザー こまる
提出日時 2020-10-13 19:59:04
言語 Haskell
(9.10.1)
結果
AC  
実行時間 28 ms / 5,000 ms
コード長 1,753 bytes
コンパイル時間 4,953 ms
コンパイル使用メモリ 203,520 KB
実行使用メモリ 15,872 KB
最終ジャッジ日時 2024-07-20 18:46:02
合計ジャッジ時間 6,155 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます
コンパイルメッセージ
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 # 

{-# LANGUAGE BangPatterns #-}

import           Control.Monad.ST
import           Data.Bool
import qualified Data.Vector.Fusion.Stream.Monadic as VFSM
import qualified Data.Vector.Unboxed               as VU
import qualified Data.Vector.Unboxed.Mutable       as VUM

main :: IO ()
main = do
  [m, n] <- map read . words <$> getLine
  s  <- VU.fromList <$> getLine
  t  <- VU.fromList <$> getLine
  print $ levenshteinDistance m n s t

stream :: Monad m => Int -> Int -> VFSM.Stream m Int
stream !l !r = VFSM.Stream step l
  where
    step x
      | x < r     = return $ VFSM.Yield x (x + 1)
      | otherwise = return $ VFSM.Done
    {-# INLINE [0] step #-}
{-# INLINE [1] stream #-}

rep :: Monad m => Int -> (Int -> m ()) -> m ()
rep n = flip VFSM.mapM_ (stream 0 n)
{-# INLINE rep #-}

rep' :: Monad m => Int -> (Int -> m ()) -> m ()
rep' n = flip VFSM.mapM_ (stream 0 (n + 1))
{-# INLINE rep' #-}

rep1 :: Monad m => Int -> (Int -> m ()) -> m ()
rep1 n = flip VFSM.mapM_ (stream 1 n)
{-# INLINE rep1 #-}

rep1' :: Monad m => Int -> (Int -> m ()) -> m ()
rep1' n = flip VFSM.mapM_ (stream 1 (n + 1))
{-# INLINE rep1' #-}

levenshteinDistance :: Int -> Int -> VU.Vector Char -> VU.Vector Char -> Int
levenshteinDistance m n s t = runST $ do
  dp <- VUM.replicate ((m + 1) * (n + 1)) (1000000007 :: Int)
  rep' m $ \i -> VUM.unsafeWrite dp (i * (n + 1)) i
  rep' n $ \j -> VUM.unsafeWrite dp j j
  rep m $ \i -> rep n $ \j -> do
    item1 <- VUM.unsafeRead dp (i * (n + 1) + j)
    item2 <- VUM.unsafeRead dp (i * (n + 1) + (j + 1))
    item3 <- VUM.unsafeRead dp ((i + 1) * (n + 1) + j)
    VUM.unsafeWrite dp ((i + 1) * (n + 1) + (j + 1)) (minimum [item1 + bool 0 1 (s VU.! i /= t VU.! j), item2 + 1, item3 + 1])
  VUM.unsafeRead dp (m * (n + 1) + n)
0