2013/01/17 (Thu) 淡路先生の力戦2局

第50期十段戦予選 高梨聖健 対 淡路修三

2013/01/17 (Thu) python: 組み合わせの計算（gmpyとscipy）

組み合わせ（n C k ）を計算する時間は、nが大きくなるにつれて大きくなる。一番時間がかかるのはkがちょうどnの半分のところ（階乗の計算が大変）。逆にk=1またはk=nの時は、n C k = 1 なので簡単。

pythonでn C k を計算する関数としては、gmpy.comb() とscipy.misc.comb() の2つが有名どころらしい（参考URL）。このうち、gmpy.comb()は正確に組み合わせの数を計算するのに対して、scipy.misc.comb()のほうは、対数ガンマ関数を利用して近似値を返すらしい。そのため、後者のほうがスピードが速いとか。

実験してみたところ、次のような結果になった。
n = 10000を固定して、k の値を1から10000まで、色々変えて、それぞれの関数でn C k の計算にかかる時間を記録した。グラフのy軸は経過時間（ミリ秒）で。対数にしてある。

予想通り、gmpyのほうは階乗の計算が難しくなる中間のkについて、計算時間は長くなり、真ん中（k=5000）のあたりで最大になる。一方、近似値を計算するscipyの方は、そういった問題が起こらないらしく、一定の速度で計算できるらしい。そもそも計算の簡単な範囲（kが1に近いかnに近いか）では、gmpyで計算したほうが速いということもグラフから分かる。

使用したファイル：comb_test.py, comb-graph.R

 

2013/01/16 (Wed) Latex: HTMLやワードに変換

文書はほぼ全てlatexで書く。構造が論理的なところが好き。

共著をする時には大きな障害になりそうな気がするが、そこは追々考えるとして、当面は、ライティングの先生に添削を依頼する時にワードに変換する必要が生じた。

文章だけが必要なので、方程式や図表は綺麗に表示出来なくても差し支えない。この目的なら、tth (HTMLへ変換)とlatex2rtf (RTFへ変換。ワードで読める)が良い。どちらもUbuntuの公式リポジトリに入っている。

特に、tthは知らないパッケージが出てきても、その部分を無視してとにかく変換を成し遂げてくれる、小粋なソフトウェアだ。つまりクラッシュしない（何も吐き出さないということはない）ので、ストレスが少ない。tthでHTMLに変換→ワードにコピペという方法もなかなか有効。

スタイルを忠実に再現したり、様々なパッケージに対応できるコンバーターもあるようだけど、今はとにかくクラッシュしないtthがベスト。

使い方：
tth <tex file>
latex2rtf <tex file>

【追記】
tthは日本語でもいける（Nice!）。アウトプットはオリジナルのファイルと同じエンコーディングになるようだ。

2013/01/16 (Wed) 情熱大陸

2013年1月13日「井山裕太」
面白かった。

http://www.mbs.jp/jounetsu/2013/01_13.shtml

2013/01/14 (Mon) 黒嘉嘉 in Ted x Taipei

どうやらこれがフルバージョンみたい。ずいぶん焦らすんだなー。