何かを書き留める何か

数学や読んだ本について書く何かです。最近は社会人として生き残りの術を学ぶ日々です。

Project Euler Problem 8

Find the greatest product of five consecutive digits in the 1000-digit number.

73167176531330624919225119674426574742355349194934
96983520312774506326239578318016984801869478851843
85861560789112949495459501737958331952853208805511
12540698747158523863050715693290963295227443043557
66896648950445244523161731856403098711121722383113
62229893423380308135336276614282806444486645238749
30358907296290491560440772390713810515859307960866
70172427121883998797908792274921901699720888093776
65727333001053367881220235421809751254540594752243
52584907711670556013604839586446706324415722155397
53697817977846174064955149290862569321978468622482
83972241375657056057490261407972968652414535100474
82166370484403199890008895243450658541227588666881
16427171479924442928230863465674813919123162824586
17866458359124566529476545682848912883142607690042
24219022671055626321111109370544217506941658960408
07198403850962455444362981230987879927244284909188
84580156166097919133875499200524063689912560717606
05886116467109405077541002256983155200055935729725
71636269561882670428252483600823257530420752963450

http://projecteuler.net/problem=8

 

1000桁の数字から連続する5つの数を取り出してその積を計算し、その中で最大値を求める問題。

今まではSageの組み込み函数で一発で計算できるというケースしかなかったが、今回はPythonの標準モジュールを用いてみることにする。用いるのは次の2つ:

 

operator 函数形式の標準演算子

os 雑多なオペレーティングシステムインターフェース

http://docs.python.jp/2.6/library/index.html を参照した。

このうち、osはファイルを取り込むためにディレクトリを指定するために用いたためプログラムにはそれほど関係が無い。operatorはsum函数の掛け算版を作るために用いた。

----

import operator
import os

def prod(lst):
    return reduce(operator.mul, lst)
os.chdir("PATH")


target = [int(line) for line in open('Problem_008.txt', 'r')]
target = map(int,list("".join(map(str,target))))
print max([prod(target[i:i+5]) for i in xrange(len(target)) if len(target[i:i+5]) ==5])
----

 

最初に書いたときはosモジュールの存在を知らず一々手動でリストにappendで取り込んでいた。20行ならそれほどでもないが行数が多くなるとすぐ破綻する。この記事を書くために相当ブラッシュアップしたがまだまだ改善の余地はあるかもしれない。