--/--/--

スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

2008/10/10

Project Euler: Problem 11 修正

Project Euler: Problem 11 - たろうbぉg が間違っていた。縦は下方向、横は右方向、斜めは右下方向だけ調べていたのだが、左下方向を調べていなかった。ので修正した。結果も変わったようである。

また、修正するにあたり、高階関数を使うようにプログラムを変えた。余計見にくくなった感も否めない。

実行結果

$ runhaskell 11.hs
70600674

プログラム

 1 input = [08,02,22,97,38,15,00,40,00,75,04,05,07,78,52,12,50,77,91,08,
 2          49,49,99,40,17,81,18,57,60,87,17,40,98,43,69,48,04,56,62,00,
 3          81,49,31,73,55,79,14,29,93,71,40,67,53,88,30,03,49,13,36,65,
 4          52,70,95,23,04,60,11,42,69,24,68,56,01,32,56,71,37,02,36,91,
 5          22,31,16,71,51,67,63,89,41,92,36,54,22,40,40,28,66,33,13,80,
 6          24,47,32,60,99,03,45,02,44,75,33,53,78,36,84,20,35,17,12,50,
 7          32,98,81,28,64,23,67,10,26,38,40,67,59,54,70,66,18,38,64,70,
 8          67,26,20,68,02,62,12,20,95,63,94,39,63,08,40,91,66,49,94,21,
 9          24,55,58,05,66,73,99,26,97,17,78,78,96,83,14,88,34,89,63,72,
10          21,36,23,09,75,00,76,44,20,45,35,14,00,61,33,97,34,31,33,95,
11          78,17,53,28,22,75,31,67,15,94,03,80,04,62,16,14,09,53,56,92,
12          16,39,05,42,96,35,31,47,55,58,88,24,00,17,54,24,36,29,85,57,
13          86,56,00,48,35,71,89,07,05,44,44,37,44,60,21,58,51,54,17,58,
14          19,80,81,68,05,94,47,69,28,73,92,13,86,52,17,77,04,89,55,40,
15          04,52,08,83,97,35,99,16,07,97,57,32,16,26,26,79,33,27,98,66,
16          88,36,68,87,57,62,20,72,03,46,33,67,46,55,12,32,63,93,53,69,
17          04,42,16,73,38,25,39,11,24,94,72,18,08,46,29,32,40,62,76,36,
18          20,69,36,41,72,30,23,88,34,62,99,69,82,67,59,85,74,04,36,16,
19          20,73,35,29,78,31,90,01,74,31,49,71,48,86,81,16,23,57,05,54,
20          01,70,54,71,83,51,54,69,16,92,33,48,61,43,52,01,89,19,67,48]
21
22 tateMax values = fn values
23                     (\p len -> (p+60>= len)
24                     (\_     -> False)
25                     (+0)
26                     20
27
28 yokoMax values = fn values
29                     (\p len -> (p+3> len)
30                     (\p     -> (p+3`mod` 20 == 0)
31                     (+3)
32                     1
33
34 migiNanameMax values = fn values
35                        (\p len -> (p+63> len)
36                        (\p     -> (p+3`mod` 20 == 0)
37                        (+3)
38                        21
39
40 hidariNanameMax values = fn values
41                             (\p len -> (p+60>= len)
42                             (\p     -> p `mod` 20 == 0)
43                             (+3)
44                             19
45
46 fn values finish skip skipNext step = f 0 0 (length values)
47     where f r p len | finish p len  = r
48                     | skip p        = f r (skipNext p) len
49                     | otherwise     = let n = product' p step values
50                                           in if n > r
51                                                 then f n (p+1) len
52                                                 else f r (p+1) len
53
54 product' p i values = (values !! p) *
55                       (values !! (p+i)) *
56                       (values !! (p+(i*2))) *
57                       (values !! (p+(i*3)))
58
59 main = print $ maximum $ map (\-> f $ input)
60                              [tateMax, yokoMax,
61                               migiNanameMax, hidariNanameMax]
スポンサーサイト

comment

post




上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。