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

Mebius関数を素因数分解しないで求めるメモ

MANUEL BENITO AND JUAN L. VARONA. RECURSIVE FORMULAS RELATED TO THE SUMMATION OF THE MOBIUS FUNCTIONにあったやつ Mebius関数とは (nがある素数の2乗で割り切れる) (nが相異なるk個の素因数に分解できる) Wikipedia メビウス関数 計算方法 Mertens fun…

Rの課題?

メモ http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1377312633課題5 対角要素が1:10の対角行列を作りなさい。答え diag(1:10,ncol=10,nrow=10);http://www.okada.jp.org/RWiki/?%B9%D4%CE%F3Tips%C2%E7%C1%B4#n42877c2課題6 matrix(1:100,10,1…

正多角形グラフの隣接行列とハミルトン閉路,パスの数を数えてみる

MathWorldに載っているけれど確かめてみた http://ideone.com/WayQn結果だけ 正四面体 // Hamiltonian Path: 6本 // Cycle: 6本 // Tetrahedral graph */正六面体 // Hamiltonian Path: 18本 // Cycle: 12本 // Cubical graph */正八面体 // Hamiltonian Pat…

Rで素数の無限リスト

以下を参考に.. SICP 3.5 Diaspar Journal 関数型(12) - 無限リスト # 無限リスト[x, fn(x), fn(fn(x)), ... ]を作る infinite.stream <- function(x, fn){ list(car=x, cdr=function() infinite.stream(fn(x),fn)) } # 無限リストlstから最初のn項のベク…

リンゴ振り分け問題をGLPKで解く

3つの八百屋から値段の異なるリンゴをいくつか仕入れた。 これを5つの袋に決められた個数ずつ振り分ける。 袋ごとの平均単価が同じくらいになるように振り分ける方法を見つけたい。元ネタはここ要するに画像の5x3の表を埋めつつ目的関数を最小化したい ただ…

無羨望分割(envy-free division) ケーキを公平に切り分ける方法のお話あとで読む http://www.brl.ntt.co.jp/people/manabe/cutMFCS.pdf