和果子一起来做题-Project Euler-08-R语言版本

这是第8题:

The four adjacent digits in the 1000-digit number that have the greatest product are 9 × 9 × 8 × 9 = 5832.

mark

Find the thirteen adjacent digits in the 1000-digit number that have the greatest product. What is the value of this product?

这里有1000个连续的数字，如果连续的四个数字相乘，那么最大的乘积是 9 × 9 × 8 × 9 = 5832

如果13个连续的数相乘，那么最大的数值是多少呢？

#复制粘贴这串数字

seq <- "73167176531330624919225119674426574742355349194934

96983520312774506326239578318016984801869478851843

85861560789112949495459501737958331952853208805511

12540698747158523863050715693290963295227443043557

66896648950445244523161731856403098711121722383113

62229893423380308135336276614282806444486645238749

30358907296290491560440772390713810515859307960866

70172427121883998797908792274921901699720888093776

65727333001053367881220235421809751254540594752243

52584907711670556013604839586446706324415722155397

53697817977846174064955149290862569321978468622482

83972241375657056057490261407972968652414535100474

82166370484403199890008895243450658541227588666881

16427171479924442928230863465674813919123162824586

17866458359124566529476545682848912883142607690042

24219022671055626321111109370544217506941658960408

07198403850962455444362981230987879927244284909188

84580156166097919133875499200524063689912560717606

05886116467109405077541002256983155200055935729725

71636269561882670428252483600823257530420752963450"

#去掉换行符

seqnon <- paste(unlist(strsplit(seq,"\n")),collapse = "")

#转换成数字向量

numberlist <- as.numeric(unlist(strsplit(seqnon,"")))

#创建空向量，用于储存计算的乘积

result <- c()

for (i in 1:988){

  result[i] = prod(numberlist[i:(i+12)])

 #prod这个函数可以计算向量内部元素的乘积

}

#找出最大值

max(result)

# 显示为 23514624000

#找出他的位置，就可以得到是哪些数字相乘的

which(result==max(result))

#得到198

numberlist[198:(198+12)]

#是由这13个数字得到的，5 5 7 6 6 8 9 6 6 4 8 9 5

最终结果:

5 X 5 X 7 X 6 X 6 X 8 X 9 X 6 X 6 X 4 X 8 X 9 X 5 =23514624000

[图片上传失败...(image-bdd813-1517135146261)]]

note：

1.which函数用于已知数值返回位置

举例子：

num <- c(5,6,7,34,45,87) 

我想知道45处于什么位置

which(num==45)

返回的是5，验证一下

num[5] 

确实是45，这个太简单了，再看看如果知道数据框或者矩阵中的一个数值，想知道他的位置怎么办呢？

data <- matrix(seq(1:20),nrow = 4)

> data

      [,1] [,2] [,3] [,4] [,5]

[1,]    1    5    9  13  17

[2,]    2    6  10  14  18

[3,]    3    7  11  15  19

[4,]    4    8  12  16  20

想知道15所在的位置

which(data==15,arr.ind = T)

返回的是：

> which(data==15,arr.ind = T)

    row col

[1,]  3  4

在第3行第4列，那就对了。

2.prod函数用于计算向量内数值的乘积

举例子：

num <- c(2,3,4) 

prod(num) 

返回的值就是2,3,4的乘积，24