kiyo_hikoのブログ

メモ+日記?

iotaを書いてパスカルの三角形を作った

iota10分ぐらい、三角形10分ぐらいでやっつけた。
大したことはやってないし多分遅い。

(defun iota (to &optional (from 0) (step 1))
  (do ((i from (incf i step))
	   (ns nil))
	  ((< to i) (nreverse ns))
	(setf ns (cons i ns))))

(defun pas (n k)
  (if (or (= n 1) (= k 1) (= k n)) 1
	(+ (pas (- n 1) (- k 1)) (pas (- n 1) k))))

(defun line (n)
  (mapcar #'(lambda (k) (pas n k)) (iota n 1)))

(defun sankaku (n)
  (mapcar #'(lambda (n) (line n)) (iota n 1)))

(defun main()
  (mapc #'(lambda (x) (format t "~a~%" x)) (sankaku 15))
  nil)

出力。

> (main)
(1)
(1 1)
(1 2 1)
(1 3 3 1)
(1 4 6 4 1)
(1 5 10 10 5 1)
(1 6 15 20 15 6 1)
(1 7 21 35 35 21 7 1)
(1 8 28 56 70 56 28 8 1)
(1 9 36 84 126 126 84 36 9 1)
(1 10 45 120 210 252 210 120 45 10 1)
(1 11 55 165 330 462 462 330 165 55 11 1)
(1 12 66 220 495 792 924 792 495 220 66 12 1)
(1 13 78 286 715 1287 1716 1716 1287 715 286 78 13 1)
(1 14 91 364 1001 2002 3003 3432 3003 2002 1001 364 91 14 1)
NIL

追記

iotaの返却値を勘違いしてたっぽい。直した。

http://kiyohikou.hateblo.jp/entry/2013/05/29/175726