Fortran Implementation of Problem 8

View source code here on GitHub!

integer Problem0008/p0008()
 1! Project Euler Problem 8
 2!
 3! This was easier to do in Fortran than I would have thought
 4!
 5! Problem:
 6!
 7! The four adjacent digits in the 1000-digit number that have the greatest product are 9 × 9 × 8 × 9 = 5832.
 8!
 9! 73167176531330624919225119674426574742355349194934
10! 96983520312774506326239578318016984801869478851843
11! 85861560789112949495459501737958331952853208805511
12! 12540698747158523863050715693290963295227443043557
13! 66896648950445244523161731856403098711121722383113
14! 62229893423380308135336276614282806444486645238749
15! 30358907296290491560440772390713810515859307960866
16! 70172427121883998797908792274921901699720888093776
17! 65727333001053367881220235421809751254540594752243
18! 52584907711670556013604839586446706324415722155397
19! 53697817977846174064955149290862569321978468622482
20! 83972241375657056057490261407972968652414535100474
21! 82166370484403199890008895243450658541227588666881
22! 16427171479924442928230863465674813919123162824586
23! 17866458359124566529476545682848912883142607690042
24! 24219022671055626321111109370544217506941658960408
25! 07198403850962455444362981230987879927244284909188
26! 84580156166097919133875499200524063689912560717606
27! 05886116467109405077541002256983155200055935729725
28! 71636269561882670428252483600823257530420752963450
29!
30! Find the thirteen adjacent digits in the 1000-digit number that have the greatest product. What is the value of this product?
31
32module Problem0008
33    use constants
34    implicit none
35contains
36    pure integer(i18t) function p0008() result(answer)
37        integer :: i, j
38        integer(i18t) :: tmp
39        integer(i1t), dimension(1000) :: digits
40        character(len=1000) :: plain_digits
41        plain_digits = ( &
42            "73167176531330624919225119674426574742355349194934" // &
43            "96983520312774506326239578318016984801869478851843" // &
44            "85861560789112949495459501737958331952853208805511" // &
45            "12540698747158523863050715693290963295227443043557" // &
46            "66896648950445244523161731856403098711121722383113" // &
47            "62229893423380308135336276614282806444486645238749" // &
48            "30358907296290491560440772390713810515859307960866" // &
49            "70172427121883998797908792274921901699720888093776" // &
50            "65727333001053367881220235421809751254540594752243" // &
51            "52584907711670556013604839586446706324415722155397" // &
52            "53697817977846174064955149290862569321978468622482" // &
53            "83972241375657056057490261407972968652414535100474" // &
54            "82166370484403199890008895243450658541227588666881" // &
55            "16427171479924442928230863465674813919123162824586" // &
56            "17866458359124566529476545682848912883142607690042" // &
57            "24219022671055626321111109370544217506941658960408" // &
58            "07198403850962455444362981230987879927244284909188" // &
59            "84580156166097919133875499200524063689912560717606" // &
60            "05886116467109405077541002256983155200055935729725" // &
61            "71636269561882670428252483600823257530420752963450")
62        answer = 0
63        do i = 1, 1000
64            digits(i) = int(ichar(plain_digits(i:i)) - ichar('0'), kind=1)
65        end do
66        do i = 1, (1000 - 12)
67            tmp = digits(i)
68            do j = (i + 1), (i + 12)
69                tmp = tmp * digits(j)
70            end do
71            answer = max(answer, tmp)
72        end do
73    end function p0008
74end module Problem0008