A hobby implementation of linearcongruential generators.

lcg.clj

(defn lcg
  "Returns a lazy seq of random numbers using linear congruential
  generators (LCGs)."
  [seed a b m]
  (map #(float (/ %))
       (drop 1
             (iterate (fn [x] (mod (+ (* a x) b) m))
             seed))))

(defn lcg-average
  "Returns average of x random numbers using LCGs with seed."
  [x seed]
  (let [ans (take x (lcg seed 997 1 65536))]
    (/ (apply + ans) x)))