Problem 10 » 履歴 » バージョン 6
Noppi, 2023/12/27 13:42
1 | 1 | Noppi | [ホーム](https://redmine.noppi.jp) - [[Wiki|Project Euler]] |
---|---|---|---|
2 | # [[Problem 10]] |
||
3 | |||
4 | 6 | Noppi | ## Summation of Primes |
5 | The sum of the primes below $10$ is $2 + 3 + 5 + 7 = 17$. |
||
6 | Find the sum of all the primes below two million. |
||
7 | |||
8 | ## 素数の和 |
||
9 | 10以下の素数の和は 2 + 3 + 5 + 7 = 17 である. |
||
10 | 200万以下の全ての素数の和を求めよ. |
||
11 | |||
12 | 1 | Noppi | ```scheme |
13 | #!r6rs |
||
14 | #!chezscheme |
||
15 | 4 | Noppi | |
16 | (import (chezscheme)) |
||
17 | 2 | Noppi | |
18 | 6 | Noppi | (define (make-prime-list num) |
19 | (fold-left |
||
20 | (lambda (lis iota-num prime?) |
||
21 | (if prime? |
||
22 | (cons iota-num lis) |
||
23 | lis)) |
||
24 | '() |
||
25 | (iota num) |
||
26 | (vector->list (eratosthenes num)))) |
||
27 | 3 | Noppi | |
28 | 4 | Noppi | (define (eratosthenes-sub! table erase-num) |
29 | 3 | Noppi | (assert (vector-ref table erase-num)) |
30 | (let loop ([current (* erase-num 2)]) |
||
31 | (cond |
||
32 | [(<= (vector-length table) current) table] |
||
33 | [else |
||
34 | (vector-set! table current #f) |
||
35 | 4 | Noppi | (loop (+ current erase-num))]))) |
36 | |||
37 | 1 | Noppi | (define (eratosthenes num) |
38 | (let ([table (make-vector num #t)] |
||
39 | [count (isqrt (sub1 num))]) |
||
40 | 4 | Noppi | (vector-set! table 0 #f) |
41 | 3 | Noppi | (vector-set! table 1 #f) |
42 | 1 | Noppi | (let loop ([index 2] [table table]) |
43 | 3 | Noppi | (cond |
44 | 1 | Noppi | [(< count index) table] |
45 | 3 | Noppi | [(vector-ref table index) |
46 | (loop (add1 index) (eratosthenes-sub! table index))] |
||
47 | [else (loop (add1 index) table)])))) |
||
48 | |||
49 | 1 | Noppi | (define answer-10 |
50 | 6 | Noppi | (apply + (make-prime-list 2000001))) |
51 | 1 | Noppi | |
52 | (printf "10: ~D~%" answer-10) |
||
53 | ``` |