用欧拉计划学Rust语言(第17~21题)

方宏富
2023-12-01

最近想学习Libra数字货币的MOVE语言,发现它是用Rust编写的,所以先补一下Rust的基础知识。学习了一段时间,发现Rust的学习曲线非常陡峭,不过仍有快速入门的办法。

学习任何一项技能最怕没有反馈,尤其是学英语、学编程的时候,一定要“用”,学习编程时有一个非常有用的网站,它就是“欧拉计划”,网址: About - Project Euler

这个网站提供了几百道由易到难的数学问题,你可以用任何办法去解决它,当然主要还得靠编程,编程语言不限,论坛里已经有Java、C#、Python、Lisp、Haskell等各种解法,当然如果你直接用google搜索答案就没任何乐趣了。

学习Rust最好先把基本的语法和特性看过一遍,然后就可以动手解题了,解题的过程就是学习、试错、再学习、掌握和巩固的过程,学习进度会大大加快。

第17题

问题描述:

1到1000用英文单词写下来,求总字符个数(空格和连字符不算),例如:342,英文单词是:three hundred and forty-two。

问题分解:

  • 数字转换成英文单词
    • 1到19的拼写
    • 20到99的拼写
    • 100到999的拼写
    • 1000的拼写
  • 单词中去掉空格和连字符
  • 取字符总数

1到19的拼写比较特殊,需要分别对待,而超过20的数,可以利用递归调用。这里可以学到String的语法知识点。Rust中的字符串有点烦人,list[n].to_string()、"one thousand".to_string()的这种写法让人非常不适应。除了String之外,还有字符串切片(slice)、字符常量,需要看基础教程慢慢理解。

fn english_number(n: usize) -> String {
    let list0_9 = vec![
        "zero", "one", "two", "three", "four", 
        "five", "six", "seven", "eight", "nine",
    ];
    if n <= 9 {
        return list0_9[n].to_string();
    }
    if n <= 19 {
        let list = vec![
            "ten", "eleven", "twelve", "thirteen", "fourteen",
            "fifteen", "sixteen", "seventeen", "eighteen", "nineteen",
        ];
        return list[n - 10].to_string();
    }
    if n <= 99 {
        let a: usize = n / 10; // 十位
        let b: usize = n % 10;
        let list = vec![
            "", "", "twenty", "thirty", "forty", 
            "fifty", "sixty", "seventy", "eighty", "ninety"
        ];
        let str = list[a].to_string();
        if b > 0 {
            return str + "-" + &english_number(b);
        }
        return str;
    }
    if n <= 999 {
        let a: usize = n / 100; // 百位
        let b: usize = n % 100;
        let str = list0_9[a].to_string() + " hundred";
        if b > 0 {
            return str + " and " + &english_number(b);
        }
        return str;
    }
    if n == 1000 {
        return "one thousand".to_string();
    }
    return "unknown".to_string();
}

从字符串里移除特定的字符,要利用函数式编程,还有filter()和collect()函数,一气呵成。filter()函数中的*c又是让人写错的地方。

fn remove_space(s: &str) -> String {
    s.chars().filter(|c| *c != ' ' && *c != '-').collect()
}

主程序就比较容易了,求和即可。

let mut sum = 0;
for n in 1..=1000 {
    let s = remove_space(&english_number(n));
    sum += s.len();
    // println!("{}: {} {}", n, english_number(n), s.len());
}
println!("{}", sum);

第18题

问题描述:

从堆成三角的数字中,找到一条路径,使其和最大,求和。一个节点的下一个点只能是下一层的左、右节点。

为了节省内存空间,用一维数组表示这些数,需要准确地计算出各个索引位置的行号,我为了方便,最上一层的行号为1。

let w = [
                                75, // row 1
                              95, 64,  // row 2
                            17, 47, 82,  // row 3
                          18, 35, 87, 10,
                        20, 04, 82, 47, 65,
                      19, 01, 23, 75, 03, 34,
                    88, 02, 77, 73, 07, 63, 67,
                  99, 65, 04, 28, 06, 16, 70, 92,
                41, 41, 26, 56, 83, 40, 80, 70, 33,
              41, 48, 72, 33, 47, 32, 37, 16, 94, 29,
            53, 71, 44, 65, 25, 43, 91, 52, 97, 51, 14,
          70, 11, 33, 28, 77, 73, 17, 78, 39, 68, 17, 57,
        91, 71, 52, 38, 17, 14, 91, 43, 58, 50, 27, 29, 48,
      63, 66, 04, 68, 89, 53, 67, 30, 73, 16, 69, 87, 40, 31,
    04, 62, 98, 27, 23, 09, 70, 98, 73, 93, 38, 53, 60, 04, 23,
];

