2016年1月8日金曜日

数列中での最大の積

Project Euler 8

次の1000桁の中で隣り合う4つの整数の積のうち最大の値は、$9\times 9\times 8\times 9 = 5832$ である。

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

では、1000桁の中で隣り合う数の積が最大になる13桁の数は何か?
また、その値は?


先頭から13桁を取得し、1桁ずつ乗算する。
なお、計算量を削減するため、数字の中に「0」が含まれる場合は乗算を中止する。
# SAMPLE_1

N = 13
max_num = 0
max_str = ""

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



def prod(sx):
    x = 1
    if "0" in sx:
        return 0
    for i in range(len(sx)):
        x *= int(sx[i])
    return x

x_str = s[:N-1]
for i in xrange(N-1, len(s)):
    x_str += s[i]
    x = prod(x_str)
    if max_num < x:
        max_num = x
        max_str = x_str
        print max_num, max_str
    x_str = x_str[1:]

0 件のコメント:

コメントを投稿