aboutsummaryrefslogtreecommitdiff
path: root/2023/14-Parabolic_Reflector_Dish/second.hs
blob: 2479daa1e8b439f5f4985cce193c0896730e3cf0 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
-- requires cabal install --lib megaparsec parser-combinators
module Main (main) where

import Control.Applicative.Permutations
import Control.Monad (void, when)
import Data.Char qualified as C
import Data.Either
import Data.Functor
import Data.List qualified as L
import Data.Map qualified as M
import Data.Maybe
import Data.Set qualified as S
import Data.Vector qualified as V
import Data.Void (Void)
import Text.Megaparsec
import Text.Megaparsec.Char

import Debug.Trace

exampleExpectedOutput = 64

data Tile = Cube | Empty | Round deriving (Eq, Ord)
instance Show Tile where
  show Cube = "#"
  show Empty = "."
  show Round = "O"
type Row = [Tile]
type Input = [Row]

type Parser = Parsec Void String

parseTile :: Parser Tile
parseTile = char '#' $> Cube
        <|> char '.' $> Empty
        <|> char 'O' $> Round

parseRow :: Parser Row
parseRow = some parseTile <* eol

parseInput' :: Parser Input
parseInput' = some parseRow <* eof

parseInput :: String -> IO Input
parseInput filename = do
  input <- readFile filename
  case runParser parseInput' filename input of
    Left bundle -> error $ errorBundlePretty bundle
    Right input' -> return input'

compute :: Input -> Int
compute input = sum $ map (fst . L.foldr weight (0, 1)) (allPossibilities L.!! theOne)
  where
    transposedInput = L.transpose input
    shift :: Int -> Row -> Row
    shift n [] = replicate n Empty
    shift n (Cube:xs) = replicate n Empty ++ Cube : shift 0 xs
    shift n (Empty:xs) = shift (n+1) xs
    shift n (Round:xs) = Round : shift n xs
    weight :: Tile -> (Int, Int) -> (Int, Int)
    weight Round (acc, i) = (acc + i, i+1)
    weight _ (acc, i) = (acc, i+1)
    theOne = start + (1_000_000_000 - start) `rem` (end - start)
    (start, end) = cycle M.empty 0 allPossibilities
    allPossibilities = iterate process transposedInput
    process = step 4 (L.transpose . map (reverse . shift 0))
    step :: Int -> (a -> a) -> a -> a
    step 0 _ x = x
    step n f x = step (n-1) f $ f x
    cycle :: M.Map Input Int -> Int -> [Input] -> (Int, Int)
    cycle m i (x:xs) = case M.lookup x m of
      Just j -> (j, i)
      Nothing -> cycle (M.insert x i m) (i+1) xs

main :: IO ()
main = do
  example <- parseInput "example"
  let exampleOutput = compute example
  when  (exampleOutput /= exampleExpectedOutput)  (error $ "example failed: got " ++ show exampleOutput ++ " instead of " ++ show exampleExpectedOutput)
  input <- parseInput "input"
  print $ compute input