習っておくと筋が良くなると言われるschemeを勉強しようと思い立ち、SICPを読み始めた。無料で全文公開だし、演習問題も豊富についてるし、非常に有難い。 wgetで全体を落としてスマホの中に突っ込んどけば、どこでも読める。emacsの上で、emacs-w3mを使って…

プレゼントに対する反応 長男は、現在小学3年生である。もうそろそろ、現実的になっても良い年頃なのではないかと思うのだが、彼は未だにサンタクロースを信じている。とはいうものの、さすがに多少の疑念が芽生え始めたようだ。例年ならば、「今年はサンタ…

10月から金利デリバティブチームに移った。そのために、Libor Market Modelの復習をしようと思った。そのために、測度変換について復習をした。その時のメモ。 Forward Measure

問題 http://code.google.com/codejam/contest/dashboard?c=889487#s=p2 方針

問題 https://code.google.com/codejam/contest/dashboard?c=889487#s=p1 方針

単体テストというものの有難味が分かり始めた。Rでの単体テストフレームワークであるRUnitの使い方をメモる。 テスト駆動開発なので、まずはテストを書く。 test.func1 <- function(){ checkEquals(func1(), 1, "func1 : NG") } test.func2 <- function(){ c…

結果は同じだけど、実行時間は大分違う。 x <- matrix(0, nrow=100, ncol=100) system.time(for(j in 1:10000)for(i in 1:100)x[0,i] <- 1) ## user system elapsed ## 2.203 0.037 2.321 system.time(for(j in 1:10000)x[0,] <- 1) ## user system elapsed …

という関数を考える。1回微分したらヘビサイド関数、2回微分したらデルタ関数になりそうだ。で、実際、 となるそうだ。

連続性と絶対連続性 区間上で定義された関数に関して、普通の連続性よりも強い連続性の概念に絶対連続性というのがある(cf. wikipedia)。が絶対連続であるのは、関数が存在して がで成り立つのと同値らしい。この時、はa.e.で微分出来て、になる。 連続だけ…

その1 の続き。

恐竜博に行ってきた。子供が、学校でチラシを貰って来て、行ってみたいと言うので。本物の恐竜の化石って、子供の頃にやっぱりどこかの特別展で見たような記憶があるくらい。子供だったので、自分の中に特別な印象は、殆ど無い。しかし、大人になって見ると…

Rのリストでは、要素に名前が付いているとき、list$nameでその要素にアクセス出来る。それを使っていて、たまたま変な挙動を見つけたので、メモ。 > l <- list(x=1, xxx=2, y=3, yyy=4) > l$x [1] 1 > l$xx [1] 2 > l$y [1] 3 > l$yy [1] 4 > l[["xx"]] NULL…

メモ1 2本のベクトルx, yがある時に、x, y両方に含まれている要素を調べる方法。 > x <- 2 * (1:10) > y <- 3 * (1:10) > x[apply(outer(x, y, '=='), 1, any)] [1] 6 12 18 メモ2 Rの付置演算子をemacsのalignで揃える。 ;; align for R (require 'align) (…

解析教程 上 新装版作者: E.ハイラー,G.ワナー,蟹江幸博出版社/メーカー: シュプリンガー・ジャパン(株)発売日: 2006/09/15メディア: 単行本 クリック: 22回この商品を含むブログ (5件) を見る解析教程を読んでいて、何となく開集合とか位相とか思い出した。…

方針 DP。というかメモ化再帰。 f(x)を長さxの時に必要な最小手数とすると、 f(x) = min{min{f(x/p) + 1| p : xの素因数}, min(m-1, x-m+1)} が成り立つ。但し、m = M % x。 この漸化式の気持ちは、 長さNの時の最小手数は、 1回 Divide 操作をやってさらに …

方針 素直に調べる。 雑感 何で Round 1 なの？

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11096&rd=14514 下1桁が 4 or 7 の正の自然数を lucky number と呼ぶ 正の自然数 n を lucky number の和で表したい。この時、なるべく summand の数を少なくしたい。いくつの数が必要か？ 方針 f…

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11432&rd=14514

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11279&rd=14435 (要ログイン) 方針 やるだけ。

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11154&rd=14435 (要ログイン) 初期状態の都市のサイズが population[i] (i = 0, .., N-1) で与えられる。 任意の順番で、都市を2つづつ選び、それらを併合して1つの都市にする。最終的に1つだけが…

問題 http://code.google.com/codejam/contest/dashboard?c=975485#s=p2 方針 これは、方針を考えるのに苦労した。けど、方針が分かればコーディングは一瞬。前の2問とは逆のパターン。 求められていることを、抽象的に述べると、次のとおりになる。 数列 a[…

問題 http://code.google.com/codejam/contest/dashboard?c=975485#s=p1 方針 素直にシミュレーション

問題 http://code.google.com/codejam/contest/dashboard?c=975485#s=p0 2本の廊下にそれぞれBlue, Orangeのロボットが居る。 各廊下は長さ100で、地点 1, 2, ..., 100にボタンがある。 O 2, B 1, B 2, O 4 というようなリストでボタンを押す順番のリストが…

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11399&rd=14434 (要ログイン) 全部小文字で書かれた英語の文章があるので、文の先頭を大文字に変えろ。 方針 問題の条件をよく読むと、「文の先頭 = 文章の先頭 or "."の2文字後」、なので、文の…

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11397&rd=14434 (要ログイン) 数列 a[1], a[2], ..., a[n] が以下の2条件を満たす時、その数列は perfect であると言う。 a[1] + ... + a[n] = a[1] * ... * a[n] a[1] >= 0, ..., a[n] >= 0 適当…

キャンプ場 5/3,4と河口湖畔でキャンプをしてきた。家族でキャンプをするのは、去年のGWに引き続き今回で2回目。今回は、富士山の麓、戸沢センターキャンプ場にお世話になった。富士山の麓の河口湖の湖畔にテントを張って、一泊。

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11233&rd=14433 黒と白の玉が一列に並んでいる。 先頭から順番に取り出す。 取り出した玉が 白なら、残った玉の順番をひっくり返す 黒なら、玉の色を反転させる。 方針 基本的には、素直に調べる…

問題 http://www.topcoder.com/stat?c=problem_statement&pm=11350&rd=14433 (要ログイン) 方針 素直に調べる。

問題 http://code.google.com/codejam/contest/dashboard?c=619102#s=p0 ビルが2つ並んでいて、左のビルの高さA[i]の所から右のビルの高さB[i]の所に電線が張ってある。 電線の交点の数を求めなさい、という問題。 但し、3つの線が一点で交わる事はない事は…

問題 http://code.google.com/codejam/contest/dashboard?c=544101# 縦横Nマスに文字R, Bが配置されている。 これを右に90度回転させて、各文字を「重力」で下に落とす。 この時、R, BがK個以上縦横斜めのいずれかの方向に並んでいるかどうかをチェックしな…