题目37：找出所有11个可以双向裁剪的质数的和

mathtimes

运行结果

1. $ time target/release/main
   
2. 23
   
3. 53
   
4. 73
   
5. 37
   
6. 313
   
7. 373
   
8. 317
   
9. 797
   
10. 3137
    
11. 3797
    
12. 739397
    
13. 
    
14. real        0m0.007s
    
15. user        0m0.007s
    
16. sys        0m0.000s

复制代码

题目结果

1. sum=0;for i in `target/release/main`; do((sum+=i));done;echo $sum
   
2. 748317

复制代码

使用 primal 库，Miller-Rabin方法检测质数

1. use std::collections::HashSet;
   
2. use primal::*;
   
3. fn prime_r(x:u64) -> bool {
   
4.     if x == 0 { return true;}
   
5.     if ! is_prime(x) {return false;}
   
6.     prime_r(x / 10)
   
7. }
   
8. fn main() {
   
9.     let mut v = vec![2,3,5,7];
   
10.     let mut bit = 1;
    
11.     while ! v.is_empty() {
    
12.         let mut r = Vec::new();
    
13.         for i in v {
    
14.             for j in 1..10 {
    
15.                 let x = i + j * 10u64.pow(bit);
    
16.                 if is_prime(x) {
    
17.                     r.push(x);
    
18.                 }
    
19.             }
    
20.         }
    
21.         v = r;
    
22.         bit += 1;
    
23.         for i in &v {
    
24.             if *i > 1e10 as u64 {
    
25.                 return ();
    
26.             }
    
27.             if prime_r(*i) {
    
28.                 println!("{i}");
    
29.             }
    
30.         }
    
31.     }
    
32. }

复制代码