Hamiltonian Monte Carlo(HMC)について学んだのでメモ。 前提/対象/自身の学習前の知識 参考文献 Metropolis法 Metropolis-Hastings法 Hamiltonian Monte-Carlo 理解の手順 Hamiltonianとは何か導出してみる Lagrangian Legendre変換 Hamiltonian HMCのアル…

のときに がすぐ理解できなかったのでメモ。 手順 あるを使って と書ける。これを変形して ここで明らかにで、または よって

実に5ヶ月ぶり、今年最後のSRM。 8月はSWoPPといろいろ、9月は旅行いったりなんだり、後半から12月頭まではSCとその後処理 としばらく参加できていなかった… 結果 o--(+0/-0) 228.60pts 219th 1959->1965(+6)維持できてよかった。 250 MonstersValley 典型的…

結果 o--(+0/-0) 232.13pts(192nd) 1705->1768(+63) またもやhighest更新です。Petrとはじめて同部屋でした。 250 各bitごとに調べていった。 下から[a,b]のxorを取ったときの下からkbit(0-indexed)目は (([0,b]の下からkbit目が1の個数)-([0,a-1]の下からkb…

ソース http://judge.u-aizu.ac.jp/onlinejudge/review.jsp?rid=404152 方針 制約が大きい…JavaならBigIntegerを使えばよいかもしれない。 2^32より大きいmodを使うとlong longでもoverflowする可能性があるので 10^5((10^5)^2 >= 10^9.6≒2^32)程度のmodを素…

蟻本他を参考にした。 modの逆元 与えられたa,b,mから ax≡b(mod m) を満たすxを求めたいとする。 このとき ay≡1(mod m) なるyが求まれば y*ax≡y*b(mod m) x≡y*b(mod m)(∵ay≡1)としてxを求めることができる。 ay≡1(mod m)⇔∃k(ay=1+km)であるから変形して ay-m…

ソース http://judge.u-aizu.ac.jp/onlinejudge/review.jsp?rid=382200 解き方 継子立てのバリエーション。(ヨセフスの問題とも言うらしい) http://ja.wikipedia.org/wiki/%E3%81%BE%E3%81%BE%E3%81%93%E7%AB%8B%E3%81%A61,2,...,nのn人を円形にclockwiseに…