角谷定理
【题目】给定一个正整数 N, 如果是偶数就折半,如时奇数就 3 * N + 1,如此往复运算,最后一定得到 1。我们把一个数 N 到达 1 用了多少步运算称为它的角谷步长。
求 100 以内的数字角谷步长的最大长度。
例如,N=7时,它的运算过程:22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1
所以,它的角谷步长是17
这种折腾序列的问题是 haskell 的强项,表达起来灵活、强大。
代码:
jiao :: Int -> [Int]
jiao 1 = [1]
jiao n | even n = n : jiao (n `div` 2)
| otherwise = n : jiao (n * 3 + 1)
jiao_len :: Int -> Int
jiao_len n = maximum [length (jiao x) | x<-[1..n]] - 1
main = print $ jiao_len 100
代码十分简单、易读,就不画蛇填足解释了。
jiao 这个函数解决了根本问题,给一个n,它就把运算序列都给出来了,以后随便怎么加工都可以的。
比如,我们想知道是谁贡献了最大角谷序列。
maximumBy (compare `on` length) [jiao x | x<-[1..100]]
一行代码就解决问题了。
注意,maximumBy 在Data.List中, on 在 Data.Function 中。
on 用了中缀(infix)表达形式。如果不加反引号,普通函数默认都是前缀(prefix)形式。