about summary refs log tree commit diff
path: root/THT/C/TP-2018/func.scm
blob: a7d0af520646425d4be6b5ce3a60eee8e6ad512e (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
(define (func k r p q m)
  (define (func-iter a b c d count)
    (display (cons a b))
    (newline)
    (cond ((= count 0) (cons (remainder b m)
                             (remainder a m)))
          ((even? count) (func-iter (remainder a m)
                                    (remainder b m)
                                    (remainder (+ (* c c) (* d d)) m)
                                    (remainder (+ (* 2 c d) (* d d r)) m)
                                    (/ count 2)))
          (else (func-iter (remainder (+ (* b d) (* a d r) (* a c)) m)
                           (remainder (+ (* b c) (* a d)) m)
                           (remainder c m)
                           (remainder d m)
                           (- count 1)))))
  (let ((g (gcd p q)))
    (func-iter (/ q g) (/ p g) 0 1 (- k 1))))

(define (iter)
  (let* ((k (read))
         (r (read))
         (p (read))
         (q (read))
         (m (read)))
    (if (not (eof-object? m))
        (let ((c (func k r p q m)))
          (display (car c))
          (display " ")
          (display (cdr c))
          (newline)
          (iter)))))

(with-input-from-file "FUNC.INP"
  (lambda () (with-output-to-file "FUNC.OUT" iter)))