Question

我正在努力让中国剩余定理算法工作，所以我一直在网上寻找帮助。我试图在haskell中得到这个CRT的例子来编译，但是我遇到了这些错误。我已经实现了自己的extGCD功能。

extGCD :: Integer -> Integer -> (Integer, Integer, Integer)
extGCD a 0 = (a, 1, 0)
extGCD a b = let (q, r) = divMod a b
             (d, m, n) = extGCD b r
         in (d, n, m - n*q)

crt :: [Integer] -> [Integer] -> Integer
crt as ms = let {p = product ms
            ;ls = [extGCD x (div p x) !! 1 |x<- ms]
        }in sum [ div (x*y*p) z | (x,y,z)<- zip3 as ls ms ]

这是错误：

 Couldn't match expected type `[t0]'
            with actual type `(Integer, Integer, Integer)'
    In the return type of a call of `extGCD'
    In the first argument of `(!!)', namely `extGCD x (div p x)'
    In the expression: extGCD x (div p x) !! 1
Failed, modules loaded: none.

Answer 1

这不是Python。列表和元组不是同一类型。甚至没有关闭。

这就是错误消息告诉你的内容。 (Integer, Integer, Integer)不是列表。因此，您无法将!!运算符应用于extGCD的返回。

base包不包含使用三元组的函数，因此我可能会将列表推导更改为[ x' | x <- ms, let (_, x', _) = extGCD x (div p x) ]。

Haskell中的CRT实现

1 个答案: