# Euler problems/71 to 80

### From HaskellWiki

Line 54: | Line 54: | ||

Solution: | Solution: | ||

<haskell> | <haskell> | ||

− | + | import Data.Array (Array, array, (!), elems) | |

+ | import Data.Char (ord) | ||

+ | import Data.List (foldl1') | ||

+ | import Prelude hiding (cycle) | ||

+ | |||

+ | fact :: Integer -> Integer | ||

+ | fact 0 = 1 | ||

+ | fact n = foldl1' (*) [1..n] | ||

+ | |||

+ | factorDigits :: Array Integer Integer | ||

+ | factorDigits = array (0,2177281) [(x,n)|x <- [0..2177281], let n = sum $ map (\y -> fact (toInteger $ ord y - 48)) $ show x] | ||

+ | |||

+ | cycle :: Integer -> Integer | ||

+ | cycle 145 = 1 | ||

+ | cycle 169 = 3 | ||

+ | cycle 363601 = 3 | ||

+ | cycle 1454 = 3 | ||

+ | cycle 871 = 2 | ||

+ | cycle 45361 = 2 | ||

+ | cycle 872 = 2 | ||

+ | cycle 45362 = 2 | ||

+ | cycle _ = 0 | ||

+ | |||

+ | isChainLength :: Integer -> Integer -> Bool | ||

+ | isChainLength len n | ||

+ | | len < 0 = False | ||

+ | | t = isChainLength (len-1) n' | ||

+ | | otherwise = (len - c) == 0 | ||

+ | where | ||

+ | c = cycle n | ||

+ | t = c == 0 | ||

+ | n' = factorDigits ! n | ||

+ | |||

+ | -- | strict version of the maximum function | ||

+ | maximum' :: (Ord a) => [a] -> a | ||

+ | maximum' [] = undefined | ||

+ | maximum' [x] = x | ||

+ | maximum' (a:b:xs) = let m = max a b in m `seq` maximum' (m : xs) | ||

+ | |||

+ | problem_74 :: Int | ||

+ | problem_74 = length $ filter (\(_,b) -> isChainLength 59 b) $ zip ([0..] :: [Integer]) $ take 1000000 $ elems factorDigits | ||

</haskell> | </haskell> | ||

## Revision as of 10:58, 20 August 2007

## Contents |

## 1 Problem 71

Listing reduced proper fractions in ascending order of size.

Solution:

import Data.Ratio (Ratio, (%), numerator) fractions :: [Ratio Integer] fractions = [f | d <- [1..1000000], let n = (d * 3) `div` 7, let f = n%d, f /= 3%7] problem_71 :: Integer problem_71 = numerator $ maximum $ fractions

## 2 Problem 72

How many elements would be contained in the set of reduced proper fractions for d ≤ 1,000,000?

Solution:

Using the Farey Sequence method, the solution is the sum of phi (n) from 1 to 1000000.

See problem 69 for phi function

problem_72 = sum [phi x|x <- [1..1000000]]

## 3 Problem 73

How many fractions lie between 1/3 and 1/2 in a sorted set of reduced proper fractions?

Solution:

import Data.Ratio (Ratio, (%), numerator, denominator) median :: Ratio Int -> Ratio Int -> Ratio Int median a b = ((numerator a) + (numerator b)) % ((denominator a) + (denominator b)) count :: Ratio Int -> Ratio Int -> Int count a b | d > 10000 = 1 | otherwise = count a m + count m b where m = median a b d = denominator m problem_73 :: Int problem_73 = (count (1%3) (1%2)) - 1

## 4 Problem 74

Determine the number of factorial chains that contain exactly sixty non-repeating terms.

Solution:

import Data.Array (Array, array, (!), elems) import Data.Char (ord) import Data.List (foldl1') import Prelude hiding (cycle) fact :: Integer -> Integer fact 0 = 1 fact n = foldl1' (*) [1..n] factorDigits :: Array Integer Integer factorDigits = array (0,2177281) [(x,n)|x <- [0..2177281], let n = sum $ map (\y -> fact (toInteger $ ord y - 48)) $ show x] cycle :: Integer -> Integer cycle 145 = 1 cycle 169 = 3 cycle 363601 = 3 cycle 1454 = 3 cycle 871 = 2 cycle 45361 = 2 cycle 872 = 2 cycle 45362 = 2 cycle _ = 0 isChainLength :: Integer -> Integer -> Bool isChainLength len n | len < 0 = False | t = isChainLength (len-1) n' | otherwise = (len - c) == 0 where c = cycle n t = c == 0 n' = factorDigits ! n -- | strict version of the maximum function maximum' :: (Ord a) => [a] -> a maximum' [] = undefined maximum' [x] = x maximum' (a:b:xs) = let m = max a b in m `seq` maximum' (m : xs) problem_74 :: Int problem_74 = length $ filter (\(_,b) -> isChainLength 59 b) $ zip ([0..] :: [Integer]) $ take 1000000 $ elems factorDigits

## 5 Problem 75

Find the number of different lengths of wire can that can form a right angle triangle in only one way.

Solution: This is only slightly harder than problem 39. The search condition is simpler but the search space is larger.

problem_75 = length . filter ((== 1) . length) $ group perims where perims = sort [scale*p | p <- pTriples, scale <- [1..10^6 `div` p]] pTriples = [p | n <- [1..1000], m <- [n+1..1000], even n || even m, gcd n m == 1, let a = m^2 - n^2, let b = 2*m*n, let c = m^2 + n^2, let p = a + b + c, p <= 10^6]

## 6 Problem 76

How many different ways can one hundred be written as a sum of at least two positive integers?

Solution:

problem_76 = undefined

## 7 Problem 77

What is the first value which can be written as the sum of primes in over five thousand different ways?

Solution:

problem_77 = undefined

## 8 Problem 78

Investigating the number of ways in which coins can be separated into piles.

Solution:

problem_78 = undefined

## 9 Problem 79

By analysing a user's login attempts, can you determine the secret numeric passcode?

Solution:

problem_79 = undefined

## 10 Problem 80

Calculating the digital sum of the decimal digits of irrational square roots.

Solution:

problem_80 = undefined