2008-01-01から1年間の記事一覧

30kmLSD

この地方では、今日はフルマラソンの日のようだが、私にはまだフルを走りきる体力はないので、30kmを。別に日をあわせたわけじゃないけどね。 LSD(Long Slow Distance)というのは定義はよくわからないが、ゆっくり走ると自転車と同じで体がブレやすい、だ…

アクセス時間が一定でない配列の2分探索(1)

0〜1024までのラベルがついた箱が左から一列に並んでおり、それぞれの箱には0か1の数字が書かれたボールが一つずつ入っている。0の箱には0のボールが、1024の箱には1のボールが入っていて、右隣の箱のボールの数字は左隣の箱のボールの数字より必ず大きいか…

第29回篠山ABCマラソンエントリー

巷では、今日は東京マラソンの当落で悲喜こもごものようだが、私は今年に引き続いてこの大会で勝負しようと思う。あまり面白いコースではないけれど、日帰りできるのはコスト面で大きい。 予想タイムを3:28としておいた。これで自分にプレッシャーがかけられ…

多項式の零点(2)

零点が存在する円の半径Rに対して、複素平面上で原点を中心とした-RからRの正方形を考えて、それを縦横10分割した格子点から最急降下法で零点を求める。 これを、f(x) = 16-52x+48x2+68x3-30x4-39x5+7x6+4x7-10x8-2x9+6x10+x11-x12 に対して適用したところ、…

第27回半田市民マラソン 大会結果

まだランネットのほうのみ。 46位/181人かあ。

複素数

Pythonで複素数は次のように使う。 z = 1 + 1j print z * z # 2j虚数単位は「j」で表す。ただし、上のように単に「j」とは書けず、数値が前に必要。 各種関数はcmathモジュールにある。 from cmath import *z = 1j print sin(z) # 1.17520119364j実部と虚部…

多項式の零点(1)

f(x) = a0 + a1x + anxn の零点を探索する。 その前に、補題を考える。 f(x) = anxn - (a0 + ... + an-1xn-1) n ≧ 1 ai ≧ 0 an > 0 は、x > 0 で f(x) > 0 なら x < x' で f(x') > 0 xf'(x) = nanxn - (a1 + ... + (n-1)an-1xn-1) > (n-1)a1 + ... + an-1xn-…

第27回半田市民マラソン大会

あまりうまく走れなかった。力を出し切れなかった。 それにしても、こんなに厳しいコースだとは思わなかった。 電車で半田へ。途中、大府で乗り換え。そこで短パンを忘れたことに気がついた。何やってんだか。すぐに電車が来たのでそれに乗って家に戻ると、…

listの加算代入演算子のオーバーライド

こんなコードを実行する。 class Poly(list): def __add__(self, a): res = Poly([]) m = len(self) n = len(a) if m res[0:n] = a[0:n] for i in range(m): res[i] += self[i] else: res[0:m] = self[0:m] for i in range(n): res[i] += a[i] return resf =…

今日は1日雨降り。土日はなるべく30kmは走りたいので、雨がやむのを待っていたが、なかなかやまない。結局、夕方にだいたいやんだので、そこから走り出した。10kmの中盤にはもう暗くなっていた。暗いところを走るのは危ないからなるべくやりたくない。車が何…

ランステ

皇居外周コースを走るためにランステを利用した。いつもは皇居の近くのホテルに泊まって朝早く起きて走りに出かけて帰ってきて風呂に入って10時間際にチェックアウトする。しかし、今回は夜行で来たのでそうはいかない。 皇居の外周コースは約5kmだが、5kmに…

行列式の計算

Gaussの掃き出し法のような感じで行列式は計算できるのだが、行列が大きくなると非常に大変。かつ、疎な行列の場合は、余因子展開したほうが速い。

シューズに慣れない

先週の月曜日に履き始めたシューズにまだ慣れない。 その翌日に3kmほど走ったら、左足の中指だったかが痛くなった。長い距離を走っていたらマメになるところだった。どうもシューズのつま先の空間が狭いようだ。シューズのサイズは28cmだが、それでも小さい…

第3回鈴鹿山麓かもしかハーフマラソン結果

http://www2.town.komono.mie.jp/icity/browser?ActionCode=content&ContentID=1224044771358&SiteID=000000000000042分台が4人しかいない。ずいぶんばらけている。 同じくらいのタイムで女性が2人いたが、8kmくらいだったかで抜いた女性は28秒差、5kmくらい…

第27回半田市民マラソン大会参加通知書

7:45〜9:55 受付 10:40 スタート ずいぶんゆったりとしたスケジュールだ。 無料シャトルバスがある。詳しくはウェブで。 http://www.city.handa.aichi.jp/sports/shigoto/marathon/traffic.shtml 半田運動公園は陸の孤島にあるらしい。どこの市もひどいね。…

