Gaucheでtraceするまで

Gauche memo

dankogaiさんのエントリを見ると、どうやらtraceすれば分かりやすくなりそう。なので書いてある通りにslibをインストールしようとする

cd ~
curl -O http://swissnet.ai.mit.edu/ftpdir/scm/slib3b3.tar.gz
cd /usr/local
sudo tar zxvf ~/slib3b3.tar.gz

が、

tar: Unrecognized archive format
tar: Error exit delayed from previous errors.

というエラーが解凍時に出てきて進めない。ファイルが壊れてるのかな？


何回か試したけどやっぱり出来ないので、こっちのやり方に変更。

まずココからslibをDLして/usr/localにダウンロードファイル(slib-3b3.zip)を置いておく。
で、

cd /usr/local/
unzip slib-3b3.zip
cd slib
sudo make infoz
sudo make install

とやった後にGaucheをリビルドする。ついでにスレッドを使えるように(参照)しておく。

cd Gauche-0.9.2
 ./configure --enable-threads=pthreads --enable-multibyte=utf-8
 make && make check && sudo make install

これでtraceが使えるはず。


試してみる。

gosh> (use slib)
#<undef>
gosh> (require 'trace)
#t
gosh> (define (fact n) (if (<= n 1) 1 (* n (fact (- n 1)))))
fact
gosh> (fact 10)
3628800
gosh> (trace fact)
#<closure (debug:trace-procedure debug:trace-procedure)>
gosh> (fact 10)
CALL fact 10
  CALL fact 9
    CALL fact 8
      CALL fact 7
	CALL fact 6
	RETN fact 720
      RETN fact 5040
    RETN fact 40320
  RETN fact 362880
RETN fact 3628800
3628800

gosh> (untrace fact)
#<closure fact>
gosh> (fact 10)
3628800

できてますねー。

Gaucheをemacsで書く際のインデント設定->gauche-mode.el

emacs Gauche

※ メモ的エントリです。


let1を使うとインデント崩れるなー。letでいいやー。
というプログラマとしてあるまじき状態だったのですが、ココの内容をそのまま.emacsにコピペして解決しました。どうやら素のScheme-modeでは特殊形式や構文にあんまり対応してないようですね。


let1が使える…これで勝つる…!


(追記)
コメントでgauche-mode.elのことを教えていただいたのでそちらに変更しました。

まず、ココからgauche-mode.elをダウンロードして、パスの通ったところ(M-x describe-variable load-pathで確認可)に入れます。
で、.emacsに

(autoload 'gauche-mode "gauche-mode" nil t)

(setq auto-mode-alist (cons '("\\.scm$" . gauche-mode) auto-mode-alist))

を追加すると.scmファイルを開いたときに自動的にgaucheモードで開かれます。
(追記終)


※ 一応参照先の内容をこちらにもコピーさせていただきます。(更に元はコチラでした)

(put 'and-let* 'scheme-indent-function 1)
(put 'begin0 'scheme-indent-function 0)
(put 'call-with-client-socket 'scheme-indent-function 1)
(put 'call-with-input-conversion 'scheme-indent-function 1)
(put 'call-with-input-file 'scheme-indent-function 1)
(put 'call-with-input-process 'scheme-indent-function 1)
(put 'call-with-input-string 'scheme-indent-function 1)
(put 'call-with-iterator 'scheme-indent-function 1)
(put 'call-with-output-conversion 'scheme-indent-function 1)
(put 'call-with-output-file 'scheme-indent-function 1)
(put 'call-with-output-string 'scheme-indent-function 0)
(put 'call-with-temporary-file 'scheme-indent-function 1)
(put 'call-with-values 'scheme-indent-function 1)
(put 'dolist 'scheme-indent-function 1)
(put 'dotimes 'scheme-indent-function 1)
(put 'if-match 'scheme-indent-function 2)
(put 'let*-values 'scheme-indent-function 1)
(put 'let-args 'scheme-indent-function 2)
(put 'let-keywords* 'scheme-indent-function 2)
(put 'let-match 'scheme-indent-function 2)
(put 'let-optionals* 'scheme-indent-function 2)
(put 'let-syntax 'scheme-indent-function 1)
(put 'let-values 'scheme-indent-function 1)
(put 'let/cc 'scheme-indent-function 1)
(put 'let1 'scheme-indent-function 2)
(put 'letrec-syntax 'scheme-indent-function 1)
(put 'make 'scheme-indent-function 1)
(put 'match 'scheme-indent-function 1)
(put 'match-lambda 'scheme-indent-function 1)
(put 'match-let 'scheme-indent-fucntion 1)
(put 'match-let* 'scheme-indent-fucntion 1)
(put 'match-letrec 'scheme-indent-fucntion 1)
(put 'match-let1 'scheme-indent-function 2)
(put 'match-define 'scheme-indent-fucntion 1)
(put 'multiple-value-bind 'scheme-indent-function 2)
(put 'parameterize 'scheme-indent-function 1)
(put 'parse-options 'scheme-indent-function 1)
(put 'receive 'scheme-indent-function 2)
(put 'rxmatch-case 'scheme-indent-function 1)
(put 'rxmatch-cond 'scheme-indent-function 0)
(put 'rxmatch-if  'scheme-indent-function 2)
(put 'rxmatch-let 'scheme-indent-function 2)
(put 'syntax-rules 'scheme-indent-function 1)
(put 'unless 'scheme-indent-function 1)
(put 'until 'scheme-indent-function 1)
(put 'when 'scheme-indent-function 1)
(put 'while 'scheme-indent-function 1)
(put 'with-builder 'scheme-indent-function 1)
(put 'with-error-handler 'scheme-indent-function 0)
(put 'with-error-to-port 'scheme-indent-function 1)
(put 'with-input-conversion 'scheme-indent-function 1)
(put 'with-input-from-port 'scheme-indent-function 1)
(put 'with-input-from-process 'scheme-indent-function 1)
(put 'with-input-from-string 'scheme-indent-function 1)
(put 'with-iterator 'scheme-indent-function 1)
(put 'with-module 'scheme-indent-function 1)
(put 'with-output-conversion 'scheme-indent-function 1)
(put 'with-output-to-port 'scheme-indent-function 1)
(put 'with-output-to-process 'scheme-indent-function 1)
(put 'with-output-to-string 'scheme-indent-function 1)
(put 'with-port-locking 'scheme-indent-function 1)
(put 'with-string-io 'scheme-indent-function 1)
(put 'with-time-counter 'scheme-indent-function 1)
(put 'with-signal-handlers 'scheme-indent-function 1)

Project EulerをSchemeで(16-20)

Gauche Project Euler

今回は問16から問20の5問。処理系はGauche。
15問まとめて書いたのはちょっと長すぎたと思います。


まずはユーティリティ関数とか。前回でも使いました。

(use srfi-1)

; (integer->list 12345) => (1 2 3 4 5)
(define (integer->list i)
  (letrec ((i->rl (lambda (i)
		    (cons (modulo i 10)
			  (if (< i 10)
			      '()
			      (i->rl (quotient i 10)))))))
    (reverse (i->rl i))))

問16.

2の1000乗の各桁を足した数を求める問題。
2の1000乗が普通に計算できることにびっくり。

(define (e16)
  (apply + (integer->list (expt 2 1000))))

問17.

1から1000までを英語で書いたときの文字数の和を求める問題。
ベクタに1から1000までを英語で書きだして、最後に文字数の和を求める。

(define (e17)
  (let ((v (make-vector 1001 "")))
    (for-each (lambda (i s) (vector-set! v i s))
	      (iota 19 1) '("one" "two" "three" "four" "five" "six" "seven" "eight" "nine" "ten" "eleven" "twelve" "thirteen" "fourteen" "fifteen" "sixteen" "seventeen" "eighteen" "nineteen"))
    (for-each (lambda (i s) 
		(for-each (lambda (j) 
			    (vector-set! v (+ i j) (string-append s (vector-ref v j)))) 
			  (iota 10)))
	      (iota 9 20 10) '("twenty" "thirty" "forty" "fifty" "sixty" "seventy" "eighty" "ninety"))
    (for-each (lambda (i) 
		(for-each (lambda (j)
			    (vector-set! v
					 (+ (* i 100) j) 
					 (string-append (vector-ref v i) "hundred"
							(if (not (zero? j))
							    (string-append "and" (vector-ref v j))
							    ""))))
			  (iota 100)))
	      (iota 9 1))
    (vector-set! v 1000 "onethousand")
    (apply + (map (lambda (i) (string-length (vector-ref v i))) (iota 1000 1)))))

問18.

三角形に並んだ数字の列があって、三角形の頂点から底まで移動するときに通過した数値の和の最大値を求める問題。
最初は泥臭い実装をしていたけれど、人の解答を見てショックを受けた。
fold-rightってこういう使い方するんだ、とか(map max k (cdr k))ってめっちゃシンプルやん、とか。よって丸パクリ。

(define (e18)
  (fold-right (lambda (x m)
		(let ((k (map + x m)))
		  (if (null? (cdr k))
		      (car k)
		      (map max k (cdr k)))))
	      (iota 15 0 0)
	      '((75)
		(95 64)
		(17 47 82)
		(18 35 87 10)
		(20 04 82 47 65)
		(19 01 23 75 03 34)
		(88 02 77 73 07 63 67)
		(99 65 04 28 06 16 70 92)
		(41 41 26 56 83 40 80 70 33)
		(41 48 72 33 47 32 37 16 94 29)
		(53 71 44 65 25 43 91 52 97 51 14)
		(70 11 33 28 77 73 17 78 39 68 17 57)
		(91 71 52 38 17 14 91 43 58 50 27 29 48)
		(63 66 04 68 89 53 67 30 73 16 69 87 40 31)
		(04 62 98 27 23 09 70 98 73 93 38 53 60 04 23))))

問19.

20世紀中に月の初め(1日)が日曜日だったのは何回あったかという問題。
1世紀分の月の日数リストを作って各月が日曜日で始まっているか判定している。

(define (e19)
  (let ((sd (modulo (+ 1 365) 7)) 
	(ml (let mloop ((l '()) (yl (iota 100 1)))
	      (if (null? yl)
		  l
		  (mloop (append l (list 31 (if (zero? (modulo (car yl) 4)) 29 28) 31 30 31 30 31 31 30 31 30 31))
			 (cdr yl))))))
    (let loop ((d sd) (res (if (zero? sd) 1 0)) (ml ml))
      (if (null? ml)
	  res
	  (loop (modulo (+ d (car ml)) 7)
		(if (zero? d) (+ 1 res) res)
		(cdr ml))))))

問20.

100の階乗の各桁の和を求める問題。
まぁそのまんま。

(define (e20)
  (apply + (integer->list (apply * (iota 100 1)))))

SchemeでProject Euler(1-15 )

Gauche Project Euler

Rubyを勉強しようと思っていたらいつの間にかSchemeを勉強していた(ry

というわけでProject Eulerに挑戦しています。


といってもSchemeはやったことないので、入門としてもうひとつの Scheme 入門にお世話になりました。
また、Scheme的な書き方も分からないので、分からなくなったらとりあえず人の解答を見て写経しています。
いくつかの問題は自分一人で解ける程度には慣れてきた…はず。

(use srfi-1)

; (sum '(1 2 3)) => 6
(define (sum l)
  (fold + 0 l))

(define (e1)
  (sum (filter
	(lambda (x) (or (zero? (modulo x 5))
			(zero? (modulo x 3))))
	(iota 1000))))

(use util.stream)

(define (e2)
  (letrec ((fib (lambda (a b)
		  (stream-delay (stream-cons a (fib b (+ a b)))))))
    (sum (stream->list (stream-filter even? (stream-take-while (cut < <> (* 400 10000)) (fib 1 2)))))))

(追記)
コメントで再帰でも出来る(し早い)と教えていただいたので再帰版も。(むしろ普通はこちらですかね…)

(define (e2-recursive)
  (let loop ((a 1) (b 2) (res 2))
    (let1 c (+ a b)
      (if (>= c 4000000) ;(* 400 10000)
	  res
	  (loop b c (if (even? c) (+ res c) res))))))

(さらに追記)
上のは再帰じゃなく反復ですよね…(自信なし)
ので改めて再帰版、末尾再帰版、反復版を。

(define (e2-rec)
  (define (rec a b)
    (let1 c (+ a b)
      (if (>= c 4000000)
	0
        (+ (if (even? c) c 0)
           (rec	b c)))))
  (rec 1 1))

(define (e2-tail)
  (define (rec a b r)
    (let1 c (+ a b)
      (if (>= c 4000000)
          r
          (rec b c (if (even? c) (+ r c) r)))))
  (rec 1 2 2))

(define (e2-iter)
  (let loop ((a 1) (b 2) (r 2))
    (let1 c (+ a b)
      (if (>= c 4000000)
          r
          (loop b c (if (even? c) (+ r c) r))))))

;; Sieve of Eratosthenes
; (primes 10) => (2 3 5 7)
(define (primes n)
  (if (<= n 2)
      '()
      (let loop ((l (unfold (cut < n <>) values (cut + 2 <>) 3))
		 (prime-list '(2)))
	(let ((m (car l)))
	  (if (> (expt m 2) n)
	      (append (reverse prime-list) l)
	      (loop (remove (lambda (x) (zero? (modulo x m))) l) (cons m prime-list)))))))

(define (e3)
  (let* ((v 600851475143)
	 (l (reverse (primes (truncate (sqrt v))))))
    (if (null? l)
	v
	(let loop ((l l))
	  (if (zero? (modulo v (car l)))
	      (car l)
	      (loop (cdr l)))))))

(define (e4)
  (let ((max2 (lambda (l) (apply max l))))
    (max2 (map (lambda (x)
		 (max2 (map (lambda (y)
			      (let* ((v (* x y)) (l (integer->list v)))
				(if (every = l (reverse l))
 				    v
				    0)))
			    (iota x 1))))
	       (iota 999 1)))))

(define (e5)
  (lcm 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20))

(define (e6)
  (let ((l (iota 100 1)))
    (- (expt (sum l) 2)
       (sum (map (cut expt <> 2) l)))))

(define (e7)
  (list-ref (primes 200000) 10000))

; (integer->list 12345) => (1 2 3 4 5)
(define (integer->list i)
  (letrec ((i->rl (lambda (i)
		    (cons (modulo i 10)
			  (if (< i 10)
			      '()
			      (i->rl (quotient i 10)))))))
    (reverse (i->rl i))))

(define (e8)
  (let ((l (integer->list 7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450)))
    (apply max (map * l (cdr l) (cddr l) (cdddr l) (cddddr l)))))

(define (e9)
  (any (lambda (c)
	 (any (lambda (b)
		(any (lambda (a)
		       (and (= (+ a b c) 1000)
			    (= (+ (* a a) (* b b)) (* c c))
			    (* a b c)))
		     (iota (- b 1) 1)))
	      (iota (- c 1) 1)))
       (iota 500 1)))

(define (e10)
  (sum (primes (* 200 10000))))

; (mul '(1 2 3 4 5)) => 120
(define (mul l)
  (fold * 1 l))

(define (e11)
  (let ((n 20)
	(p '#(08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08
	      49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00
	      81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65
	      52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91
	      22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80
	      24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50
	      32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70
	      67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21
	      24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72
	      21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95
	      78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92
	      16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57
	      86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58
	      19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40
	      04 52 08 83 97 35 99 16 07 97 57 32 16 26 26 79 33 27 98 66
	      88 36 68 87 57 62 20 72 03 46 33 67 46 55 12 32 63 93 53 69
	      04 42 16 73 38 25 39 11 24 94 72 18 08 46 29 32 40 62 76 36
	      20 69 36 41 72 30 23 88 34 62 99 69 82 67 59 85 74 04 36 16
	      20 73 35 29 78 31 90 01 74 31 49 71 48 86 81 16 23 57 05 54
	      01 70 54 71 83 51 54 69 16 92 33 48 61 43 52 01 89 19 67 48)))
    (apply max (append-map (lambda (dl xl yl)
			     (append-map (lambda (y) 
					   (map (lambda (x) 
						  (mul (map (lambda (d) 
							      (vector-ref p (+ d x y)))
							    dl))) 
						xl)) 
					 yl))
			   (list (iota 4) (iota 4 0 n) (iota 4 0 (+ n 1)) (iota 4 0 (- n 1)))
			   (list (iota (- n 3)) (iota n) (iota (- n 3)) (iota (- n 3) 3))
			   (list (iota n 0 n) (iota (- n 3) 0 n) (iota (- n 3) 0 n) (iota (- n 3) 0 n))))))

;(let ((res (factorize 10 (primes 10)))) 
;  (list (hash-table-keys res) (hash-table-values res)))
;=> ((7 5 3 2) (0 1 0 1))
(define (factorize n ps)
  (let ((res (make-hash-table)))
    (for-each (lambda (p) (hash-table-put! res p 0)) ps)
    (let loop ((n n) (ps ps))
      (cond ((null? ps) res)
	    ((zero? (modulo n (car ps)))
	     (begin (hash-table-put! res (car ps) (+ 1 (hash-table-get res (car ps))))
		    (loop (/ n (car ps)) ps)))
	    (else (loop n (cdr ps)))))))

(define (e12)
  (let ((ps (primes 50))
	(tri (lambda (n) (/ (* n (+ n 1)) 2)))
	(factor-counts (lambda (n ps) (hash-table-values (factorize n ps)))))
    (let loop ((n 2))
      (if (<= 501 (mul (map (cut + <> 1) (factor-counts (tri n) ps))))
	  (tri n)
	  (loop (+ n 1))))))

(define (e13)
  (take (integer->list 
	 (sum '(37107287533902102798797998220837590246510135740250
		46376937677490009712648124896970078050417018260538
		74324986199524741059474233309513058123726617309629
		91942213363574161572522430563301811072406154908250
		23067588207539346171171980310421047513778063246676
		89261670696623633820136378418383684178734361726757
		28112879812849979408065481931592621691275889832738
		44274228917432520321923589422876796487670272189318
		47451445736001306439091167216856844588711603153276
		70386486105843025439939619828917593665686757934951
		62176457141856560629502157223196586755079324193331
		64906352462741904929101432445813822663347944758178
		92575867718337217661963751590579239728245598838407
		58203565325359399008402633568948830189458628227828
		80181199384826282014278194139940567587151170094390
		35398664372827112653829987240784473053190104293586
		86515506006295864861532075273371959191420517255829
		71693888707715466499115593487603532921714970056938
		54370070576826684624621495650076471787294438377604
		53282654108756828443191190634694037855217779295145
		36123272525000296071075082563815656710885258350721
		45876576172410976447339110607218265236877223636045
		17423706905851860660448207621209813287860733969412
		81142660418086830619328460811191061556940512689692
		51934325451728388641918047049293215058642563049483
		62467221648435076201727918039944693004732956340691
		15732444386908125794514089057706229429197107928209
		55037687525678773091862540744969844508330393682126
		18336384825330154686196124348767681297534375946515
		80386287592878490201521685554828717201219257766954
		78182833757993103614740356856449095527097864797581
		16726320100436897842553539920931837441497806860984
		48403098129077791799088218795327364475675590848030
		87086987551392711854517078544161852424320693150332
		59959406895756536782107074926966537676326235447210
		69793950679652694742597709739166693763042633987085
		41052684708299085211399427365734116182760315001271
		65378607361501080857009149939512557028198746004375
		35829035317434717326932123578154982629742552737307
		94953759765105305946966067683156574377167401875275
		88902802571733229619176668713819931811048770190271
		25267680276078003013678680992525463401061632866526
		36270218540497705585629946580636237993140746255962
		24074486908231174977792365466257246923322810917141
		91430288197103288597806669760892938638285025333403
		34413065578016127815921815005561868836468420090470
		23053081172816430487623791969842487255036638784583
		11487696932154902810424020138335124462181441773470
		63783299490636259666498587618221225225512486764533
		67720186971698544312419572409913959008952310058822
		95548255300263520781532296796249481641953868218774
		76085327132285723110424803456124867697064507995236
		37774242535411291684276865538926205024910326572967
		23701913275725675285653248258265463092207058596522
		29798860272258331913126375147341994889534765745501
		18495701454879288984856827726077713721403798879715
		38298203783031473527721580348144513491373226651381
		34829543829199918180278916522431027392251122869539
		40957953066405232632538044100059654939159879593635
		29746152185502371307642255121183693803580388584903
		41698116222072977186158236678424689157993532961922
		62467957194401269043877107275048102390895523597457
		23189706772547915061505504953922979530901129967519
		86188088225875314529584099251203829009407770775672
		11306739708304724483816533873502340845647058077308
		82959174767140363198008187129011875491310547126581
		97623331044818386269515456334926366572897563400500
		42846280183517070527831839425882145521227251250327
		55121603546981200581762165212827652751691296897789
		32238195734329339946437501907836945765883352399886
		75506164965184775180738168837861091527357929701337
		62177842752192623401942399639168044983993173312731
		32924185707147349566916674687634660915035914677504
		99518671430235219628894890102423325116913619626622
		73267460800591547471830798392868535206946944540724
		76841822524674417161514036427982273348055556214818
		97142617910342598647204516893989422179826088076852
		87783646182799346313767754307809363333018982642090
		10848802521674670883215120185883543223812876952786
		71329612474782464538636993009049310363619763878039
		62184073572399794223406235393808339651327408011116
		66627891981488087797941876876144230030984490851411
		60661826293682836764744779239180335110989069790714
		85786944089552990653640447425576083659976645795096
		66024396409905389607120198219976047599490197230297
		64913982680032973156037120041377903785566085089252
		16730939319872750275468906903707539413042652315011
		94809377245048795150954100921645863754710598436791
		78639167021187492431995700641917969777599028300699
		15368713711936614952811305876380278410754449733078
		40789923115535562561142322423255033685442488917353
		44889911501440648020369068063960672322193204149535
		41503128880339536053299340368006977710650566631954
		81234880673210146739058568557934581403627822703280
		82616570773948327592232845941706525094512325230608
		22918802058777319719839450180888072429661980811197
		77158542502016545090413245809786882778948721859617
		72107838435069186155435662884062257473692284509516
		20849603980134001723930671666823555245252804609722
		53503534226472524250874054075591789781264330331690))) 10))

(use gauche.sequence)

(define (e14)
  (let ((ht (make-hash-table))
	(range (* 100 10000)))
    (hash-table-put! ht 1 1)
    (for-each (lambda (n)
		(let loop ((n n) (rec '()))
		  (if (hash-table-exists? ht n)
		      (let ((c (+ 1 (hash-table-get ht n))))
			(for-each-with-index (lambda (i k) 
					       (hash-table-put! ht k (+ c i)))
					     rec))
		      (loop (if (even? n) 
				(/ n 2) 
				(+ (* 3 n) 1))
			    (cons n rec)))))
	      (iota range 1))
    (hash-table-fold ht (lambda (k v p)
			  (if (and (<= k range)
				   (> v (hash-table-get ht p)))
			      k p)) 1)))

(define (e15)
  (/ (mul (iota 20 21))
     (mul (iota 20 1))))

「全裸で学ぶMVC事始め」をruby + Sinatraでやってみた

勉強のため、ゆーすけべーさんの「全裸で学ぶMVC事始め」をruby + Sinatraで書いてみました。
MVCとは何か、どのような処理をするのかについては元エントリで分かりやすく解説されているのでそちらを参照してください。(なぜ全裸なのかも)

つくるもの = zenrize

入力されたテキストを改変して表示する簡単なWebサービスを作成します。
例えば「これから歯磨きをするよ」→「これから全裸で歯磨きをするよ」のようになります。
【全裸で】を差し込む箇所については、全裸botの作者sugyanさんが色々調整を行なっていますが、
今回は単純に「文章内の動詞の前に【全裸で】を必ずつける」とします。

なお、テキスト中の品詞を特定するのにはYahoo!形態素解析APIを使用します。

実装

まずはモデル部分です。(こちらのページを大いに大いに*1参考にさせていただきました)
Zenraクラスを定義し、その中にzenrizeメソッドを定義します

class Zenra

  def initialize(position='動詞', text='全裸で')
    @appid    = 'YOUR_APP_ID'
    @position = position
    @text     = text
    @base_url = 'http://jlp.yahooapis.jp/MAService/V1/parse'
  end

  def zenrize(sentence)
    raise "appid is necessary!" if @appid.nil? || @appid.empty?
    raise "Japanese sentence is necessary!" if sentence.nil? || sentence.empty?

    req = "#{@base_url}?results=ma&appid=#{@appid}&sentence=#{URI.encode(sentence)}"
    res = open(req)
    xml = REXML::Document.new(res.read)

    return sentence if xml.elements['ResultSet/ma_result/total_count'].text == '0'

    result = ''
    xml.elements.each('ResultSet/ma_result/word_list/word') do |e|
      result << @text if e.elements['pos'].text == @position
      result << e.elements['surface'].text
    end
    result
  end

end

次にコントローラ部分です。
get, postメソッドに応じたルーティングを設定してリクエストを処理し、ビューを用いてHTMLをレンダリングしています。

get '/' do
  erb :index
end

post '/result' do
  @model = Zenra.new
  erb :result, :locals => {:res => @model.zenrize(params[:text])}
end

最後にビュー部分です。
今回はerbテンプレートを使用しています。

@@ layout
<!DOCTYPE html>
<html>
<head><title>zenrize</title></head>
<body style="width:500px;margin:0 auto;"><h1>zenrize</h1><%= yield %></body>
</html>

@@ index
<form action="/result" method="post">
<textarea rows="3" cols="60" name="text"></textarea>
<br />
<input type='submit' value='zenrize' />
</form>

@@ result
<p><%= h res %></p>
<a href='/'>戻る</a>

以上の3つを繋ぎ合わせたコードは以下のようになります。

# coding: utf-8
require 'rubygems'
require 'sinatra'
require 'open-uri'
require 'rexml/document'

helpers do
  include Rack::Utils; alias_method :h, :escape_html
end

class Zenra

  def initialize(position='動詞', text='全裸で')
    @appid    = 'YOUR_APP_ID'
    @position = position
    @text     = text
    @base_url = 'http://jlp.yahooapis.jp/MAService/V1/parse'
  end

  def zenrize(sentence)
    raise "appid is necessary!" if @appid.nil? || @appid.empty?
    raise "Japanese sentence is necessary!" if sentence.nil? || sentence.empty?

    req = "#{@base_url}?results=ma&appid=#{@appid}&sentence=#{URI.encode(sentence)}"
    res = open(req)
    xml = REXML::Document.new(res.read)

    return sentence if xml.elements['ResultSet/ma_result/total_count'].text == '0'

    result = ''
    xml.elements.each('ResultSet/ma_result/word_list/word') do |e|
      result << @text if e.elements['pos'].text == @position
      result << e.elements['surface'].text
    end
    result
  end

end

get '/' do
  erb :index
end

post '/result' do
  @model = Zenra.new
  erb :result, :locals => {:res => @model.zenrize(params[:text])}
end

__END__

@@ layout
<!DOCTYPE html>
<html>
<head><title>zenrize</title></head>
<body style="width:500px;margin:0 auto;"><h1>zenrize</h1><%= yield %></body>
</html>

@@ index
<form action="/result" method="post">
<textarea rows="3" cols="60" name="text"></textarea>
<br />
<input type='submit' value='zenrize' />
</form>

@@ result
<p><%= h res %></p>
<a href='/'>戻る</a>


これをzenra.rbとして保存し、

ruby zenra.rb

と実行するとWebブラウザからzenrizeを試すことができます

こちらが入力画面

こちらが結果画面

すごく…全裸です…

まとめ的な

"MVCという概念を学ぶにはこうした一連の流れを網羅する小さなサンプルをつくるのが一番取っ付きやすい(元エントリより)"気がします。
Sinatraいいですね。初めて触った(rubyも初心者状態)にも関わらず、すんなり理解出来ました。
全裸寒い。

最後に

上でも述べていますが、ruby初心者+今回がSinatraデビューな未熟者ですので
ここがおかしいとか、こうした方がいい等あればコメント、もしくはTwitterで教えていただけると幸いです。

*1:丸パクリなレベルで。すいません…

さくらのVPSを借りたらとりあえずやっておくべき設定はこれでいいのかしら

memo

「さくらのVPSを借りたらとりあえずやっておくべき設定」みたいなタイトルのブログをいくつか読んだけど
最低限sshは鍵認証にしてiptablesをちょっといじっておけ、という感じみたい。
ここにも一応手順を書いておく(メモ程度だけど)。

rootのパスワード変更
passwd
sshのポート変更

sshデーモンの設定を開いて

vi /etc/ssh/sshd_config

とりあえず22番ポートから変更しておく

Port 10022 # 適当な数字

終わったらsshdの再起動

service sshd restart
作業ユーザの追加

ユーザを作成してパスワードを設定する

useradd nomnel
passwd nomnel

作成したユーザがsudo出来るように、まずwheelグループに追加する

usermod -G wheel nomnel
visudo

visudoで開いたファイルから下の行を探してコメントアウトを外す(%wheelの行)

## Allows people in group wheel to run all commands
%wheel  ALL=(ALL)    ALL
RSA鍵認証でのSSH接続

作業用端末でRSA鍵を作成

ssh-keygen

RSA公開鍵を作業用端末からVPSに転送

scp -P 10022 ~/.ssh/id_rsa.pub nomnel@[IPアドレス]:~/

接続 (以降VPSで)

ssh -p 10022 nomnel@[IPアドレス]

.sshフォルダを準備

cd
mkdir .ssh

公開鍵の内容をauthorized_keysにコピー

cat id_rsa.pub >> .ssh/authorized_keys

パーミッションを変更する

chmod 700 .ssh
chmod 600 .ssh/authorized_keys

sshdの設定を変更する

sudo vi /etc/ssh/sshd_config

下の行を探し、それぞれ"no"に変更する

PermitRootLogin no

PasswordAuthentication no

PermitEmptyPasswords no

sshdを再起動

sudo service sshd restart

iptablesの設定

sudo vi /etc/sysconfig/iptables

ssh, httpのみ開けておく

*filter
:INPUT   ACCEPT [0:0]
:FORWARD ACCEPT [0:0]
:OUTPUT  ACCEPT [0:0]
:RH-Firewall-1-INPUT - [0:0]

-A INPUT -j RH-Firewall-1-INPUT
-A FORWARD -j RH-Firewall-1-INPUT
-A RH-Firewall-1-INPUT -i lo -j ACCEPT
-A RH-Firewall-1-INPUT -p icmp --icmp-type any -j ACCEPT
-A RH-Firewall-1-INPUT -p 50 -j ACCEPT
-A RH-Firewall-1-INPUT -p 51 -j ACCEPT
-A RH-Firewall-1-INPUT -p udp --dport 5353 -d 224.0.0.251 -j ACCEPT
-A RH-Firewall-1-INPUT -p udp -m udp --dport 631 -j ACCEPT
-A RH-Firewall-1-INPUT -p tcp -m tcp --dport 631 -j ACCEPT
-A RH-Firewall-1-INPUT -m state --state ESTABLISHED,RELATED -j ACCEPT

# SSH, HTTP
-A RH-Firewall-1-INPUT -m state --state NEW -m tcp -p tcp --dport 10022  -j ACCEPT
-A RH-Firewall-1-INPUT -m state --state NEW -m tcp -p tcp --dport 80        -j ACCEPT

-A RH-Firewall-1-INPUT -j REJECT --reject-with icmp-host-prohibited

COMMIT

iptablesの再起動とデフォルトで起動するように

sudo /etc/rc.d/init.d/iptables restart
sudo /sbin/chkconfig iptables on
yumのアップデート
sudo yum -y update


とりあえずはこんな感じでいいのだろうか。
後はCentOSをサーバーとして活用するための基本的な設定 - さくらインターネット創業日記で紹介されていたように不要なサービスを落としておこうかな。

(3/12追記)

内部エンコードをUTF-8に変えておく
sudo vi /etc/sysconfig/i18n

で開いて、LANGのところを

LANG="ja_JP.UTF-8"

に書き換える。(設定の反映方法が分かんなかったのでとりあえず再起動…XD)

Gaucheの環境構築メモ

Gauche emacs

SICPが読みたくなってもうひとつの Scheme 入門でSchemeをちょっと勉強してたんですが、「もうすぐ無職の穀潰しがそんな余裕ぶっこいてていいの…？」とか弟に言われたのでちょっと棚上げ中。


一応、環境構築のメモだけ残しておく。
ちなみにOS X Lion。

Gaucheのインストール

ココからDLする。
後は、

$ gzcat Gauche-0.9.2.tgz | tar xf -
$ cd Gauche-0.9.2
$ ./configure --enable-threads=pthreads --enable-multibyte=utf-8  # スレッドを使えるように
$ make
$ make install

参考：本家サイト

追記(2012/03/09)
homebrewでもいけるみたい

$ brew install gauche

参考：MQL5リファレンス - MetaTrader 5 のためのアルゴリズムの/自動化されたトレーディング言語の使い方

Emacs23のインストール

$ curl -O http://ftp.gnu.org/pub/gnu/emacs/emacs-23.4.tar.gz
$ tar xvfz emacs-23.4.tar.gz
$ cd emacs-23.4
$ ./configure --with-ns --without-x
$ make
$ make install
$ open nextstep/Emacs.app

make installすればEmacs.appが生成されるのでApplicationsフォルダに(移動しても問題無い)。


※ GUIが必要なければ

brew install emacs -cocoa

して、.bashrcなどに

#brewでインストールしたemacsに変更
alias emacs='/usr/local/Cellar/emacs/23.4/bin/emacs'

でおｋ。

参考：Emacs23 (Cocoa Emacs) 入門から中毒まで : 紹介マニア

おまじない

.emacs(~/emacs.d/init.el) に以下を記述

(setq scheme-program-name "gosh")
(require 'cmuscheme)

(defun scheme-other-window ()
  "Run scheme on other window"
  (interactive)
  (switch-to-buffer-other-window
   (get-buffer-create "*scheme*"))
  (run-scheme scheme-program-name))

(define-key global-map
  "\C-cS" 'scheme-other-window)

覚えておくこと。

C-c S：Emacsの画面が分割されて、goshのインタープリタが起動
C-c C-l：インタープリタにファイルをロード
C-x C-e：式を評価(ex. (+ 2 3)の直後にカーソルを合わせてC-x C-eするとインタープリタに5が表示)
C-x o：次のウィンドウへ

参考：http://khead.kahua.org/kahua-web/show/eg/emacs


さてこっそりやるか