euler / clojure / 21.clj

; Let d(n) be defined as the sum of proper divisors of n (numbers less than n
; which divide evenly into n).
; If d(a) = b and d(b) = a, where a ≠ b, then a and b are an amicable pair and
; each of a and b are called amicable numbers.
; 
; For example, the proper divisors of 220 are 1, 2, 4, 5, 10, 11, 20, 22, 44, 55
; and 110; therefore d(220) = 284. The proper divisors of 284 are 1, 2, 4, 71
; and 142; so d(284) = 220.
; 
; Evaluate the sum of all the amicable numbers under 10000.
;
; Answer: 31626


(defn factors [n]
  (let [max (inc (int (Math/sqrt n)))]
   (loop [fs [] i 2]
     (if (= i max)
       (set (cons 1 fs)) ; Unique items
       (if (zero? (mod n i))
         (recur (concat fs [i (/ n i)]) (inc i))
         (recur fs (inc i)))))))

(defn d [n]
  (apply + (factors n)))

(defn amicable? [n]
  (let [i (d n)]
    (and (not (= n i)) (= (d i) n))))

(prn (apply + (filter amicable? (range 1 10001))))
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.