在 Idris 中实现 Eratosthenes 筛法
Implementing the Sieve of Eratosthenes in Idris
我正在努力将我对埃拉托色尼筛法的定义翻译成伊德里斯。这是到目前为止的功能:
%default total
eratos : Nat -> (l : List Nat) -> { auto ok: NonEmpty l } -> List Nat
eratos limit (prime :: rest) =
if prime * prime > limit -- if we've passed the square root of n
then prime :: xs -- then we're done!
-- otherwise, subtract the multiples of that prime and recurse
else prime :: (eratos limit (rest \ [prime^2,prime^2+prime..limit]))
main : IO ()
main = printLn $ eratos [2..100]
不幸的是,我遇到了一个奇怪的编译器错误:
idris --build euler.ipkg
./E003.idr:18:18: error: expected: ")",
dependent type signature
else prime :: (eratos n (xs \ [prime^2,prime^2+prime..n]))
为什么编译器要查找类型签名?
我能够按如下方式实现它:
eratos : Nat -> (l : List Nat) -> List Nat
eratos _ [] = []
eratos limit (prime :: rest) =
if prime * prime > limit -- if we've passed the square root of n
then prime :: rest -- then we're done!
-- otherwise, subtract the multiples of that prime and recurse
else prime :: eratos limit (rest \ [(prime*prime),(prime*prime+prime)..limit])
通过此实现,类型检查器会考虑此函数 "covering"。
最佳情况下,我们不需要第一种情况,可以将输入限制为列表长度 >= 1 的情况。但是,很难向编译器表明列表永远不会为空或这个函数的第二个参数在每次递归调用时在结构上变得更小。如果有人有建议,请添加为评论或其他答案!
我正在努力将我对埃拉托色尼筛法的定义翻译成伊德里斯。这是到目前为止的功能:
%default total
eratos : Nat -> (l : List Nat) -> { auto ok: NonEmpty l } -> List Nat
eratos limit (prime :: rest) =
if prime * prime > limit -- if we've passed the square root of n
then prime :: xs -- then we're done!
-- otherwise, subtract the multiples of that prime and recurse
else prime :: (eratos limit (rest \ [prime^2,prime^2+prime..limit]))
main : IO ()
main = printLn $ eratos [2..100]
不幸的是,我遇到了一个奇怪的编译器错误:
idris --build euler.ipkg
./E003.idr:18:18: error: expected: ")",
dependent type signature
else prime :: (eratos n (xs \ [prime^2,prime^2+prime..n]))
为什么编译器要查找类型签名?
我能够按如下方式实现它:
eratos : Nat -> (l : List Nat) -> List Nat
eratos _ [] = []
eratos limit (prime :: rest) =
if prime * prime > limit -- if we've passed the square root of n
then prime :: rest -- then we're done!
-- otherwise, subtract the multiples of that prime and recurse
else prime :: eratos limit (rest \ [(prime*prime),(prime*prime+prime)..limit])
通过此实现,类型检查器会考虑此函数 "covering"。
最佳情况下,我们不需要第一种情况,可以将输入限制为列表长度 >= 1 的情况。但是,很难向编译器表明列表永远不会为空或这个函数的第二个参数在每次递归调用时在结构上变得更小。如果有人有建议,请添加为评论或其他答案!