; 357 numbers
(define-syntax while
(syntax-rules ()
((while pred? body ...)
(do () ((not pred?)) body ...))))
(define (range . args)
(case (length args)
((1) (range 0 (car args) (if (negative? (car args)) -1 1)))
((2) (range (car args) (cadr args) (if (< (car args) (cadr args)) 1 -1)))
((3) (let ((le? (if (negative? (caddr args)) >= <=)))
(let loop ((x(car args)) (xs '()))
(if (le? (cadr args) x)
(reverse xs)
(loop (+ x (caddr args)) (cons x xs))))))
(else (error 'range "unrecognized arguments"))))
(define (filter pred? xs)
(let loop ((xs xs) (ys '()))
(cond ((null? xs) (reverse ys))
((pred? (car xs))
(loop (cdr xs) (cons (car xs) ys)))
(else (loop (cdr xs) ys)))))
(define (357? n)
(while (zero? (modulo n 3)) (set! n (/ n 3)))
(while (zero? (modulo n 5)) (set! n (/ n 5)))
(while (zero? (modulo n 7)) (set! n (/ n 7)))
(= n 1))
(define (list357 n) (filter 357? (range 1 n 2)))
(time (display (length (list357 1000000))))
(newline)
(define pq-empty '())
(define pq-empty? null?)
(define (pq-first pq)
(if (null? pq)
(error 'pq-first "can't extract minimum from null queue")
(car pq)))
(define (pq-merge lt? p1 p2)
(cond ((null? p1) p2)
((null? p2) p1)
((lt? (car p2) (car p1))
(cons (car p2) (cons p1 (cdr p2))))
(else (cons (car p1) (cons p2 (cdr p1))))))
(define (pq-insert lt? x pq)
(pq-merge lt? (list x) pq))
(define (pq-merge-pairs lt? ps)
(cond ((null? ps) '())
((null? (cdr ps)) (car ps))
(else (pq-merge lt? (pq-merge lt? (car ps) (cadr ps))
(pq-merge-pairs lt? (cddr ps))))))
(define (pq-rest lt? pq)
(if (null? pq)
(error 'pq-rest "can't delete minimum from null queue")
(pq-merge-pairs lt? (cdr pq))))
(define (list357 n)
(let ((pq (pq-insert < 1 pq-empty)))
(let loop ((357s (list)))
(if (pq-empty? pq) (reverse 357s)
(let ((x (pq-first pq)))
(while (and (not (pq-empty? pq))
(= (pq-first pq) x))
(set! pq (pq-rest < pq)))
(when (< (* 3 x) n)
(set! pq (pq-insert < (* 3 x) pq)))
(when (< (* 5 x) n)
(set! pq (pq-insert < (* 5 x) pq)))
(when (< (* 7 x) n)
(set! pq (pq-insert < (* 7 x) pq)))
(loop (cons x 357s)))))))
(time (display (length (list357 1000000))))