|
| 1 | +import time |
| 2 | +import numpy as np |
| 3 | + |
| 4 | +#start timer |
| 5 | +start = time.time() |
| 6 | + |
| 7 | +#read the numbers into a grid |
| 8 | +numbers = """08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08 |
| 9 | +49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00 |
| 10 | +81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65 |
| 11 | +52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91 |
| 12 | +22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80 |
| 13 | +24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50 |
| 14 | +32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70 |
| 15 | +67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21 |
| 16 | +24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72 |
| 17 | +21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95 |
| 18 | +78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92 |
| 19 | +16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57 |
| 20 | +86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58 |
| 21 | +19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40 |
| 22 | +04 52 08 83 97 35 99 16 07 97 57 32 16 26 26 79 33 27 98 66 |
| 23 | +88 36 68 87 57 62 20 72 03 46 33 67 46 55 12 32 63 93 53 69 |
| 24 | +04 42 16 73 38 25 39 11 24 94 72 18 08 46 29 32 40 62 76 36 |
| 25 | +20 69 36 41 72 30 23 88 34 62 99 69 82 67 59 85 74 04 36 16 |
| 26 | +20 73 35 29 78 31 90 01 74 31 49 71 48 86 81 16 23 57 05 54 |
| 27 | +01 70 54 71 83 51 54 69 16 92 33 48 61 43 52 01 89 19 67 48""" |
| 28 | + |
| 29 | +#puts the numbers into a vector to be converted into a array |
| 30 | +numbersVector = [] |
| 31 | +for i in range(400): |
| 32 | + i = i*3 |
| 33 | + numbersVector.append(int(numbers[i:i + 2])) |
| 34 | + |
| 35 | +#converts vector to array |
| 36 | +numbersArray = np.array(numbersVector) |
| 37 | +numbersArray.resize((20,20)) |
| 38 | + |
| 39 | +#finds and returns the greatest product of all horizontal products |
| 40 | +def horizontalProductMax(array, length): |
| 41 | + matrixDim = len(array) |
| 42 | + greatestProduct = 0 |
| 43 | + for y in range(matrixDim - (length-1)): |
| 44 | + for x in range(matrixDim): |
| 45 | + tempProduct = 1 |
| 46 | + for number in range(length): |
| 47 | + tempProduct = tempProduct * array[x][y+number] |
| 48 | + if tempProduct > greatestProduct: |
| 49 | + greatestProduct = tempProduct |
| 50 | + return greatestProduct |
| 51 | + |
| 52 | +#finds and returns the greatest product of all vertical products |
| 53 | +def verticalProductMax(array, length): |
| 54 | + matrixDim = len(array) |
| 55 | + greatestProduct = 0 |
| 56 | + for y in range(matrixDim): |
| 57 | + for x in range(matrixDim - (length-1)): |
| 58 | + tempProduct = 1 |
| 59 | + for number in range(length): |
| 60 | + tempProduct = tempProduct * array[x+number][y] |
| 61 | + if tempProduct > greatestProduct: |
| 62 | + greatestProduct = tempProduct |
| 63 | + return greatestProduct |
| 64 | + |
| 65 | +#finds and returns the greatest product of a diagonal product |
| 66 | +def diagonalProductMax(array, length): |
| 67 | + matrixDim = len(array) |
| 68 | + greatestProduct = 0 |
| 69 | + for y in range(matrixDim - (length-1)): |
| 70 | + for x in range(matrixDim - (length-1)): |
| 71 | + tempProduct = 1 |
| 72 | + for number in range(length): |
| 73 | + tempProduct = tempProduct * array[x+number][y+number] |
| 74 | + if tempProduct > greatestProduct: |
| 75 | + greatestProduct = tempProduct |
| 76 | + for y in range(matrixDim - (length - 1)): |
| 77 | + for x in range(matrixDim - (length - 1)): |
| 78 | + tempProduct = 1 |
| 79 | + for number in range(length): |
| 80 | + tempProduct = tempProduct * array[x + length - 1 - number][y + number] |
| 81 | + if tempProduct > greatestProduct: |
| 82 | + greatestProduct = tempProduct |
| 83 | + return greatestProduct |
| 84 | + |
| 85 | +#returns the max value in a list |
| 86 | +def returnMax(list): |
| 87 | + max = list[0] |
| 88 | + for item in list[1:len(list)]: |
| 89 | + if item > max: |
| 90 | + max = item |
| 91 | + return max |
| 92 | + |
| 93 | +maxProduct = returnMax([verticalProductMax(numbersArray, 4), |
| 94 | + horizontalProductMax(numbersArray, 4), |
| 95 | + diagonalProductMax(numbersArray, 4)]) |
| 96 | + |
| 97 | +#end timer |
| 98 | +end = time.time() |
| 99 | + |
| 100 | +#prints results |
| 101 | +print("Time: " + str(end - start)) |
| 102 | +print("Solution: " + str(maxProduct)) |
0 commit comments