看到各位对“笨办法系列”的东西还比较感兴趣,我也很乐意继续写下去。今天的示例将会用到数据框(data.frame)这种数据类型,并学习如何组合计算两个向量,以及如何排序。我们将用所学的东西来解决Project Euler的,就是找出一个集合中最大的回文数。回文数是指一个像1534351这样“对称”的数,如果将这个数的数字按相反的顺序重新排列后,所得到的数和原来的数一样。开始啦! 得到的结果是906609,本例是先将乘积排序后再判断是否回文数,找到的第一个就是答案,所以速度会快一点。如果不用expand.grid函数的话,可以利用嵌套for来组合计算。另外还有一种作法是利用R本身的rev函数,先将数字转为字符,再切开成一串向量,用rev反转后判断。 from 数据科学与R语言 http://xccds1977.blogspot.com/
# 预备练习 x <- y <- 1:9 <- (x=x,y=y) () z <- $x * $y # 一个九九乘法表 z <- (z,=9) (1) x <- ((10),2) (x) (x) x[(x)[1]] (x) x[(x)] x[(x)] y <- 1:10 <- (x,y) () () [1,] [,1] $x [($x),] # Project Euler 4 # 在两个三位数字的乘积中,找出最大的回文数 # 先建立一个将数字顺序进行反转的函数 reverse <- (n) { reversed <- 0 while (n > 0) { reversed <- 10 * reversed + n %% 10 n <- n%/%10 } (reversed) } # 从大到小搜索回文数
x <- y <- 999:100 <- (x=x,y=y) $prod <- $x * $y <- [($prod,decreasing=T),] ()
value <- $prod for (i in 1:(value)) { isequal <- (value[i] == reverse(value[i])) if (isequal) { ([i,]) break }}