プロジェクトオイラー問8 - sinapusu2002 @ ウィキ

http://odz.sakura.ne.jp/projecteuler/index.php?cmd=read&page=Problem%208
Problem 8 「数字列中の最大の積」 †
以下の1000桁の数字から13個の連続する数字を取り出して その積を計算する. そのような積の中で最大のものの値はいくらか.
略

1000桁の数字はpe8.txtとして保存
Haskell解

take (length xs-(n-1)) . map (take n) . tails $ xs

は移動平均線などを求めるときのイディオムとして使えるかと。

import System.IO
import System.Directory
import Data.Char
import Data.List

f :: String -> [Int]
f ""=[]
f ('\n':xs)=f xs
f (x:xs)  =(digitToInt x) : (f xs) 

g::Int->[Int]->Int
g n xs=maximum $ map product $ take (length xs-(n-1)) . map (take n) . tails $ xs
main = do
     con<-readFile "pe8.txt"
     putStrLn $ show $ (g 13) $ f con

1000桁の数字はＰｒｏｌｏｇの文字列としてpe8.txtに保存して解。

max(B,B,C):-B>C,!.
max(C,_,C):-!.

mult([],C,0):-C>0,!.
mult(_,0,1):-!.
mult([X|Xs],C,Result):-
	C1 is C-1,
	mult(Xs,C1,Re),
	Result is Re*(X-48).
search([],Ans):-
	!,
	write(Ans).
 

search([X|Xs],Ans):-
	mult([X|Xs],13,Ans1),
 	max(Ans2,Ans,Ans1),
	search(Xs,Ans2).

main:-
	see('pe8.txt'),
	read(X),
	seen,
	search(X,0).

下記はアペンドを使った場合。
データはソースに直書き。
実はＰｒｏｌｏｇ素人とプロの差の第一歩はアペンドを自由自在に使えるか使えないかだったりする。
今回のコードの注目点はここ
findall(E,(append(_,B,A),calc(B,0,1,E)),Es),
手前から積のスタート地点を一つずつずらしているがＢは参照に過ぎないので処理的に言って何の負荷もない
ＣＰＵ的には軽い処理なのである。
まあfindallしちゃってるのはメモリ的には重い処理なのでそこは考え物です。

calcが１３で終わるというのは汎用性に乏しい。
積の長さは述語の引数として渡せるほうが柔軟だな、そこは下記コードは失敗している。
失敗した探索はコードが勝手に捨ててくれる感覚がわかるとＰｒｏｌｏｇが柔軟に動くようになる。

a("7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450").

calc(_,13,M,M):-!.
calc([],_,_,_):-!,fail.
calc([X|Xs],L,M,Result):-
	M1 is M*(X-48),
	L1 is L+1,
	calc(Xs,L1,M1,Result).

main:-
	a(A),
	findall(E,(append(_,B,A),calc(B,0,1,E)),Es),
	sort(Es,Es1),
	append(_,[Ans],Es1),
	write(Ans).