[論文レビュー] Deterministic Random Walks for Rapidly Mixing Chains

Takeharu Shiraga, Yukiko Yamauchi+2 more|arXiv (Cornell University)|Nov 15, 2013

Markov Chains and Monte Carlo Methods参考文献 12被引用数 10

ひとこと要約

本稿では、マーカフ連鎖の決定的アナログである関数ルーターモデルを導入し、無理数を含む任意の遷移確率に対しても一般化されたローター・ルーターモデルを提示する。決定的モデルにおける頂点のトークン数と、それに対応するマーカフ連鎖における期待値との差異について上界を確立し、マッチングサンプリング連鎖ではO(m²n² log n)、線形遷移連鎖ではO(n⁴ log n)の上界を得た。この上界は混合時間および連鎖の構造に依存する。

ABSTRACT

The rotor-router model is a deterministic process analogous to a simple random walk on a graph. This paper is concerned with a generalized model, functional-router model, which imitates a Markov chain possibly containing irrational transition probabilities. We investigate the discrepancy of the number of tokens at a single vertex between the functional-router model and its corresponding Markov chain, and give an upper bound in terms of the mixing time of the Markov chain.

研究の動機と目的

• 関数ルーターモデルを用いて、無理数や任意の遷移確率を持つマーカフ連鎖を一般化し、決定的シミュレーションを可能にする。
• 関数ルーターモデルにおける決定的トークン分布と、それに対応するマーカフ連鎖における期待値との差異を分析する。
• 混合時間およびマーカフ連鎖の構造的性質に依存する、単一頂点における差異の上界を提供する。
• 複雑な問題、たとえばマッチングの数え上げやナップサック解の数え上げにまで、決定的ランダムウォーク技術を拡張する。
• #P完全問題におけるMCMC手法の決定的化のための理論的基盤を確立し、誤差の保証付きで成立する上界を示す。

提案手法

• 固定され事前に計算された遷移シーケンスに基づいてトークンをルーティングすることで、マーカフ連鎖をシミュレートする決定的プロセスとして関数ルーターモデルを導入する。
• 非一様かつ無理数の遷移確率を扱えるように、『関数的』ルーティングの概念を用い、標準的なローター・ルーターモデルを一般化する。
• 混合時間およびスペクトルギャップの概念を用いて、決定的プロセスと確率的プロセスとの差異を制限する。
• 無理数の遷移確率を近似するために、最短残存時間（SRT）ルールの変種およびヴァン・ダール・コルプト列を適用する。
• 遷移行列の2番目に大きな固有値および固有ベクトル構造を用いて境界を導出する。特に可逆連鎖に対して有効である。
• 線形遷移モデル（P_LIN）やマッチングサンプリング連鎖（P_MAT）といった特定の連鎖を分析し、nおよびmを用いた明示的な差異上界を導出する。

実験結果

リサーチクエスチョン

• RQ1 無理数の遷移確率を持つマーカフ連鎖を、差異が有界に保たれる決定的プロセスでシミュレートすることは可能か？
• RQ2 関数ルーターモデルにおける頂点のトークン数と、それに対応するマーカフ連鎖における期待値との最大差異は何か？
• RQ3 マーカフ連鎖の混合時間は、決定的シミュレーションにおける差異にどのように影響するか？
• RQ4 差異上界は、最大次数Δや定常確率の比に依存せずに表現可能か？
• RQ5 このフレームワークは、マッチングの数え上げやナップサック解の数え上げといった#P完全問題に、保証付きの誤差を伴って適用可能か？

主な発見

• 線形遷移モデル（P_LIN）では、差異がO(n⁴ log n)で有界であり、nは状態数である。
• マッチングサンプリング連鎖（P_MAT）では、差異がO(m²n² log n)で有界であり、mは辺の数、nは頂点数である。
• 差異上界はグラフ構造だけでなく、マーカフ連鎖の混合時間およびスペクトル的性質に依存する。
• 上記の上界は可逆かつ再帰的マーカフ連鎖に成立し、関数的ルーティングにより無理数の遷移確率を持つ連鎖へも拡張可能である。
• 結果から、決定的ランダムウォークが#P完全問題におけるMCMC手法の決定的化ツールとして利用可能である可能性が示唆される。
• 解析から、差異が入力サイズに対して多項式的であることが示され、実用的な近似アルゴリズムの実現可能性を支持する。