; phone numbers and prime factors
 
(define rand
  (let* ((a 6364136223846793005)
         (c 1442695040888963407)
         (m 18446744073709551616)
         (seed (current-time)))
    (lambda args
      (when (pair? args) (set! seed (modulo (car args) m)))
      (set! seed (modulo (+ (* a seed) c) m))
      (/ seed m))))
 
(define randint
  (case-lambda
    ((n) (floor (* (rand) n)))
    ((first past) (+ (floor (* (rand) (- past first))) first))))
 
(define (uniq-c eql? xs)
  (if (null? xs) xs
    (let loop ((xs (cdr xs)) (prev (car xs)) (k 1) (result '()))
      (cond ((null? xs) (reverse (cons (cons prev k) result)))
            ((eql? (car xs) prev) (loop (cdr xs) prev (+ k 1) result))
            (else (loop (cdr xs) (car xs) 1 (cons (cons prev k) result)))))))
 
(define (last xs)
  (if (or (null? xs) (null? (cdr xs))) xs
    (last (cdr xs))))
 
(define (cycle . xs) (set-cdr! (last xs) xs) xs)
 
(define wheel
  (cons 1 (cons 2 (cons 2 (cycle 4 2 4 2 4 6 2 6)))))
 
(define (factors n)
  (let loop ((n n) (f 2) (w wheel) (fs (list)))
    (if (< n (* f f)) (reverse (cons n fs))
      (if (zero? (modulo n f)) (loop (/ n f) f w (cons f fs))
        (loop n (+ f (car w)) (cdr w) fs)))))
 
(define (omega n)
  (let loop ((fs (factors n)) (prev 0) (count 0))
    (if (null? fs) count
      (if (= (car fs) prev)
          (loop (cdr fs) prev count)
          (loop (cdr fs) (car fs) (+ count 1))))))
 
(define (rand-phone) (randint #e1e9 #e1e10))
 
(define (phones n)
  (do ((n n (- n 1))
       (xs (list) (cons (omega (rand-phone)) xs)))
      ((zero? n) (uniq-c = (sort xs <)))))
 
(display (phones 100)) (newline)

OyBwaG9uZSBudW1iZXJzIGFuZCBwcmltZSBmYWN0b3JzCgooZGVmaW5lIHJhbmQKICAobGV0KiAoKGEgNjM2NDEzNjIyMzg0Njc5MzAwNSkKICAgICAgICAgKGMgMTQ0MjY5NTA0MDg4ODk2MzQwNykKICAgICAgICAgKG0gMTg0NDY3NDQwNzM3MDk1NTE2MTYpCiAgICAgICAgIChzZWVkIChjdXJyZW50LXRpbWUpKSkKICAgIChsYW1iZGEgYXJncwogICAgICAod2hlbiAocGFpcj8gYXJncykgKHNldCEgc2VlZCAobW9kdWxvIChjYXIgYXJncykgbSkpKQogICAgICAoc2V0ISBzZWVkIChtb2R1bG8gKCsgKCogYSBzZWVkKSBjKSBtKSkKICAgICAgKC8gc2VlZCBtKSkpKQoKKGRlZmluZSByYW5kaW50CiAgKGNhc2UtbGFtYmRhCiAgICAoKG4pIChmbG9vciAoKiAocmFuZCkgbikpKQogICAgKChmaXJzdCBwYXN0KSAoKyAoZmxvb3IgKCogKHJhbmQpICgtIHBhc3QgZmlyc3QpKSkgZmlyc3QpKSkpCgooZGVmaW5lICh1bmlxLWMgZXFsPyB4cykKICAoaWYgKG51bGw/IHhzKSB4cwogICAgKGxldCBsb29wICgoeHMgKGNkciB4cykpIChwcmV2IChjYXIgeHMpKSAoayAxKSAocmVzdWx0ICcoKSkpCiAgICAgIChjb25kICgobnVsbD8geHMpIChyZXZlcnNlIChjb25zIChjb25zIHByZXYgaykgcmVzdWx0KSkpCiAgICAgICAgICAgICgoZXFsPyAoY2FyIHhzKSBwcmV2KSAobG9vcCAoY2RyIHhzKSBwcmV2ICgrIGsgMSkgcmVzdWx0KSkKICAgICAgICAgICAgKGVsc2UgKGxvb3AgKGNkciB4cykgKGNhciB4cykgMSAoY29ucyAoY29ucyBwcmV2IGspIHJlc3VsdCkpKSkpKSkKCihkZWZpbmUgKGxhc3QgeHMpCiAgKGlmIChvciAobnVsbD8geHMpIChudWxsPyAoY2RyIHhzKSkpIHhzCiAgICAobGFzdCAoY2RyIHhzKSkpKQoKKGRlZmluZSAoY3ljbGUgLiB4cykgKHNldC1jZHIhIChsYXN0IHhzKSB4cykgeHMpCgooZGVmaW5lIHdoZWVsCiAgKGNvbnMgMSAoY29ucyAyIChjb25zIDIgKGN5Y2xlIDQgMiA0IDIgNCA2IDIgNikpKSkpCgooZGVmaW5lIChmYWN0b3JzIG4pCiAgKGxldCBsb29wICgobiBuKSAoZiAyKSAodyB3aGVlbCkgKGZzIChsaXN0KSkpCiAgICAoaWYgKDwgbiAoKiBmIGYpKSAocmV2ZXJzZSAoY29ucyBuIGZzKSkKICAgICAgKGlmICh6ZXJvPyAobW9kdWxvIG4gZikpIChsb29wICgvIG4gZikgZiB3IChjb25zIGYgZnMpKQogICAgICAgIChsb29wIG4gKCsgZiAoY2FyIHcpKSAoY2RyIHcpIGZzKSkpKSkKCihkZWZpbmUgKG9tZWdhIG4pCiAgKGxldCBsb29wICgoZnMgKGZhY3RvcnMgbikpIChwcmV2IDApIChjb3VudCAwKSkKICAgIChpZiAobnVsbD8gZnMpIGNvdW50CiAgICAgIChpZiAoPSAoY2FyIGZzKSBwcmV2KQogICAgICAgICAgKGxvb3AgKGNkciBmcykgcHJldiBjb3VudCkKICAgICAgICAgIChsb29wIChjZHIgZnMpIChjYXIgZnMpICgrIGNvdW50IDEpKSkpKSkKCihkZWZpbmUgKHJhbmQtcGhvbmUpIChyYW5kaW50ICNlMWU5ICNlMWUxMCkpCgooZGVmaW5lIChwaG9uZXMgbikKICAoZG8gKChuIG4gKC0gbiAxKSkKICAgICAgICh4cyAobGlzdCkgKGNvbnMgKG9tZWdhIChyYW5kLXBob25lKSkgeHMpKSkKICAgICAgKCh6ZXJvPyBuKSAodW5pcS1jID0gKHNvcnQgeHMgPCkpKSkpCgooZGlzcGxheSAocGhvbmVzIDEwMCkpIChuZXdsaW5lKQ==