Project Euler No. 1 ~ No.10 を解いた

No.1 ~ 10までのコードをペタペタしていきます.

No.1 Multiples of 3 and 5

1000以下のxについて,x mod 5 == 0, x mod 3 == 0の数の合計を求める

check :: Int -> [Int]
check a =
  if (a `mod` 3 == 0 || a `mod` 5 == 0)
  then a : (check $ a + 1)
  else check $ a + 1

main = do
  print $ sum . takeWhile (\x -> x < 1000) $ check 1

やってることは,checkで無限リストを作成して,takeWhileで1,000以下のものを取ってくる.その後,合計を求めて終わりです.

No.2 Even Fibonacci numbers

フィボナッチから偶数を取り出して,合計を求める.

main = do
  let fibo = 0 : 1 : zipWith (\x y -> x + y) (tail fibo) fibo
  print $ sum $ filter even $ takeWhile (\x -> x < 4000000) fibo

fiboというフィボナッチの無限リストを作成し,それからtakeWhileしてsumするだけ.

No.3 Largest prime factor

600851475143の最大の素因数を求める.

factors :: Integer -> [Integer]
factors n = [x | x <- [1..n], n `mod` x == 0]

factorization :: Integer -> [Integer]
factorization 1 = []
factorization x = v : factorization (x `div` v)
  where
    v = (factors x) !! 1

main = do
  print $ factorization 600851475143

以下のブログのものを利用しました.ありがとうございます🙇‍♂️

tnomura9.exblog.jp

No.4 Largest palindrome product

100..999までの数字をかけて,回文となる数値を求める.

reverseInt :: Integer -> Integer
reverseInt x = read . reverse . show $ x

check :: Integer -> Integer
check a = if a == reverseInt a then a else 0

main = do
  let initList = [100,101..999]
  print . maximum . fmap check $ (*) <$> initList <*> initList

数字を文字列にして反転して数値化し,元の数値と比較する.等しければ,その数値を返す.あとはmaxとって終わり.

No.5 Smallest multiple

1..20までの最大公約数を求める問題.

check :: [Integer] -> Integer
check [] = 1
check (x:xs) = lcm (x) (check xs)

main = do
  print $ check [1..20]

Listを回して,lcm(最大公約数を求める関数)するだけ.

No. 6 Sum square difference

(a ^ 2 + b ^ 2... ) - (a + b ...) ^ 2 の解を求める.範囲は1..100

main = do
  let sumNum = sum [1..100]
  let a = (*) sumNum sumNum
  let b = sum $ fmap (\x -> x * x) [1..100]
  print $ a - b

aが各要素を足したあと最後に二乗した数値.そしてbが各要素を二乗していき,最後に合計を求めた数値.最後にa - bをして終わり.

No.7 10001st prime

10001番目の素数を求める.

import Data.Numbers.Primes
main = do
  print . last . take 10001 $ wheelSieve 1

素数を求める時には,Data.Numbers.Primesパッケージを利用しました.

インストールする必要があるので,stack install primesでインストールしてください.

No.8 Largest product in a series

与えられた数値の隣り合う13桁の数を取り出していき,その13個の数の積の最大を求める.

calc :: [Integer] -> Integer
calc (a:b:c:d:e:f:g:h:i:j:k:l:m:[]) = a * b * c * d * e * f * g * h * i * j * k * l * m
calc otherwise = 0

check :: [String] -> [Integer]
check xss@(a:b:c:d:e:f:g:h:i:j:k:l:m:xs) = (calc $ fmap read [a,b,c,d,e,f,g,h,i,j,k,l,m]) : (check $ tail xss)
check otherwise = []

main = do
  print $ maximum $ check $ fmap (:[]) $ show 7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450

なかなかゴリ押ししている感じがあったりします.まず与えられた数値をリストにして,一つずつ右にずらしていき13個ずつ取って行く. あとはそれらの最大をとって終わり.

No.9 Special Pythagorean triplet

ピタゴラスの定理が成り立つa, b, cをとり,a + b + c = 1000になるものを求める.

main = do
  print $ [z | z <- (\x y -> let c = sqrt ((x * x) + (y * y))
      in if x < y && x + y + c == 1000.0
        then Just (x,y,c)
        else Nothing) <$> [1..999] <*> [1..999], not (z == Nothing)] !! 0

ゴリ押しました.[1..999]のリストを二つ与えて,条件にあったもののみをJustに包んで返します.そしてリスト内包表記でNothingを弾きます.(多分もっといい方法ありそう)

No.10 Summation of primes

2,000,000以下の素数の合計を求める.

import Data.Numbers.Primes
main = do
  print $ sum $ takeWhile (\x -> x < 2000000) $ wheelSieve 1

特に言うことがない

まとめ

Haskellの練習がてら少しやったのですが,Haskellのすごさを色々と思い知り感動させてもらいました.また隙間時間に解いて行きます.🙇‍♂️