男女各8人で合コンしたとする。 そこで女がそれぞれ気に入った男を2人選んだとする。 その選び方が全くのランダムなとき、 8組のカップルができる確率を求める。 最初に考えた問題は、数学的には同じだが、上と違う例だった。 しかし、グラフ理論ではこうい…

雨天中止。

最後ぐらい歩いて帰りたかったが、 なかなかそうもいかず。 関ヶ原以来まともに歩いてない。

なかなか机が片付かない

手続きのため住民票取ったりと。

6:08 刈谷 〜 大垣 7:12 大垣 〜 揖斐 〜 （揖斐川マラソンコース） 〜 大垣 Dst 80.93 Mx 47.0 Tm 4.43.34 Av 17.1 まあ色々あったのだが、 「鮎は美味しかった」 とだけ記しておこう。 16:55 大垣 〜 刈谷 17:42 実際には3分遅れ

その話はお断りだな。

5ヶ月ぶりに輪行袋にしまおうとしたら、なかなかうまくいかない。 結局12時過ぎに就寝。 明朝は5時半起きなのに。

予定通り蟹江までなんだとか。

一番高いところの標高は170mくらいなんだそうだ。 折り返し地点が違ってたらしい。

なんかロボットがどーとか。

熊谷さんから土曜日の計画が。 大垣駅から揖斐川を登る。 標高500ｍくらいあるらしい。 そっから下って、蟹江まで。 117kmとか。 大垣まででいいかも。

Dst 59.74 Mx 42.0 Tm 2.54.11 Av 20.6 まだ直していない。

Dst 56.89 Mx 47.5 Tm 2.50.59 Av 20.0 調節していないので、3%くらい距離長めのような気がする。

http://d.hatena.ne.jp/inamori/20050914で書いたのは練習問題。 本当に知りたかったのは次の確率。 8つのコップがあり8人がボールを2つずつその中に入れる。 ただし、一人が同じコップに2つボールを入れてはいけない。 そのとき全てのコップに2個ずつボール…

手続き開始。

1:05 米原 〜 関ヶ原 5:25 徒歩 約21km？ 5:55 関ヶ原 〜 刈谷 7:19

13:39 刈谷 〜 米原 15:10 15:13 米原 〜 彦根 16:20 徒歩 約6.5km？ 16:30 彦根 〜 大阪 17:44 17:54 大阪 〜 本町 18:15 軽くジョギング 約2.7km 20:50 本町 〜 梅田 21:30くらい 22:36 大阪 〜 京都 23:31 23:55 京都 〜 米原 1:02

昨日は擬似乱数で失敗したという話をした。 0〜7をランダムに発生させるのに、 cstdlibのrand()を使って、 rand()&7 と下の3ビットを使ったら、 ぱっと見はそんなに問題ないけど、 やっぱり多少ホントのランダムよりばらついちゃうんだね、 確率が大きくなっ…

渋滞の中、軽く走って駅へ。

昨日そのプログラムを作って、 念のため乱数を使ったプログラムを作って検証してみたら、 かなり違う値が出たので更新を保留していた。 よくよく見てみたら、 どうやら擬似乱数の使い方が間違っていたらしい。 擬似乱数は恐い。 結局、昨日の確率は、0.029%。

8つのコップがあって、ボールをランダムに16個コップに入れる。 そのとき全てのコップに2個ずつボールが入る確率は？ 難しそうだが、ここはコンピュータの力を借りて計算しよう。 といってもしらみつぶしは無理だから、 多少工夫してみる。 状態を次のように…

23:05 刈谷 〜 知立 23:11 23:13 知立 〜 東岡崎 23:24 0:20 東岡崎 〜 東京駅丸の内南口 5:13 夜行バス 約25km 東京駅 〜 石川島公園 〜 月島 〜 晴海 〜 豊洲 〜 夢の島 〜 葛西 〜 浦安 〜 市川 〜 二俣新町駅 11:52着 11:58 二俣新町 〜 海浜幕張 12:08 …

ルジャンドル多項式の次にチェビジェフ多項式が出てきた。 これは、 と表せるもので、 例えば、 だから、 cosnθの計算って高校のとき誰もがやったと思う。 なかなか計算が大変で、5次くらいまでしかやった記憶がない。 ここではなるべく高次まで計算したいと…

10:29 刈谷 〜 豊橋 11:00 実際には3分遅れ 11:10 豊橋 〜 目的地 11:07 6kmくらい？ 食事してるうちに自転車用意してもらって、 裏に回って、おやっさんいるかなと思ってのぞいてたら、 みえて、 お礼を言って、出発。 とにかく手が痛い。 左がまだいいので…

1500m 6分15秒 サボっているとこんなもんか。

と書くと、 これで計算すると、やっぱり31次までしかルジャンドル多項式の係数が出なかった。 なんだ、苦労したのに。 これを元に多項式をmimeTeXで表示するプログラムを書くと、

海に向かって逃げる

昔作った64ビット整数を用いた多項式の計算を行うライブラリを使って、ルジャンドル多項式を計算してみよう。 単純に、 n乗してn回微分して割り算して、 だとすぐにオーバーフローして、13次までのルジャンドル多項式しか出てこない。 ちょっと工夫して、 1…

というもので数式が書けるらしい。 これはルジャンドルの多項式。 数値計算の本を読んでいたら出てきた。 懐かしいね。 確か、高校でもちょこっと出てきた。 ヒマがあったら計算してみよう。