日志

倒数之和

热度 2已有 1202 次阅读2015-11-1 23:16 |个人分类:LISP|系统分类:应用| 自然数, 倒数

本帖最后由 aimisiyou 于 2015-10-14 23:28 编辑
1/a+1/b+1/c=1/n,a<=b<=c,且都为自然数,已知自然数n求所有的解？
(defun ff (n)
(setq i (+ n 1) lst nil)
(while (<= i (* 3 n))
   (setq n1 (/ (* n i 1.0) (- i n)))
   (setq j (+ (fix n1) 1))
   (while (<= j (* 2 n1))
      (setq n3 (/ (* n i j 1.0) (- (* i j) (* n i) (* n j))))
      (if (= n3 (fix n3))
            (progn
               (setq qlst (mapcar 'fix (vl-sort (mapcar '(lambda (x) (* x 1.0)) (list i j n3)) '<) ))
               (if (= nil (member qlst lst))(setq lst (cons qlst lst)))
            )
         )
      (setq j (+ j 1))
      )
   (setq i (+ i 1))
   )
(princ (length lst))
(princ)
(princ lst)
(princ)
)
(ff 100)