[論文レビュー] The combinatorics of the colliding bullets problem

Nicolas Broutin, Jean‐François Marckert|arXiv (Cornell University)|Sep 4, 2017

Stochastic processes and statistical mechanics参考文献 13被引用数 3

ひとこと要約

この論文は、[0,1] 内で i.i.d. に一様分布する速度を持つ n 支の弾丸を逐次発砲した際の生存弾丸数の正確な分布を、ランダム置換のサイクル構造と関連付ける独自の組合せ的アプローチにより、有限の衝突弾丸問題を解いている。主な結果は、生存分布に対する再帰的関係式であり、中心化およびスケーリングを行った後、Gumbel 型の極限分布に収束する。その平均および分散は、漸近的にともに ~ (1/2) log n に比例する。

ABSTRACT

The finite colliding bullets problem is the following simple problem: consider a gun, whose barrel remains in a fixed direction; let $(V_i)_{1\le i\le n}$ be an i.i.d.\ family of random variables with uniform distribution on $[0,1]$; shoot $n$ bullets one after another at times $1,2,\dots, n$, where the $i$th bullet has speed $V_i$. When two bullets collide, they both annihilate. We give the distribution of the number of surviving bullets, and in some generalisation of this model. While the distribution is relatively simple (and we found a number of bold claims online), our proof is surprisingly intricate and mixes combinatorial and geometric arguments; we argue that any rigorous argument must very likely be rather elaborate.

研究の動機と目的

• n 支の弾丸が [0,1] 内で i.i.d. に一様分布する速度で発砲される有限の衝突弾丸問題における、生存弾丸数の正確な確率分布を特定すること。
• モデルを、ランダムな発砲間隔、固定速度とランダムな発砲時刻、一般の速度関数に従う加速弾丸を含むように一般化すること。
• 弾丸の衝突プロセスとランダム置換のサイクル構造との間の深い組合せ的関係を確立すること。
• n → ∞ のときの生存弾丸数の極限分布を導出し、適切な中心化およびスケーリングの後、Gumbel 型の分布に収束することを示すこと。

提案手法

• 最初の弾丸の速度とそれに続く部分問題に基づいて、生存分布 q_N(k) に対する再帰的関係式を定義する。
• メモリが 2 のマルコフ連鎖モデルを用い、状態が直前の2つの生存数に依存する。
• 弾丸衝突の結果とランダム置換のサイクル分解の間の双対写像を確立する。特に、奇数長のサイクルに注目する。
• 2変数の母関数と特異点解析を用いて、分布の漸近的挙動を導出する。
• 解析的組合せ論における準パワー定理を適用し、中心化およびスケーリング後のガウス型極限分布の存在を証明する。
• 生存弾丸数が、サイズ n の一様ランダム置換における奇数長サイクルの数と同一分布に従うという事実を活用する。

実験結果

リサーチクエスチョン

• RQ1 n 支の弾丸が [0,1] 内で i.i.d. に一様分布する速度で発砲された際、生存弾丸数の正確な分布は何か？
• RQ2 発砲間隔がランダムまたは固定である場合、生存弾丸の分布はどのように変化するか？
• RQ3 n → ∞ のとき、生存弾丸数の極限分布は何か？
• RQ4 生存プロセスに置換のサイクル構造に基づく組合せ的解釈は存在するか？
• RQ5 生存数の漸近的挙動は、平均および分散の増大度として特徴付けられるか？

主な発見

• 生存分布は再帰式 q_N(k) = (1/N) q_{N-1}(k-1) + (1 - 1/N) q_{N-2}(k) を満たし、初期条件は q_1(1)=1, q_1(0)=0, q_0(0)=1 である。
• 生存弾丸数は、サイズ n の一様ランダム置換における奇数長サイクルの数と同一分布に従う。
• 生存弾丸の期待値は漸近的に (1/2) log n に比例し、分散についても同様に ~ (1/2) log n に比例する。
• 中心化を (1/2) log n で、スケーリングを √(log n) で行った後、生存弾丸の分布は Gumbel 型の極限分布に収束する。
• 極限分布は特異点解析および準パワー定理を用いて導出され、対数スケールにおけるガウス型フラクチュエーションが確認された。
• 証明は、サイクルへの挿入と2段階の有向木構造を含む非自明な組合せ的分解に依存しており、これが再帰的構造を生み出している。