1
0
Fork 0
projecteuler/69.hs

12 lines
391 B
Haskell

import Math.NumberTheory.Primes.Factorisation
import Data.List
ratio :: Integer -> Double
ratio n = foldr (\(prime, _) acc -> let p = fromIntegral prime in p / (p - 1) * acc) 1.0 $ factorise n
answer :: Integer -> (Integer, Double)
answer n = let range = [2, 4..n]
in maximumBy (\(i1, r1) (i2, r2) -> compare r1 r2) $ zip range $ map ratio range
main = print $ answer 1000000