第3回鈴鹿山麓かもしかハーフマラソン(2)

近鉄四日市から湯の山線に乗り換える。単線で、朝は四日市に出る人が多いだろうから、どうしてもこちらが待つ側になる。中菰野駅から歩いて町役場へ。約1kmをだらだらと歩いて。途中で嫌になったので軽くジョギングする。 役場では、はやくも体操をやってい…

第3回鈴鹿山麓かもしかハーフマラソン

でも、走ったのは10km。 前半はベストのペースの22分、あるいはもう少し速く、後半は余裕があればペースアップ、あわよくば42分台を、というおおざっぱな計画。 前のほうにポジションを取って、スタート。スタートラインは数m前なのに、そこまで6秒かかった…

20kmビルドアップ

距離 タイム キロ当たり 5.398 28:33.06 5:17 1.882 10:09.09 5:24 1.882 09:46.55 5:12 1.882 09:19.94 4:58 1.882 09:15.87 4:55 1.882 09:15.81 4:55 5.398 26:30.08 4:55 20.206 1:42:50.40 5:05 前半キロ5分30秒、余裕があれば後半5分という計画。 しか…

多項式の最大公約数

というか、最大公約式? 整係数の範囲でLCDを求めるプログラムを作った。ヒマがあれば解説も後ほど。

月間記録

今日少し走れば140kmを超えるはずだったのに、雨がえらく降っていてちょっと無理。140kmを超えたのは確かまだ1回しかない。こういうことにそれほど意味はないとは思うが。

行列式の計算(6)

行列はリストのリストで表現し、行列式の関数を作った。 def swap(Mat, i): n = len(Mat); for k in range(i + 1, n): if Mat[k][i] != 0: tmp = Mat[i] Mat[i] = Mat[k] Mat[k] = tmp return 1 return 0def det(Mat): result = 1 div = 1 n = len(Mat) for …

行列式の計算(5)

多項式の行列式を計算するために、多項式クラスを作ってみた。 こんな風に使う。 f = Poly([1, 2, 0, 3]) g = Poly([1, -2]) print(f) print(g) print(f + g) print(-g)

行列式の計算(4)

掃き出し法で行列式を求める。これならO(n3)で求められる。その前に、定義を再掲。 まず、 は定義より明らか。 も明らか。 (k行とl行が同じ) これは、σとτ=σ(k l)がペアになっていて、符合が逆でキャンセルする。 これらより、 これにより掃き出し法が可能…

長居公園長距離走路

長居公園内に、1周2.813kmのジョギングコースが整備されている。すなわち、15周すると42.195kmを走ったことになる。ここを走ってみた。 地下鉄長居駅のコインロッカーに荷物を預けて、そこからわずかに歩くと外周コースに出た。1車線よりはだいぶ広い道路の…

行列式の計算(3)

定義の説明が済んだので、行列式の計算がやっとできる。 の行列式を定義に沿って計算する。 S3の要素は、 e, (1 2), (1 3), (2 3), (2 3 1), (3 1 2) 互換の符号は-1、(2 3 1) = (1 3)(1 2)、(3 1 2) = (2 3)(1 2)だから、残りは-1。eは何も変えないから、こ…

薄暮(ナイター)陸上記録会

3ヶ月連続で、今回は5000m。5000をトラックでまじめに走るのは初めてだ。 17時過ぎに受付を済ませて、アリーナのほうで少しくつろいで、それから2kmほど周りを軽く走ったが、汗がやたらと出る。記録会の日はいつも暑いような気がする。 そのあとスタジアムの…

行列式の計算(2)

「全ての置換は互換の合成であらわせる」は、バブルソートを考えれば明らかだろう。 互換の数の偶奇は決まっているかは、証明を忘れた。 ので、ちょっと考えてみた。 もし置換σが偶数と奇数のどちらでも表せたとすると、 ただし、si, tjは互換で、m+nは奇数…

行列式の計算(1)

Pythonで多項式の行列の行列式が計算できるようにする。 行列式の定義 とすると、 対称群 Snはn次対称群で、σはその要素である置換。n次対称群というのは、n個の要素の置換全ての集合である。例えば、対称式 x+y+z があったとき、(x y z) → (y z x)という変…

2008名古屋シティマラソン

今日、地下鉄に広告が出ていた。 早速申し込んでおいた。 申込みに予定タイムが必要。秒単位まで入力させなくても。1時間36分00秒としておいた。これはこれからのトレーニング次第だが。これくらい出しておきたいところ。あまり時間がない。 一色も今日から…

第44回一色マラソン

大会要綱が届いていた。 ハーフのコースが変わった。今年のレースで、最後に失速したあの海岸の吹きさらしの道はコースからなくなった。あそこさえなければ予定通りに走れたのにというあの道が。ここのコースが最悪なのは、順風のときは少し陸地に入ったとこ…