17122007 life
昨日、富士のオートバックス(スーパーオートバックスっていうの)にいってスタッドレスに履き替えようとして、
「タイヤの付け替えしたいんですけど」
「はぁ～、付け替えってなんですか?」
って受付に言われた上、料金が倍、さらに3時間待ちって言われたので、やめて帰ってきた。そのときはスーパーだけに工賃高いのかなーって思ったけど、単に間違えてただけだろう。
あと、付け替えってオートバックス用語じゃろが～ってちょっとむっとしたけど、富士は雪なんて降んないからスタッドレス履かないからわかんないんじゃないの?って富士の人に言われてちょっと納得。
というわけで、タイヤを詰んだまま会社。帰りにいつも使ってたオートバックスで付け替えした。
作業を待っている間、数学ガールを読んでた。4章のフィボナッチ数列と母関数の話はおもろいね。ちょっとのめりこみそう。
17122007 富士
義母が来たのでみんなで焼肉。前回失敗した店は避け新たな焼肉屋の開拓ということで、ふじみっていう、若干高級感のある外装の店を攻めてみた。
全室個室らしい。実際に部屋に案内されるとくつろぎの空間で満足。
肉にはいちいち花が添えてあって、とかいうまぁそんな感じ。たん塩と鶏軟骨が美味かったね。カルビとかは食後もたれるのであんま駄目です、むしろハラミをちょっと食べれば十分なんだけど。で、今回カルビを？って人数頼んだので結局もたれたうえ、ご飯ものに手がでなかった。
それが残念。
Mastering Algorithms With PerlModularとか。Chinese Remainder Theoremってなにがおもろいんだろうとか思ったが、RSA暗号がらみ？そのうちちゃんと読んでみる。
あとMiller-Rabinとかコラッツ予想とかそんな感じ。
perlで複素数とか極座標とか扱うって話 あまり興味がなかったので流した。
最後の例のリーマンゼータ関数は面白かった。
まだ2章くらいしか読んでない数学ガールにも出てくるらしい
これも早く読みたいなぁ。
15122007 chemoinformatics perl MAWP
さわりだけなので、さらっとしている。
Mastering Algorithms With Perlちゃんと学びたければC++とかの本を探したほうがよいのかな？pymolとかAvogadroのソースを読むのがいいのかもしれんなぁ。
ドッキングシミュレーションとか色々気に入らない部分があるので、手を加えたい部分があったりとか、Caverみたいなの自分で書いてみたかったりするんだけどそういう書籍はないもんかな。
後半はparserとか
Mastering Algorithms With PerlHOPも読んだし、再帰下降パーサもいじったことあるのでさらさらと。あとハフマン符号化とか。
結局1週間近くかけて読んだが文字列検索アルゴリズムがなかなか面白かった。
perlは写経でいまいち理解度があやしいため、pythonでBoyer-Moore法を。パターンとの不一致の場合には何文字スキップできるか、一致した場合に何文字スキップできるかのテーブルを最初に作っておいて、スキップ幅の大きいほうでスキップしながらパターンマッチを進めていく。
def boyer_moore_bad_character (P): m = len(P) bc = [m] * 256 j = m for char in list(P): bc[ord(char)] = j - 1 j -= 1 return m, bc def boyer_moore_good_suffix (m,P): j = m + 1 gs = [0] * j k = [0] * j k[m] = j for i in reversed(range(1,j)): while j <= m and P[i-1] != P[j-1]: if gs[j] == 0: gs[j] = j - i j = k[j] j -= 1 k[i-1] = j kn = k[0] for j in range(m+1): if gs[j] == 0: gs[j] = kn if j == kn: kn = k[kn] gs.pop(0) return gs def boyer_moore (T,P): (m,bc) = boyer_moore_bad_character(P) gs = boyer_moore_good_suffix(m,P) (i, last_i, first_j) = (0, len(T)-m, m-1) for i in range(last_i+1): j = first_j while j >= 0 and T[i+j] == P[j]: j -= 1 if j < 0: return i else: bcn = bc[ord(T[i+j])] - m + j + 1 gsn = gs[j] i += bcn if bcn > gsn else gsn return -1 if __name__ == '__main__': print boyer_moore("perlpypythonruby","python")
for文でカウンタ変数使いたかったのでrangeをつかったけどPython流はどうやるのがよろしいんだろうか、、、カウンタ変数減算するのにfor i in reversed(range(1,j))だしな。
パターンマッチングとか。
Mastering Algorithms With Perl色々と検索アルゴリズムが紹介されていてなかなか進まない。
Boyer-Moore法を写経して、ついでにPythonで書き直してみたけどgood suffixのアルゴリズムがうまくかけません。perlのfor文とpythonのfor文は大分違うなと。というか、文字列の位置を記憶しながら逆向きに移動していってごにょごにょする方法が分からん。
for(;;)みたいなこと(カウンタ変数を使うっていうの?)をpythonでやる場合にはどういうことをすればいいんじゃと調べてみたら、そういう場合にはfor文は使わないのがお約束としか書いてなかったが、どうもwhile文を使えということらしい。
ふむ。
perlのコード見ながらpythonのコードに移植しようとすんのが間違ってるのかな。
数学ガール (数学ガールシリーズ 1)
Mastering Algorithms With Perl
Mastering Algorithms With Perl