在数组中的第n个数,行号是多少?利用了第r行一定有r个数的性质。

fn row(n: usize) -> usize {
    let mut s = 0;
    for r in 1.. {
        s += r;
        if s > n {return r;}
    }
    return 0;
}

根据上面行号的性质,可以得出节点n的下一层的左节点的编号是 n + r,右节点是 n + r + 1。

求路径的和的时候可以利用递归,终止条件是遇到最底一层的时候,由于原题只让求路径长度,这里没有记下来所走路径的编号。

fn max_path_weight(w: &[u32], n: usize) -> u32 {
    if n >= w.len() {return 0;} //越界判断
    let r = row(n);
    let bottom_row = row(w.len() - 1);
    if r == bottom_row { // 递归的退出条件
        return w[n];
    }
    let left = max_path_weight(w, n + r);
    let right = max_path_weight(w, n + r + 1);
    let max = if left > right {left} else {right};
    return w[n] + max;
}

主程序调用只需要一行,数组的总层数不多,复杂度不高,没再做进一步的性能优化。

println!("{}", max_path_weight(&w, 0));

第19题

问题描述:

求20世纪(1901年1月1日到2000年12月31日)有多少个月的一号是星期日?

本题当然可以利用闰年的性质,只用数学公式就能算出来,这里用编程办法,熟悉一下日期时间的函数库。

关于日期的库用chrono,网上有些资料比较老,建议直接参考官网上的帮助,写得非常详细,少走一些弯路。

在https://docs.rs 网站上搜索chrono即可。

use chrono::prelude::*;
use time::Duration;

代码简单粗暴,每次加一天,用到Duration;判断星期几用到weekday(),判断几号用了day(),逻辑很简单。

let mut count = 0;
let mut d = Utc.ymd(1901, 1, 1);
while d <= Utc.ymd(2000, 12, 31) {
    if d.weekday() == Weekday::Sun && d.day() == 1 {
        println!("{}", d);
        count += 1;
    }
    d = d + Duration::days(1);
}
println!("{}", count);

第20题

问题描述:

求100的阶乘中所有数字之和。

本题与第16题非常相似,稍微修改就出来,不解释。

let mut prod = BigUint::from(1 as u64);
for i in 1..=100 {
    prod *= BigUint::from(i as u64);
}

let full_str = prod.to_str_radix(10);
let s = full_str
    .chars()
    .map(|c| c.to_digit(10).unwrap())
    .sum::<u32>();
println!("{}", s);

第21题

问题描述:

10000之内的所有亲和数之和。所谓亲和数,是指两个正整数中,彼此的全部约数之和(本身除外)与另一方相等。比如,220的因子有1, 2, 4, 5, 10, 11, 20, 22, 44, 55 和 110,因子之和是 284,而284的所有因子是 1, 2, 4, 71 和 142,因子之和是220。

问题分解:

  • 求所有因子
  • 因子求和
  • 找出亲和数
  • 亲和数再求和

在第12题里已经求出了一半的因子,函数是:

fn half_factors(num: u32) -> Vec<u32> {
    let s = (num as f32).sqrt() as u32;
    (1..=s).filter(|x| num % x == 0).collect::<Vec<u32>>()
}

很容易补上后面的一半因子。

fn proper_divisors(num: u32) -> Vec<u32> {
    let mut v = half_factors(num);
    for i in (1..v.len()).rev() { //不要num自身,所以从1开始
        v.push(num / v[i]);
    }
    v
}

所有因子求和:

fn proper_divisors_sum(num: u32) -> u32 {
    let divs = proper_divisors(num);
    divs.iter().sum::<u32>()
}

主程序暴力循环即可,10000之内只有5对亲和数:

let mut sum = 0;
for a in 1u32..10000 {
    let b = proper_divisors_sum(a);
    if a != b && proper_divisors_sum(b) == a {
        sum += a;
        println!("{} {}", a, b);
    }
}
println!("{}", sum);

因为亲和数是成对出现的,还可以优化性能,引入一个数组,这里不再展开了。


在projecteuler中注册一个账号,可以添加好友,一起讨论学习,我的Key是:

1539870_KBNiIXymh4SnmDEDZmUTg7tu1MTBVlLj

 类似资料: