Tsuchiya Yoshihiro

オイラー解けてない問題

2012-02-04T21:46:00.004+09:00

このネタ。定期ポストです。

2008年の今頃からはじめたProject Eulerだが、今現在解けていない問題が41問。まる4年になるのか。40問というのが定数になっていて、残りの問題数はいつもそのくらい。問題は週1問出題なので、要するに週1ペースでしか解けていないということ。もう2009年の今頃からそんな感じ。そしてモチベーションを保つのも大切です。

分割統治法よろしく、40問を分割してみることにする。3つに分けよう。300番まで、350番まで、あと残り、で同じくらいづつある。まずは300番までを片付けようか。

163 Cross-hatched triangles 数え上げ。どう手をつけたら良いのかよくわからないのだよね。
194 Coloured Configurations 問題の意味がよくわからないし、絵の意味もよくわからないとおもっていたのだけれど、よく読んでみると、面倒くさいにおい。頑張ればできるんじゃないかという位置づけに変わった。
236 Luxury Hampers 事情はわかった。さてどうやって解こうか。
241 Perfection Quotients　σの記号がなんとなく恐ろしくて、考えていない。
248 Numbers for which Euler’s totient function equals 13!　これもなんとなく恐ろしくて考えていない。
254 Sums of Digit Factorials　これ、いままでに何度か挑んだのだけど、問題の最後を読む頃には最初の方を忘れていると言った感じ。最近集中して考えてみたが、すぐに解ける感じではなかった。
269 Polynomials with at least one integer root　これ、どういう問題なのかなあ。ちゃんと考えていない。
270 Cutting Squares　こういう問題、なんとなく嫌い、だった。今考えるとわかるのだろうか。
275 Balanced Sculptures　これも。まったく取り組んでいない。
280 Ant and seeds　この問題が出たとき、サンフランシスコのホテルにいて、問題の意味がわからんなあと、問題を音読したりした。その後、自分なりに取り組んだが妙な答えになって落胆した。
282 The Ackermann function　これ、chaemonさんの問題。最近読み直したけど、わりとストレートな問題なのね。どうやってやるのかなあ。
289 Eulerian Cycles　これ、たしかniinoさんの問題。ちゃんと考えたことなし。難しいという評判だったため。
292 Pythagorean Polygons　これは頑張ればできるんじゃないかなあ。
298 Selective Amnesia　これもniinoさんの問題。よくわからないのだよなあ。

こうしてみると食わず嫌いな問題が多いなあ。350番あたりからの問題はすごく難しい感じだけど、このあたりの問題はそれに比べると楽なのではないかと思って、やってみることにする。

Xbox360

2012-01-29T22:02:00.000+09:00

というわけで、Xbox360+Kinectがうちにあります。（前の項を参照）

ゲームのことは良く知らないです。つい先日まではうちにゲーム機はなく、クリスマスにサンタさんがニンテンドーDSを持ってきてしまったのが文明開化でした。

Xbox360と書きましたが、360の他にモデルがあるかどうかすら知りません。Kinectが付くモデルと付かないモデルがあるかどうかすら知りません。

このXboxにはソフトが２つついてました。1つはフィットネス。もう1つはアドベンチャーゲーム。

この一週間、フィットネスの方に取り組んでます。脚、筋肉痛です。スクワットとかその他の脚腰をいじめる系のメニューにいじめられてます。怠けているとお姉さんの声で「もっと低く！」とか言われるので頑張ります。30分くらいで体がポカポカになります。

しかしXboxには不満が少々。テレビにつなげるまでは簡単でした。ちょっとそのフィットネスとやらをやってみようとしたのですが、「アップデート」するのが先だ、というようなメッセージが出ます。Xboxとの押し問答を経て、結局Xboxを家のwifiにつなぎました。何やらXboxはやっていたのですが、私は「winrows update」のようなことをやっていたのではないかとにらんでいます。なにやらリブートらしきことを繰り返した後、使えるようになりました。

変な感じだったのはプロフィール作成。人の動きが分かって、スクワットの高さまで監視できるわけですから、いちいちプロフィール作成とかログインとかおかしくないか？

とにかくログインしなければならないのです。判断しろよ。機械でしょ、できるでしょあなた。スティーブジョブだったらこういうの、許さなかっただろうなあ。

文句は言ったけど、面白いです。画面がきれいです。というか、きれいとか思わないくらい自然です。筋肉痛です。インストラクターと並んで運動します。アドベンチャーゲームは子供達が喜んで使ってます。DSで体が固まっているのに比べると、だいぶ良いように思えます。

samurai coding

2012-01-23T23:02:00.001+09:00

samurai codingに行ってきました。早稲田の鷲崎研究室と、ゲームのGREEのコンテストです。会社の後輩の@chaemon_さんから「一緒に出ましょう」と言われたのですが、話を聞くと僕は何もしないで良さそうだったので賛成しました。オンラインで予選があって、決勝戦はGREEに招待されて会場で試合がありました。会場はなごやかなムードで楽しかったです。

結果は優勝で、Xboxを2台もらいました。私は貰わないつもりだったのですが、「テレビがない」というメンバーが多かったので、うちにも1台きました。

これは表彰式の映像。

政治について考える

2012-01-15T22:25:00.000+09:00

日本の首相は覚えきれないくらいに変わっているけど、まあどこの国に行っても、政治家の悪口を言っていればOKな状況のようです。どこもうまく行っていないですね。そこで、なんでこうなっているのかなあと考えた。100年単位くらいで考えれば何か見えるのではないか。

18世紀までは封建国家の時代。封建領主が人民を支配していた。
19世紀は国民国家の時代。国民が選んだ政治家と、国民が参加する軍隊によって国家が運営された。
20世紀はグローバル化の時代。前半は帝国主義でそれを実現しようとして失敗。冷戦が終った後半は、企業を中心とした力によってそれが行われた。なぜこうなったのか。鍵となったのは工業化、通信と情報の技術。

21世紀。これからはどうなるのか。どうあるべきなのか。どこの国に住んで、どこに行こうが自由の時代が来る。もう来ている。人は地理的な要因での支配は受けない。移動は簡単なほどよい。行った先でも同じ制度があるのが便利だ。まだそこまで出来ていないですね。これからのことです。そのときに、国境の意味は？税金や選挙の意味は？

共通化ということからすると、国に国民の口座があるような制度は不便だと思うんですよね。税金や年金の。保険が積み立て方式になっているより、掛け捨ての方がはっきりしていて良い、のと似ていると思います。

Haskellメモ

2012-01-07T17:33:00.000+09:00

思うところあって、Haskellをいじっている。いじり始めると案外面白い。しかしとっつき悪いですね。でもネットの情報は充実していて、本を買う必要なさそうです。

まずはFizzbuzzをやってみよう。3の倍数ならFizz, 5の倍数ならBuzz, 15の倍数ならFizzBuzz, そうでない場合はその数字を表示するというもの。

tsuchiya:haskell tsuchiya$ cat fizzbuzz.hs
fizzbuzz :: Integer -> Integer -> [[Char]]

fizzbuzz x y | x > y = []
| mod x 15 == 0 = "FizzBuzz": fizzbuzz (x+1) y
| mod x 3 == 0 = "Fizz": fizzbuzz (x+1) y
| mod x 5 == 0 = "Buzz": fizzbuzz (x+1) y
| otherwise = show x: fizzbuzz (x+1) y

main = print $ fizzbuzz 1 20

tsuchiya:haskell tsuchiya$ runhaskell fizzbuzz.hs
["1","2","Fizz","4","Buzz","Fizz","7","8","Fizz","Buzz","11","Fizz","13","14","FizzBuzz","16","17","Fizz","19","Buzz"]

こんな感じ。特にヒネリはありません。mod x 15はx `mod` 15と書くことができて、そっちの方が直感的ということになっているらしいけれど、私はキーボード上でバッククオートが見つけられない人なので、上のように書きました。

次、素因数分解をやってみよう。

tsuchiya:haskell tsuchiya$ cat factor.hs
factor :: Integer -> [Integer]

factor n | (mod n 2 ) == 0 = 2 : factor (div n 2)
| otherwise = fac n 3
where
fac n x | n*n < x = []
| (mod n x == 0) = x : fac (div n x) x
| otherwise = fac n (x+2)

main = print (factor 720)

tsuchiya:haskell tsuchiya$ runhaskell factor.hs
[2,2,2,2,3,3,5]

modの周りに付けた括弧はいらないんだろうなあ。まあいいや。otherwiseではなくて、Trueと言っても同じですね。そっちの方がスペルミスで怒られることがなくていいかも。

学生の頃、関数型言語とかλ計算とか言われても何も感じなかったなあ、よく覚えていないけど。

LED照明

2012-01-07T11:52:00.002+09:00

年末に、居間の蛍光灯をLED電灯に変えました。元は30Wくらいのサークラインが4つ入っているもので(120W)、今度のは60Wくらい（適当）なので、半分になります。電球を蛍光灯に変えたときほどではないと思うけど、一番使っている場所なので、効果を期待。1日6時間つけているとして、1ヶ月で60W*6*30=10KWhか。そこそこですかね。電気料金が変わるほどは大きくないか。

買ったのはパナソニック。他は東芝、日立、AGLEDっていう会社のもあった。AGLEDはロームのブランドだそうです。「自分で取り付けできますよ」と店で言われたので自分でやってみました。まあ結構大変なのですが、サークラインの取り替えはもうしないでよいのが良い点。

この一年をふりかえる

2012-01-02T23:19:00.000+09:00

去年はどんな年だったかなと思いました。普通は年末にこういうのをやるんですよね。まあいいや。

まず地震があった。3月の頭からナポリからホームステイの人が来ていたのです。金曜に地震があって、次の月曜には混乱の中、飛行機に乗ってイタリアに帰ってゆきましたので、その後のことを考えると賢明だったなと思います。

地震の時、家に連絡がつかなかったり、3時間歩きで家に帰ったりといったことは不思議な体験でした。うちの近辺ではその程度で済んだのですが。その後のトイレットペーパー不足と電力不足、この野菜は安全なんだろうかと思いつつの買い物など。

電力というのはコンピュータの仕事にとってはとても大事なのですが（しかもそれを沢山つかう）、それだけに我々は一体何をやろうとしているのかなあという気分になります。まさに地震のときに通話ができなかったiPhoneだって、元をただせば電力に依存してるのです。経済成長とか便利というのはもちろん重要なわけですが、原発みたいに一歩間違えば人命に危険を及ぼすものの上に成り立っているものと引きかえにしていいのかなあ。

本職の仕事では2月にUSENIX FAST (San Jose) でポスター発表をやっていい気分に、続いて5月にIEEE MSST (Denver) で発表してドキドキしました。まあ収穫期なんでしょうか。さらに会社の展示会で11月にミュンヘンで同じ内容を説明しました。こちらはビジネスの人向けなので、やり方はかなり変えましたが。もう何度も人に説明しているので同じ話を説明するのはちょっと飽きた。そろそろ新ネタを作りたいですね。

Project Eulerは週に一問ペースでダラダラ解いているような感じですが、私はこれ以上速くなりません。2008年に始めましたが、ここ2年くらい常時残り40問を残す感じです。いつもながら苦しいです。2度目になりましたが、twitterで声をかけて10人ほどの方々と飲み会をしました。遠い人は京都や浜松から来られまして感激しました。また春ごろやるつもりです。

「ジョブス本」を読みました。彼は「自分が納得できる、素晴らしいモノを作る」ことが目標だったようです。そう有りたいですね。

Parallels for Mac

2011-12-31T01:23:00.002+09:00

調子に乗っています。

家でWindowsを使いたいこともあるので、ParallelsとWin7を買いました。

Parallelsのインストールは解くに問題なし。Windowsはちょっと難しかったです。

「無人応答ファイルに無効なプロダクトキーが含まれています」

という意味不明なメッセージで止まってしまうのです。ネットの情報だと、「共有」が良くないのだと。

Windowsの設定から共有と名のつくものを手当たり次第にoffに（あとでonにできます）。インストールできました。

Windowsのインストールパターンで、Mac向き（？だったかな）を選びました。VMwareみたいにOSの窓がでてくるんじゃなくて、MacOSXのドックの中にWindowsのスタートメニューボタンが出来てる。なんだかすごいです。Internet Exploreもドックに入ってきます。

このwindowsを子供に使わせたいのですが、共有で使うには「共有フォルダ」にWin7を置かねばならず、そこにはWindowsを入れる64GBもないので今日のところは仕方なく共有でないパターンでインストールしました。Windowsは昨日作った外付けRAIDに入れたいのです。どうするれば良いのかなあ。子供用Win7は別にするべきなのでしょうか？

-----
(追記)
わかりました。少々ムリなやり方ですが、Parallelsのディレクトリで、パーミッションを全員に出す作戦です。Terminalからchmodコマンドで設定を変えただけです。子供が使っても私が使っても同じWindowsが立ち上がりますが問題ありません。

おうちRAID

2011-12-29T21:23:00.000+09:00

このように考えて、家にRAIDを組むことにしました。買ったのはラトックのRAIDケース↓と、AmazonオススメのWestern digitalの WD20EARX。タイの洪水のせいでディスクは高騰しているようです。安いときの2倍しているんじゃないかなあ。とはいえ、せいぜい1万円の違いなので、買っちゃいました。

ディスクは2TBが2本で、ミラーしているので　2TB。RAID1の設定したときちょっと危うい感じだった。ディスク故障の表示。ごちゃごちゃやっているうちに設定できてしまった。なんだかあやしい感じ。

これって、ディスク一個が壊れてもOKなのがミラーなんだけど、故障のディスクを交換したときにちゃんと修復するんだろうなあ。当然。少し不安です。

MacOSX Snow Leopardのマシンにつなげました。ついでにTime Machineの設定も。そこらにあった写真やビデオのDVDやらCDやらを読み込んでも余裕です。

ちょうどよい大掃除という感じで、要らないDVDが捨てられます。

うちのデータに関するあれこれ

2011-12-25T07:49:00.000+09:00

いまうちでは2008年に買ったMacbookに外付けディスプレイとHappy Hacking Keyboardをつないで子供やおくさんが使っている。データのバックアップはDVD。ハードディスクレコーダーも2008年頃のDIGAで外部とのインタフェースはなく、DVDを焼く以外の方法でデータの交換はできない。ハードディスクレコーダーの中の、家族が映ったビデオなどはDVDにも焼いてある。大事なデータは写真とビデオが主で、最近ではあまり増えなくなった。つまりデータの行き先としてはDVDである。他にはDropboxとPicasawebを使っている。私が使っているMacbookAirのデータはDropboxに入っている。

DVDを中心にしたバックアップのスキームはDVDの信頼性に依存してるわけだが、それもそろそろ心配するようになってきた。そろそろMacbookかDVDのいずれかがダメになる可能性があり、ハードディスクレコーダーも一杯になっているし、なにしろアナログなのでそのうち買い替えるだろう。DVDをいっぺん読み取っておく必要があるかなと。

MacbookにつなげられるディスクにDVDを読み取ってみようか。DVDの枚数は少なくはないが、それほど多くもない。2,30枚といったところである。一枚4GBとして120GB程度。データ全体はMacbookに入っているのが50GBとしても全体で200GBである。MacbookがダメになったらMac Miniを買おう。うちに200GB入れる場所はないので、いずれにしても新しいディスクは必要。

200GB預かってくれるクラウドのサービスってあるのかな。その場合でも手元にディスクは必要だなあ。

ためしにビデオのDVDを一枚読み取ってみたところ、30分ほどかかった。4GBを30分で読み取るのは2MB/sec程度か。ゆっくりだなあ。コンテンツは2時間ほどあるだろうからそれでも4倍速くらいで読んだことになる。ビデオのスピードってこんな遅くっていいのだっけ。動画再生に必要なディスクのスループットって100Mbpsで十分なんですね。20枚読み取るのに10時間か。ちょっと面倒だけど、何日かかければ可能な範囲。

Macユーザとしては、Airmac Extremeとそれにつながるディスクに全部入れてTimeMachineにするのが良いのだろうか。ディスクの信頼性も考えると、RAIDになっているのが良いのだけどなあ。単純にミラーが良い。

あんまり高いものを買うのもイヤなので、ワイヤレスにするのはあきらめてMacbookにつながるUSBでRAIDになってるディスクか（RAID HDDケース）、有線のLANでつながるNASにしようかな（省電力なサーバでFreeNASもいいかも）。

AirMac Expressでプリンタをネットワーク共用

2011-12-24T17:21:00.000+09:00

みんなには一切関係ない自分のメモ。

うちのプリンタFAX複合機「見楽る」だが、いままで古い方のMacにLAN直結でつないでいた。その前はYahooBBから独立したIPアドレス（グローバルなやつ！）をもらってネットワークプリンタになっていたのだった。LAN直結になったのは、NTTの光になって、まともに家に一個のIPアドレスしか貰えなくなったため。マシンが一台しかないので別に問題なくLANでつないでいたのだった。

AirMac Expressの口にUSBのコネクタがあって、そこにつないでいた時期もあったんだよなあ。どうしてやめちゃったんだろう。とにかくプリンタから垂れ下がっていたUSBケーブルをAirMacにつなぎ直して「見楽る」は再びネットワークプリンタになった。

「スティーブ・ジョブズ」

2011-12-23T20:33:00.000+09:00

「ジョブス本」を読み終わりました。ジョブスの考えは、一番良いもの (great product)を作ることにあったようです。そのためにend-to-end controlが必要だったということです。つまりハードからソフトまで一貫してアップルが提供すること。さらには製品をアップルストアで
売ることが必要だった。会社が大きくなると、金儲けに執着してしまうものだけれど、ジョブスは製品にこだわり続けのだそうです。

私が読んだのは英語のKindle版のもの。頑張って読みました。ジョブスはたぶんKindleデバイスが気に入らなかっただろうなあ。

ミュンヘン

2011-11-19T12:07:00.001+09:00

11/6-12の間、仕事で行ってきました。空港からホテルに行くまでの道の暗いこと。ハイウエイらしき通りに街灯はありません。原発事故以来、東京が暗くなったとお嘆きの方々は、いちどミュンヘンをご覧になってはいかがか。町中もいたって質素です。光る看板はほとんどなく（古い建物を大事にしているのでしょう、おそらく付けることができないのでは）、落ち着いた町並みです。人が着ているコートの色も黒か茶なのですよね。男も女も。これはちょっと地味すぎるかなあと思いました。

ビールはおいしかった。ビヤホールは2カ所ゆきました。普通のビールと白ビールをためしたけれど、普通のものの方が私には良かったです。

空いた時間にドイツ博物館に行って、昔の計算機を見ました。IBM System 360だとかCRAY-1などがあるのです。ところで、博物館に時計台があって、時計なのかなと思ってみると、針は一本で、「バロメーター」って書いてあるんですよね。調べてみると、説明があった。

2平方数の和（２）

2011-10-15T09:45:00.001+09:00

nがあたえられたときに、x^2+y^2=nを解くにはどうしたらよいか。nが素数pの場合については以前書いた。では合成数の場合はどのようにすればよいか。

n=pqのとき、以前の方法でp=u^2+v^2, q = a^2+b^2ということが分かっているものとして、つぎの恒等式をつかってnを平方数の和に分ける。

n = pq = (u^2+v^2)(a^2+b^2) = (ua+vb)^2 + (va-ub)^2

これでnはua+vbとabs(va-ub)に分けられた（どちらも正数だけ採用することにする）。これを繰り返してゆけばnが2個以上の素数で分解出来る場合も対処できる。

元々、pが4で割ったとき1余る場合、x^2+y^2=pの解はあり、pの余りが3の場合に解はない。nにそのようなpが含まれている場合、pの個数が奇数のとき、x^2+y^2=nに解がなくなる。偶数の場合には解ありである。例えば、n=9の時、p=3だが3が2個あるのでx^2+y^2=1の解 (1,0)を3倍した(3, 0)が解となる。また、n=3やn=27には解がない。

コードにしたのが以下。素因数分解は状況によって違うものを使いたい（試し割りのものと、フルイをつかったもの）ので、初期化で渡すようにした。以下ではfacs()がそれで、これに期待しているのは、n=p^r * q^sのとき、[p,r],[q,s]を渡してくれること。fscsの中で使っているfactors()が素因数分解のコア部分。

名前をprocess_primeに変えたが、これが以前書いた、素数のx^2+y^2=pを解くところ。comb()が上に書いた恒等式で展開するもの。

class square_sum(object):
def __init__(self, func):
self.func = func
def square_mod(self, p):
from random import randint
while True:
a = randint(2, p - 1)
b = pow(a, (p - 1) / 4) % p
if (b*b) % p == p - 1:
return b

def process_prime(self, p):
if p == 2: return [1, 1]
if p % 4 != 1: return [0, 0]

a = self.square_mod(p)
b = 1
k = (a * a + b * b) / p

while k > 1:
u, v = a % k, b % k
# abs(u) <= k/2
# abs(v) <= k/2
if u > k / 2: u -= k
if v > k / 2: v -= k

# u * u + v * v = 0 (mod k)
q = (u * u + v * v) / k

# (u * u + v * v) * (a * a + b * b) = (u * a + v * b)^2 + (v * a - u * b)^2 = (q * k) * (k * p)
c, d = (u * a + v * b) / k, (v * a - u * b) / k

# next turn
a, b = c, d
k = q
return [abs(a), abs(b)]

def comb(self, X, Y):
u, v = X
a, b = Y

res = set()
r1 = [u * a + v * b, abs(v * a - u * b)]
r1.sort()
res.add(tuple(r1))

r2 = [u * b + v * a, abs(v * b - u * a)]
r2.sort()
res.add(tuple(r2))

return res

def times(self, M, seq):
return map(lambda x: x*M, seq)

def solv(self, n):
fac = self.func(n)

# p%4==3 prime
M = 1
for p,r in fac:
if p % 4 == 3:
if r % 2 == 0:
M *= p**(r/2)
else:
return []

# p==2 or p%4==1 prime
seq = []
for p,r in fac:
if p == 2 or p % 4 == 1:
sq = self.process_prime(p)
seq += [sq]*r

if len(seq) == 0:
return [self.times(M, [0,1])]

res = set()
res.add(tuple(seq[0]))
for k in xrange(1, len(seq)):
new = set()
for e in res:
new |= self.comb(e, seq[k])
res = new

return map(lambda x: self.times(M,x), res)

# if n=p1^r1 x ... x pm^rm, returns [p1,r1]...[pm,rm]
def facs(n):
seq = factors(n)
s = set(seq)
r = []
for p in s:
r.append([p, seq.count(p)])
return r

使い方は、
sq = square_sum(facs)
sq.solv(1234)
のような感じ。
比べてみると、しらみつぶしでx^2+y^2=nを解くよりはかなりはやい。

禅について

2011-10-14T23:59:00.000+09:00

私と禅について書いてみる。始まりは15年ほど前。禅というものにはちょっと興味があって鈴木大拙の岩波新書の禅の本などを読んでみたが何のことか全然わからなかった。茶道にも興味があったので表千家の先生について習った。茶道の掛け軸と言えば、禅僧のものが定番なのでここでも禅にたいする疑問がわいた。禅とは何なのか？しかしどうやら本を読んでも分からないだろうなと思っていた。

松門寺

松門寺に行ったのは1998年のこと。インターネットで見つけた。松門寺は曹洞宗のお寺。ここの坐禅会は月一回、2,30人の人が集まり、坐禅の時間、正法眼蔵の講話の時間があった。お寺は住職である老師が一人で取り仕切っている。坐禅堂はモダンな作りで新しく清潔だった。坐禅の時間には、独参（どくさん）がある。独参では老師の部屋に一人一人入り対座し、坐禅についての質問や見解を述べる。坐禅がはじまりしばらくすると、独参開始の合図が鐘の音で知らされ、坐禅をやめて独参の部屋の前に並ぶ。

坐禅とは仏教である。インドからやってきた達磨が唐の時代の中国に伝えたということになっている。老師の言い方では、坐禅とは「自分自身になる」ことであった。これは指導者の説明の仕方により、いろいろな言い方がある。ある人は「即心是仏」と言った。言葉で何と言っても当たらない。言葉にはとらわれないようにとのこと。

「自分自身」、それ以上のことはない。また、従うべき指針や善悪の基準があるわけではない。つまり禅には教義がないのである。「自分自身になる」にはどうすればよいのか？「何もしなくてよい」と老師は言う。何もしなくてよい？一体どういうことなんだ？「坐禅のときに、集中ができません」、それが悪いことだと誰が決めたのか？集中すべきと思ったのは自分。そこから離れてごらん。集中できない自分を自分として受け入れる。それが坐禅の一部。という調子。一筋縄ではないのです。

老師は東工大の化学を出た人。答えが分かるのは答えが分かったとき、という事の性質上、先生が信用出来るかどうかが鍵となる。こちらも理系だからということもあったのかもしれないが、この人が有ると言ったものは本当に有るのだろうし、出来ると言ったことは出来るのだろうと思った。

それにしても、ちょっとした疑問を解消しようとしているのにしては、何だか遠くにきてしまったなあというのが私の気持ちだった。

私の独参

私の独参は、最初のうち質問。松門寺に行けば毎回行っていた。3回くらい行ったところで、もう質問することがなくなった。「公案をやってみるか」と言われた。公案とは坐禅の先生から弟子への質問である。禅宗でも曹洞宗は公案をやらないとされていて、臨済宗では公案をやる。しかし曹洞宗も臨済宗も本は同じ、違いはない、曹洞宗ではあるが、うちでは公案をやりますとのこと。次に来るときにはこれに答えてくださいという意味のことを言われた。これができたら「悟ったことになるのかな」などと思った。一ヶ月後、何と答えたのか忘れたが、それは違うね、という感じのことを言われた。公案はテストであり、悟っているかどうかがたちどころにわかるのが公案だとのこと。何なのだろうか。その答えは、本の後ろの方に「解答」が書いてあるような類いのものでない。

「誰にでもでき」て、それをするのには「何もしなくてよく」、「いつ分かってもおかしくない」と言われればいてもたってもいられないだろう。絶対に分かってやる、と思っていた。

そういう独参が2年くらい続いた。一ヶ月ごとに答えを持っていって、否定されるまでに10秒とかからない。方向は良くなってきているけれどね、など。その後ちょっと雑談するなどして退席である。一体オレはなにをやっているのかなあというのが正直なところだった。老師が言うところによると「100メートルの棒の上を歩いてきた。その先っぽにいる。お前はそこから一歩を進めないでいる」とか、「柿は熟すれば自然に木から落ちる。それなのに、お前は自分の力で木にしがみついているだけだ。なんで自分の力を使おうとするのか」など。これらの言い回しは、古来の公案にある言い方である。

さて、およそのことは分かった。しかし一体なんと言えば通してくれるのだろうか。今にして思えば、そう思っていること自体に問題があったのだろうなあ。あるとき、「一人通った人がいる」、と言うことだった。うらやましいなあと思った。「他にも通りそうな人がここには何人かいる。」それ、私のことなんだろうなあ。

最後には、この他には何もない、棒の先にいるなら先にすすもう。通ろうという気持ちで独参に向かった。いくつかのやりとりのあと、老師は言った「これからもしっかり坐ってください。」

これが2001年の夏のこと。この最後の独参のあと、老師のところには行っていない。

肝心なところは、禅は紙に書いてある知識ではなく、自分自身との対話である点。本をいくら読んでも分からなかっただろうと思う。また、チェックしてくれる指導者が必要である点です。

メガネ

2011-10-10T22:31:00.002+09:00

前回はメガネ屋さんにお世話になったのはこのブログの検索によると、3年ちょっと前。駅前の個人経営のお店でした。子供が別のお店でメガネを作ったので、私もそちらに行ってみました。だいたい、あのお店に行く時点で客層は限定されます。今時ははるかに安い店があるので、あの手の店に行く人は、メガネにならお金を出すよという人なんでしょう。フレームも高いし、店員が提案してくるレンズの値段もビックリです。説明は慇懃な感じ、こちらが選んだものについて、向こうは何も言いません。店員によるかもしれませんが、もしかしたらその店員に技術はないのかも。どうなんでしょうね。10通りくらいのオプションの中から下1/3くらいの位置にある選択肢を選びました。前の駅前の店では、こちらが何も言わなくても、あらかじめ決まっている感じで、フレームを決めれば後は何も言わなくても不思議な割引で値段が決まっている感じでした。値段は前の方が安かった。

Google code jam, Japan 2011

2011-10-02T23:14:00.001+09:00

Google code jam, Japanというイベントがあって、そういうのがあるのは知っていたのだけど、それってすごく出来る人向けのものだと思っていたので、やってみようなどと思ったことがなかったのでした。今回のこれは日本在住の人向けとのことで、だから比較的まともな時間（夜中とかじゃない）に開催されるという話は聞いていました。土曜が予選で、土曜の昼過ぎにtwitterを見ていたら、どうやら始まっているらしいことがわかり、問題を見てみようと言う気になりました。時間は1時から6時間の時間制限、一問正解（smallとlarge)で予選通過とのこと。見ると、問題は3問。Aはトランプのシャッフルの問題。Bは何かを最大化する問題。Cは2進数で1が立っている数の問題。見た段階で、AとCは正解者がいました。Aはデバッグが面倒くさい匂いが、Cは単純そうな感じがしたので、Cを考えてみることにしました。いたずらでpythonでプログラムしてみるとできた様子だったので、まじめに考えることに。オイラー以外に参加したことがなかったのでルールを読みます。small、largeの順に解く。それぞれ小問からなっていて、smallよりlargeの方が小問数が多いという関係。largeの方が時間がかかるので心してとりかかれということのようです。問題の入力をダウンロードするとカウンタが動き出してsmallは3分以内、largeは8分以内に解答しないとダメよ、とのこと。どんなサイズか見当がつかずpythonでは不安なのでC++で実装しました。smallをダウンロードして実行。ノータイムで答えはでます。解答すると、答えが間違っていると。答えそのものではなく、小問それぞれにCase #1とかCase #2とかつけて出力するところを、私はCase #0と0始まりで書いていた。それはそうだろうと直して再度提出。プログラムも提出するんですよ。間違えるとペナルティがあるらしい。次、largeをダウンロードして実行。これもノータイムで答えが出る。今度は正解。特に結果を味わうこともなく、3問中1問おわったので私は満足して家事に。買い物など。A問題も考えてみると面白いかも。皆さん全問正解をめざして頑張っておられた様子。これが予選で、本選は次の土曜とのこと。時間は同じなのかな。でも来週は運動会の日。私は小学校のグラウンドにいる筈です。

台風

2011-09-22T22:25:00.000+09:00

台風15号が関東を直撃しました。
311地震に続いて、またしても帰宅難民になってしまいました。反省してます。

12時頃はひどい雨で、外に出るのはためらわれました。そんな日にかぎって、仕事する気満々で、16時くらいになったら帰ろう、などと思っていました。会社から帰宅するように言われたのが15時ころ。駅に行くと、ホームから人がこぼれるほどの混みかた。電車は間引きで、来る電車は既に相当に混んでいるのです。並んでいる人の3列め位までしか乗れません。一本やりすごして、乗ったのは2本目の電車。30分待ちです。これ、震災後と全く同じです。

乗り換え駅に行くと、乗り換える私鉄は既に止まっていました。家に帰るにはこの電車がたよりなので、万事休すです。そのあとビール飲んだり、立ち読みしたり、食事したりと時間をつぶして、8時半くらいに動き出した電車で帰宅。9時くらいでした。

台風と地震の違い。
地震と違って、台風が来る時間は予測できる。もっと早く帰れば良かったのです。それだけです。私は何を待っていたのでしょう？
台風の場合、屋内にいれば良いのです。ビルの中にいれば安全です。歩いて帰ろうなどと思わない方が良い。

やはり、人がたくさんいる場所がこわい感じがします。ホームで電車待ちしてるときが一番あぶない感じでした。

半年前の経験とか教訓が全然生きていないなあと反省してます。会社の帰宅指示のタイミングが最悪だったということはあると思いますが、特に拘束されていなかったので、自分の責任です。もっと早く帰るべきでした。

説明のつかない日本語#2

2011-09-07T22:14:00.000+09:00

ここに引き続き、「説明できない日本語」を収集しています。
今回のは、twitterで教えてもらったもの。NHKの日本語講座の題材だそうです。

「飲みに行きませんか？」に対し、「結構です」だとノー、「結構ですね」だとイエスの意味になる

うん、ここまで来ると、日本語ってクレイジー、という感じがしますね。

食塩水の問題

2011-08-30T09:06:00.002+09:00

子供の算数を見ていて、ちょっと困った問題。

6%の食塩水100gに、15%の食塩水を混ぜたところ、濃さが10%になりました。15%の食塩水を何g混ぜましたか？

中学生以上なら、

6+0.15*x = 0.1*(100+x)

として、

x=80

とします。私は小学生のときからこうやっていましたが。

小学生の解答は以下のように考えるようです。
薄い方の食塩水の水と、食塩を取り替えて、濃さが10%になるようにするには？
逆に濃い方の食塩と、水を取り替えて、濃さが10%にするには？

つまり

(0.1-0.06)*100=4

と

(0.15-0.1)*x = 0.05*x

が等しいので

x =80

と考えるのです。

少し前に、方程式をたてないのもいいなあと思ったばかりですが、今回のこれはやっぱり嫌だなあ。問題によってやり方を変える感じが気にいらない。これだったら最初から方程式を教えた方がいいんじゃないかと思います。

夏風邪

2011-08-30T09:00:00.002+09:00

3週連続で風邪になって、週1、2日休んでいます。風邪といっても咳がでるわけではなく、ノドがはれて高熱がでるもの。何年もこのパターンはなかったのだけど、久しぶりにやってしまいました。前回は5年くらいまえだったかもしれない。お医者さんの診察室で抗生物質の点滴をされて、目の前が暗くなって倒れた。

今回は、「治った」と思って仕事に行くと風邪をぶり返す、の繰り返しで、3回目になったのでさすがにまずいと思って、今日は家で静かにしています。

MacOSX LionにScipyを入れる

2011-08-13T16:25:00.002+09:00

LionにScipyを入れようとして、macportsを使ってみたのですがいろいろなパッケージをインストールしてくれている様子でしたが、うまく行きませんでした。macportsは/opt/の下にパッケージをおくはずなのだけれど、それらしいのがない。matplotlibは入ったかに見えたけれど、matplotlib.pyplotがないとか。何がおきたのでしょう。

numpyはlionに入っていたのかなあ。それともこの前gnuplotを入れたついでにpython2.7とともに入ったのか不明でございます。

教えてくれた人がいました。scipy-superpack。誰かがパッケージにしてくれてるんでしょうね。ダウンロード元がどなたかのdropboxのフォルダでビックリしました。こういうもののメンテってどうなっているんだろう？まあいいやとりあえずグラフがかけるようになった。

全角で入力をお願いします

2011-07-31T17:52:00.002+09:00

会員登録のような画面で、住所を入力するときに「全角で」といわれることがあるのですが、あれは何とかならないのでしょうかね。まず理由がわからない。全角で入れようが半角（というのもおかしい、Windows用語か？）で入れようが、数字は数字で、文字として入力されているのならばシステムで変換可能なはずです。またWindowsには入力方法で全角文字半角文字という選択があるけれど、Macにはそれがなく、全角英数を要求されると困ってしまうのです。この不便な仕様がはびこっている理由を是非知りたいと思っています。

場合によっては半角カタカナを要求される場合があるのですが、一体何のつもりなのでしょうか。フリガナは、ひらがなでもカタカナでもローマ字でも問題なく処理でき、データベースの内部で都合の良いものにおきかえれば良いのではないでしょうか。人間が入力するのが半角である必要など全然ないと思います。

約数の個数

2011-07-29T22:32:00.003+09:00

子供が、100までの数で約数が一番多いのは何だろうかと聞いたので、考えてみた。

6の約数は1,2,3,6の4つ。12の約数は1,2,3,4,6,12の6個。これは素因数分解が関係していて、6=2¹3¹、12=2²3¹であり、それぞれの乗数に1を足したものを掛け合わせたのが約数の個数になっている。2を選択する個数（0個も含めて）は乗数に1を足したものだからだ。6の場合は(1+1)(1+1)=4、12は(2+1)(1+1)=6になっている。

さて100までの数で、乗数の多いものといったら、64=2⁶だが、これは(6+1)=7個にしかならない。2x3の倍数だと、96=2⁵3¹で(5+1)(1+1)=12個が一番多そうだ。素数が3つのものでは60=2²3¹5¹で、おなじく(2+1)(1+1)(1+1)=12。

あとは面倒なのでコンピュータを使って試し割りで確かめると、12が最大で、12になるものは他に3つあって
72=2³3²で、4x3=12
84=2²3¹7¹で、3x2x2=12
90=2¹3²5¹で、2x3x2=12
だった。

というところまでわかったのだが、子供はもう寝てしまった。

Macbook Air(3)

2011-07-28T05:40:00.001+09:00

Xcodeもインストールした。前はXcodeにcscopeが入っていたように思うけど、今はないみたい。