平方採中法をDr.Racketで表現する方法はありますか？

平方採中法は、生成する乱数の桁数を n としたとき、初期値 s の2乗を計算し、その数値を 2n 桁の数値とみて、（下の例のように 2 乗した桁数が足りないときは、0 を補います。）その中央にある n 個の数字を最初の乱数とします。次にこの乱数を 2 乗して、同じ様に、中央にある n 個の数字をとって、次の乱数とします。例えば、123 を初期値とすると

1232 = 00015129 → 0151
1512 = 00022801 → 0228
2282 = 00051984 → 0519
5192 = 00269361 → 2693
26932 = 07252249 → 2522
の様になります。この方法を用いて、初期値を入力とし乱数を２個生成し出力するプログラムを作成しようと考えています。そのあと、テストもします。

(define (square x) (* x x))
(define (integer x) (+ x x))

(test-case
"Case 2 - Previous Number: 3456" (displayln "Previous: 3456")
(display "Square of Previous: ")(displayln (expt 3456 2))
(display "Expected: 9439")
(check-equal? (von-Neumann 3456) 9439))

ただ、この後、2乗したものから真ん中の桁だけを抽出する方法がおもいつきません。

dodox86

2021/02/26 23:54

お使いの「Dr.Racket」とはプログラミング言語Schemeの実行環境を言われているのだと思いますが、タグにつけている「rkt」はそれとは全然別モノなので、「Scheme」に変更しましょう。そうすることでSchemeに詳しい方の回答を得られる可能性が高くなります。

cha-

2021/02/27 02:53

ありがとうございます、修正しました。

dodox86

2021/02/27 06:10

私も誤記がありましたので、訂正しておきます。正：DrRacket 誤：Dr.Racket https://download.racket-lang.org/

行動規範の内容に同意します

回答2件

具体的なコードで書くとこんな感じですかね。

scheme
1#!r6rs
2(import (rnrs)
3        (rackunit)
4        (only (racket base) displayln))
5
6;; 整数 num の中央 column 桁を抽出する手続
7(define (middle num column)
8  (assert (integer? num))
9  (assert (integer? column))
10  (let ((lhs (expt 10 (div column 2)))
11        (rhs (expt 10 column)))
12    (mod (div num lhs) rhs)))
13
14;; 連続して乱数生成するジェネレータを生成する手続き
15(define (make-random-generator seed column)
16  (lambda()
17    (set! seed (middle (* seed seed) column))
18    seed))
19
20;; 123 を種とする乱数を 5 個生成して表示
21(do ((count 5 (- count 1))
22     (random-generator (make-random-generator 123 4)))
23    ((zero? count))
24  (display (random-generator))
25  (newline))
26
27;;; テストケースを通す
28(let ((random-generator (make-random-generator 3456 4)))
29  (test-case
30   "Case 2 - Previous Number: 3456"
31   (displayln "Previous: 3456")
32   (display "Square of Previous: ")
33   (displayln (expt 3456 2))
34   (displayln "Expected: 9439")
35   (check-equal? (random-generator) 9439)))