项目Euler 8 - Haskell

时间:2014-08-16 21:32:34

标签: haskell

浏览Project Euler我正在将我的解决方案与here进行比较。

对于问题8,我的代码会产生正确的答案(通过网站上的支票金额确认)23514624000。

module Main where

import Data.List

main = do
    print $ last (sort eulerEight)


eulerEight = doCalc [ x | x <- toDigits 7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450]
    where
        doCalc (_:[]) = []
        doCalc (x:xs) = if 0 `notElem` take 13 (x:xs) && length (x:xs) > 12
                then product (take 13 (x:xs)) : doCalc xs
                else doCalc xs

toDigits n 
 | n < 1 = []
 | otherwise = toDigits (n `div` 10) ++ [n `mod` 10]

我意识到,检查解决方案here可能要好得多,但这似乎不正确。

import Data.Char (digitToInt)
import Data.List

problem_8 = do
        str <- readFile "number.txt"
        -- This line just converts our str(ing) to a list of 1000 Ints
        let number = map digitToInt (concat $ lines str)
        print $ maximum $ map (product . take 13) (tails number)

我已根据Project Euler网站上的问题将take 5的值更改为13,但上面的代码产生的错误答案为2091059712.我已经检查了number.txt中的数字是否正确它有两个例子的完整1000位数字。任何人都可以阐明产出不同的原因吗? (我想可能与它使用tails而不是tail的事实有关,但我不确定)

2 个答案:

答案 0 :(得分:6)

正在发生的事情是digitToInt返回Int,这在32位系统上太短,无法在5增加到13时保存测试数字。将其更改为(fromIntegral . digitToInt)它工作正常。

答案 1 :(得分:2)

问题已被识别为Int溢出,但维基代码本身不正确。它没有正确地截断列表,并且可能产生不正确的结果,这取决于输入数据(即,它通过幸运的机会产生正确的结果)。

想象一个数字字符串,以0后跟12 9 s结尾。在计算最大值时,代码将错误地考虑9^12。更简单的是,对于1000个零,它将产生1作为答案。

由于压缩的属性,我们可以实现自动截断:

import Data.Char
import Data.List

euler_8 = do
   str <- readFile "numbers.txt"
   print . maximum . map product
         . foldr (zipWith (:)) (repeat [])   -- here
         . take 13 . tails . map (fromIntegral . digitToInt)
         . concat . lines $ str

您的代码虽然 正确,但有一些问题:

  • [x | x <- xs]只是xs
  • last (sort xs)只是maximum xs,更快
  • 附加在递归调用右侧的
  • 是效率低下的已知来源。最好甚至在左边追加并在最后反转,但下面的转换更多是哈斯克利安:
toDigits n xs  -- to be called as `toDigits n []`
 | n < 1 = xs
 | otherwise = toDigits (n `div` 10) ((n `mod` 10) : xs)