ランダムウォークの原点への復帰

同義語:1次元のランダムウォーク(高校数学)return of a random walk to the origin

概要

ランダムウォークの原点への復帰(return of a random walk to the origin)とは、原点から出発し、確率 $\frac12$ ずつで独立に $+1$ か $-1$ 動く数直線上の点(1次元のランダムウォーク)が、再び原点に戻る現象である。$2n$ 回目に原点にいる確率は $u_{2n}=\binom{2n}{n}\frac{1}{4^n}$ で、$2n$ 回目まで一度も原点に戻らない確率も $u_{2n}$ に等しい。$\frac{1}{2\sqrt n}\le u_{2n}\le\frac{1}{\sqrt{2n+1}}$ から、$2n$ 回目までに原点に戻る確率 $1-u_{2n}$ は $n\to\infty$ で $1$ に近づく。$2n$ 回目までに戻る回数の期待値は $(2n+1)u_{2n}-1$ で、限りなく大きくなる。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 反復試行の確率, 最短経路の数え上げと鏡像原理, 確率変数の期待値と分散

高校での出発点:硬貨を投げて数直線を歩く

数直線の原点に点を置き、硬貨を投げる。表が出たら右へ $1$、裏が出たら左へ $1$ だけ点を動かす。これを何回もくり返すと、点は右へ行ったり左へ行ったりしながら、ときどき原点に戻ってくる。図 1 は、$20$ 回ずつ投げた 3 回の試しの記録で、横に回数、縦に点の位置をとり、原点に戻った回に丸を付けた。
硬貨を 20 回投げて、表なら +1、裏なら −1 動いた点の位置の記録を 3 回分並べた図。丸は原点に戻った回で、1 本目は 5 回、2 本目は 3 回、3 本目は 5 回戻っている 硬貨を 20 回投げて、表なら +1、裏なら −1 動いた点の位置の記録を 3 回分並べた図。丸は原点に戻った回で、1 本目は 5 回、2 本目は 3 回、3 本目は 5 回戻っている
2 本目の記録では、$13$ 回目まで一度も原点に戻っていない。表と裏は同じ確率で出るのに、点が片側に長くとどまることがある。まず、硬貨を 4 回投げる場合を全部数えてみる。

4 回投げる場合を全部数える

表を $+$、裏を $-$ と書く。4 回の結果の並びは $2^4=16$ 通りで、どれも確率 $\dfrac1{16}$ で起こる。
(1) 4 回目に原点にいるのは、$+$ と $-$ がちょうど 2 個ずつのときで、並び方は $\dbinom42=6$ 通りある:
$$ {+}{+}{-}{-},\quad {+}{-}{+}{-},\quad {+}{-}{-}{+},\quad {-}{+}{+}{-},\quad {-}{+}{-}{+},\quad {-}{-}{+}{+} $$
よって、4 回目に原点にいる確率は $\dfrac6{16}=\dfrac38$ である。
(2) 1 回目から 4 回目まで、一度も原点に戻らない並びを数える。1 回目が $+$ のとき、位置は $1$ で、2 回目に $-$ が出ると原点に戻るので 2 回目は $+$ で位置 $2$ である。3 回目は $+$(位置 $3$)でも $-$(位置 $1$)でもよい。位置 $3$ からは 4 回目に何が出ても原点に戻らず、位置 $1$ からは 4 回目に $+$ でなければならない。よって
$$ {+}{+}{+}{+},\quad {+}{+}{+}{-},\quad {+}{+}{-}{+} $$
の 3 通りである。1 回目が $-$ のときも、符号を入れ替えた 3 通りがある。合わせて 6 通りで、確率は $\dfrac6{16}=\dfrac38$ である。
(1) と (2) の確率は、どちらも $\dfrac38$ で等しい(図 2)。

上の段は 4 回目に原点にいる 6 本、下の段は 4 回目まで一度も原点に戻らない 6 本の記録。どちらも 6 本ずつある 上の段は 4 回目に原点にいる 6 本、下の段は 4 回目まで一度も原点に戻らない 6 本の記録。どちらも 6 本ずつある
「4 回目に原点にいる」と「4 回目まで一度も原点に戻らない」は、まったく違う出来事である。それなのに確率が同じになった。この記事では、これが偶然ではなく、何回投げても成り立つことを示し、それを使って次の問いに答える。

  1. $2n$ 回目に原点にいる確率はいくつか。→ thm-rwr-main の (1)
  2. $2n$ 回目まで一度も原点に戻らない確率はいくつか。→ thm-rwr-main の (2)
  3. 投げ続ければ、点はいつか必ず原点に戻るか。→ thm-rwr-eventual
  4. $2n$ 回までに原点に戻る回数の期待値はいくつか。→ cor-rwr-count
  5. 硬貨の確率や 1 回の動き方を変えると、何が変わるか。→ ex-rwr-cx-biased、ex-rwr-cx-lazy、ex-rwr-cx-odd
    高校の計算この記事の言葉大学の言葉
    表と裏の並べ方を数える長さ $n$ の歩み確率変数の列 $X_1,X_2,\ldots$
    碁盤の目の最短経路歩みを格子の道とみるランダムウォークの経路
    鏡像原理・投票定理原点に戻らない歩みの数鏡像原理
    中央の二項係数 $\binom{2n}n$ の大きさ$u_{2n}$ の大きさStirling の公式による漸近

言葉の準備:歩みと原点への復帰

1 次元のランダムウォーク

次の仮定のもとで、数直線の上を動く点を考える。
(R1) $k=1,2,3,\ldots$ について、$k$ 回目の動き $X_k$ は $+1$ か $-1$ で、それぞれ確率 $\dfrac12$ をとる。$X_1,X_2,\ldots$ は互いに独立である(硬貨を独立に投げる)。
(R2) 点ははじめ原点にあり、$n$ 回動いた後の位置を $S_n=X_1+X_2+\cdots+X_n$ とする($S_0=0$)。
この点の動きを 1 次元のランダムウォーク という。$n$ 回目までの動きの並び $(X_1,\ldots,X_n)$ を 長さ $n$ の歩み という。長さ $n$ の歩みは $2^n$ 通りあり、条件 (i) から、どれも確率 $\left(\dfrac12\right)^n$ で起こる。$k\ge1$ について $S_k=0$ となるとき、点は $k$ 回目に原点に戻る という(出発の $S_0=0$ は戻ったと数えない)。

確率がすべて等しいので、$n$ 回目までの出来事の確率は「その出来事が起こる長さ $n$ の歩みの数」を $2^n$ で割れば求まる。この記事の計算はすべてこの形で、反復試行の確率 と同じく、並び方を数える問題になる。

位置の計算と、原点に戻る回

長さ $6$ の歩み $(+1,-1,-1,+1,+1,+1)$ では、位置は
$$ S_1=1,\quad S_2=0,\quad S_3=-1,\quad S_4=0,\quad S_5=1,\quad S_6=2 $$
と動き、2 回目と 4 回目に原点に戻る。長さ $6$ の歩み $(+1,+1,-1,+1,-1,-1)$ では、位置は $1,2,1,2,1,0$ と動き、6 回目に初めて原点に戻る。

原点に戻れるのは偶数回目だけである。

奇数回目には原点にいない

$n$ 回動いた後、$+1$ の動きが $a$ 回、$-1$ の動きが $b$ 回あったとすると、$S_n=a-b$、$n=a+b$ である。特に、$S_n=0$ となるのは $n$ が偶数で $a=b=\dfrac n2$ のときに限る。したがって、奇数回目に原点にいる確率は $0$ である。

$S_n$ は $+1$ を $a$ 個、$-1$ を $b$ 個足したものなので $S_n=a-b$ であり、動きの個数を数えて $n=a+b$ である。$S_n=0$ なら $a=b$ で、$n=a+b=2a$ は偶数、$a=b=\dfrac n2$ である。$n$ が奇数なら $a=b$ となることはないので、$S_n=0$ となる歩みは 1 つもなく、確率は $0$ である。$\square$

そこで、偶数回目に原点にいる確率に名前を付ける。

$u_{2n}$

$n=1,2,3,\ldots$ について、$2n$ 回目に原点にいる確率を
$$ u_{2n}=P(S_{2n}=0) $$
と書く。また $u_0=1$ と決める(出発の時点では必ず原点にいる)。

歩みを格子の道とみる

歩みは、最短経路の数え上げと鏡像原理 の格子の道で表せる。$+1$ の動きを右へ 1(R)、$-1$ の動きを上へ 1(U)に対応させると、長さ $n$ の歩みは、原点から出発する R と U の $n$ 歩の道になる。$k$ 歩目の後の点を $(a_k,b_k)$ とすると、$a_k$ は $+1$ の回数、$b_k$ は $-1$ の回数なので
$$ S_k=a_k-b_k $$
である。したがって

  • $S_k=0$(原点にいる)は、点 $(a_k,b_k)$ が直線 $y=x$ の上にあること、
  • $S_k>0$ は、点 $(a_k,b_k)$ が直線 $y=x$ の下側($y< x$)にあること
    と同じである。「原点に戻らない」ことは「直線 $y=x$ に触れない」ことに言いかえられ、最短経路の数え上げと鏡像原理 の鏡像原理と投票定理が使える。

主定理 1:原点にいる確率と、原点に戻らない確率

ex-rwr-start では、4 回目まで一度も原点に戻らない確率と、4 回目に原点にいる確率が等しかった。6 回でも数えてみる。

$n=3$ で、正の側にとどまる歩みを数える

長さ $6$ の歩みで、$S_1,\ldots,S_6$ がすべて正のものを、最後の位置 $S_6$ で分けて数える。$S_6$ は偶数で(prop-rwr-parity と同じ理由で $S_6=a-b$、$a+b=6$)、正なので $2,4,6$ のどれかである。
(1) $S_6=6$:$+1$ が 6 回の 1 通り。
(2) $S_6=4$:$+1$ が 5 回、$-1$ が 1 回。$-1$ を 1 回目に置くと $S_1=-1$、2 回目に置くと $S_2=0$ になるので、$-1$ は 3〜6 回目のどこかで、4 通り。この 4 通りでは、$-1$ の直前の位置が $2$ 以上なので、位置は $1$ 以上のままである。
(3) $S_6=2$:$+1$ が 4 回、$-1$ が 2 回の並び $\dbinom62=15$ 通りのうち、位置がずっと正のものは
$$ {+}{+}{+}{+}{-}{-},\ \ {+}{+}{+}{-}{+}{-},\ \ {+}{+}{+}{-}{-}{+},\ \ {+}{+}{-}{+}{+}{-},\ \ {+}{+}{-}{+}{-}{+} $$
の 5 通りである。
合わせて $1+4+5=10$ 通りある。すべて負のものも、符号を入れ替えて $10$ 通りあるので、6 回目まで一度も原点に戻らない歩みは $20$ 通りで、$\dbinom63=20$ と一致する。確率は $\dfrac{20}{64}=\dfrac5{16}$ で、$u_6=\dbinom63\dfrac1{4^3}=\dfrac{20}{64}=\dfrac5{16}$ に等しい。

ex-rwr-n3 の (3) の $5$ 通りは、最短経路の数え上げと鏡像原理 の投票定理の場面そのものである。$+1$ を候補者 A の票、$-1$ を候補者 B の票とみると、「位置がずっと正」は「開票の間、A の票数がつねに B より多い」ことになる。投票定理によれば、A が $a$ 票、B が $b$ 票($a>b$)のとき、その確率は $\dfrac{a-b}{a+b}$ である。$a=4$、$b=2$ なら $\dfrac26\cdot\dbinom62=\dfrac13\cdot15=5$ 通りで、数え上げと一致する。一般の場合の証明は、この見方で最後の位置ごとに数えて足す。

原点にいる確率と、原点に戻らない確率

$n$ を正の整数とする。
(1) $2n$ 回目に原点にいる確率は
$$ u_{2n}=\binom{2n}n\frac1{4^n} $$
である。
(2) 長さ $2n$ の歩みのうち、$S_1,S_2,\ldots,S_{2n}$ がどれも $0$ でないもの($2n$ 回目まで一度も原点に戻らないもの)の数は $\dbinom{2n}n$ である。したがって
$$ P(S_1\ne0,\ S_2\ne0,\ \ldots,\ S_{2n}\ne0)=u_{2n} $$
である。

方針:(1) は並べ方を数えるだけである。(2) は、まず「原点に戻らない歩みは、ずっと正か、ずっと負か」に分ける。ずっと正の歩みを最後の位置 $S_{2n}=2j$ ごとに投票定理で数え、その和が望遠鏡のように消えて $\dfrac12\dbinom{2n}n$ になることを示す。
段 1((1) の証明)。prop-rwr-parity により、$S_{2n}=0$ となるのは、$2n$ 回の動きのうち $+1$ がちょうど $n$ 回、$-1$ が $n$ 回のときである。$+1$ を置く $n$ か所の選び方は $\dbinom{2n}n$ 通りある(格子の道で言えば、$(0,0)$ から $(n,n)$ への最短経路の数である)。長さ $2n$ の歩みは $2^{2n}=4^n$ 通りで、どれも同じ確率なので、$u_{2n}=\dbinom{2n}n\Big/4^n$ である。
段 2(ずっと正か、ずっと負か)。$S_1,\ldots,S_{2n}$ がどれも $0$ でない歩みを考える。$S_1$ は $+1$ か $-1$ である。$S_1=1$ なのに、ある $k$ で $S_k<0$ となったとする。位置は整数で、1 回に $1$ ずつしか変わらないので、$1$ 以上の値から $-1$ 以下の値に移る途中で、必ず $0$ を通る。これは「どれも $0$ でない」に反する。よって $S_1=1$ なら $S_1,\ldots,S_{2n}$ はすべて正である。同じように、$S_1=-1$ ならすべて負である。
また、歩みのすべての動きの符号を入れ替える($X_k$ を $-X_k$ にする)と、位置もすべて符号が入れ替わる。この入れ替えは、ずっと正の歩みとずっと負の歩みを 1 対 1 に対応させる(2 回入れ替えると元に戻る)。したがって
$$ (\text{原点に戻らない歩みの数})=2\times(\text{ずっと正の歩みの数}) $$
である。
段 3(ずっと正の歩みを、最後の位置ごとに数える)。ずっと正の歩みでは $S_{2n}>0$ で、prop-rwr-parity と同じく $S_{2n}=a-b$、$a+b=2n$ なので $S_{2n}$ は偶数である。$S_{2n}=2j$($j=1,2,\ldots,n$)とおくと、$+1$ の回数は $a=n+j$、$-1$ の回数は $b=n-j$ である。$+1$ を A の票、$-1$ を B の票とみると、「$S_1,\ldots,S_{2n}$ がすべて正」は「最初の 1 票から最後まで、A の票数がつねに B の票数より多い」ことと同じである。$a>b\ge0$ なので、最短経路の数え上げと鏡像原理 の定理「投票定理」が使えて、A と B の並べ方 $\dbinom{2n}{n+j}$ 通りのうち、その割合は $\dfrac{a-b}{a+b}=\dfrac{2j}{2n}=\dfrac jn$ である。よって、$S_{2n}=2j$ となるずっと正の歩みの数は
$$ \frac jn\binom{2n}{n+j} $$
である。
段 4(差の形に書きかえる)。$m=n+j$ とおく($n+1\le m\le2n$)。階乗で書くと
$$ \binom{2n-1}{m-1}=\frac{(2n-1)!}{(m-1)!\,(2n-m)!}=\frac m{2n}\cdot\frac{(2n)!}{m!\,(2n-m)!}=\frac m{2n}\binom{2n}m $$
である(分母と分子に $2n$ と $m$ を掛けて、$(2n)!=2n\cdot(2n-1)!$、$m!=m\cdot(m-1)!$ を使った)。同じように、$m\le2n-1$ なら
$$ \binom{2n-1}{m}=\frac{(2n-1)!}{m!\,(2n-1-m)!}=\frac{2n-m}{2n}\cdot\frac{(2n)!}{m!\,(2n-m)!}=\frac{2n-m}{2n}\binom{2n}m $$
であり、$m=2n$ のときは両辺とも $0$ である($\dbinom{2n-1}{2n}=0$ と約束する)。2 つの式を引くと
$$ \binom{2n-1}{m-1}-\binom{2n-1}{m}=\frac{m-(2n-m)}{2n}\binom{2n}m=\frac{2j}{2n}\binom{2n}{n+j}=\frac jn\binom{2n}{n+j} $$
となる($m-(2n-m)=2m-2n=2j$)。つまり、段 3 の数は
$$ \frac jn\binom{2n}{n+j}=\binom{2n-1}{n+j-1}-\binom{2n-1}{n+j} $$
と差の形に書ける。
段 5(足すと消える)。段 4 の式を $j=1,2,\ldots,n$ について足すと、
$$ \sum_{j=1}^n\left\{\binom{2n-1}{n+j-1}-\binom{2n-1}{n+j}\right\}=\left\{\binom{2n-1}{n}-\binom{2n-1}{n+1}\right\}+\left\{\binom{2n-1}{n+1}-\binom{2n-1}{n+2}\right\}+\cdots+\left\{\binom{2n-1}{2n-1}-\binom{2n-1}{2n}\right\} $$
となり、隣り合う中かっこの間で項が打ち消し合って、最初の項と最後の項だけが残る:
$$ \binom{2n-1}{n}-\binom{2n-1}{2n}=\binom{2n-1}{n} $$
段 4 の 2 つ目の式で $m=n$ とおくと $\dbinom{2n-1}n=\dfrac{2n-n}{2n}\dbinom{2n}n=\dfrac12\dbinom{2n}n$ なので、ずっと正の歩みの数は $\dfrac12\dbinom{2n}n$ である。
段 6(結論)。段 2 より、原点に戻らない歩みの数は $2\times\dfrac12\dbinom{2n}n=\dbinom{2n}n$ である。これを $4^n$ で割ると、確率は $\dbinom{2n}n\Big/4^n=u_{2n}$ で、段 1 の値に等しい。$\square$

  1. は「$2n$ 回目まで一度も原点に戻らない確率は、$2n$ 回目に原点にいる確率に等しい」と言っている。2 つの出来事は違うのに、確率は何回投げても等しい。この等式は、GS06 の第 12.1 節の練習問題 2 (d) と練習問題 7 で、読者が示す問題として挙げられている。この記事の後半と、初めて原点に戻る時刻・最後に原点にいた時刻 は、すべてこの等式を出発点にする。
段 3〜段 5 を $n=3$ で追う

$n=3$ では、段 3 の数は $j=1,2,3$ について
$$ \frac13\binom64=\frac{15}3=5,\qquad\frac23\binom65=\frac{2\cdot6}3=4,\qquad\frac33\binom66=1 $$
で、ex-rwr-n3 の (3)・(2)・(1) の数と一致する。段 4 の差の形では
$$ 5=\binom53-\binom54=10-5,\qquad4=\binom54-\binom55=5-1,\qquad1=\binom55-\binom56=1-0 $$
で、足すと $10-5+5-1+1-0=10=\binom53=\dfrac12\binom63$ となる。

$u_{2n}$ の値は、次の漸化式で順に計算できる。

$u_{2n}$ の漸化式と積の形

$n\ge1$ について
$$ u_{2n}=\frac{2n-1}{2n}\,u_{2n-2} $$
である。したがって
$$ u_{2n}=\frac12\cdot\frac34\cdot\frac56\cdots\frac{2n-1}{2n} $$
である。

thm-rwr-main の (1) から
$$ \frac{u_{2n}}{u_{2n-2}}=\frac{\binom{2n}n}{\binom{2n-2}{n-1}}\cdot\frac{4^{n-1}}{4^n}=\frac{(2n)!\,\{(n-1)!\}^2}{(n!)^2\,(2n-2)!}\cdot\frac14=\frac{2n(2n-1)}{n^2}\cdot\frac14=\frac{2n-1}{2n} $$
である($(2n)!=2n(2n-1)\cdot(2n-2)!$、$n!=n\cdot(n-1)!$ で約分した)。$n=1$ でも $u_2=\dfrac12=\dfrac12u_0$ で成り立つ。この式を $n,n-1,\ldots,1$ と順に使い、$u_0=1$ とすると積の形を得る。$\square$

$u_{2n}$ の値
$n$$1$$2$$3$$4$$5$
$\binom{2n}n$$2$$6$$20$$70$$252$
$4^n$$4$$16$$64$$256$$1024$
$u_{2n}$$\frac12$$\frac38$$\frac5{16}$$\frac{35}{128}$$\frac{63}{256}$
小数$0.500$$0.375$$0.313$$0.273$$0.246$

漸化式でも、$u_4=\dfrac34\cdot\dfrac12=\dfrac38$、$u_6=\dfrac56\cdot\dfrac38=\dfrac5{16}$、$u_8=\dfrac78\cdot\dfrac5{16}=\dfrac{35}{128}$ と確かめられる。$u_{2n}$ は $n$ とともに減っていくが、減り方はゆっくりで、$u_{100}$($n=50$)でも約 $0.0796$ ある。

主定理 2:いつかは原点に戻る

thm-rwr-main の (2) により、「$2n$ 回目までに少なくとも 1 回原点に戻る」確率は、余事象の確率を $1$ から引いて $1-u_{2n}$ である。$n=2$ なら $1-\dfrac38=\dfrac58$、$n=5$ なら $1-\dfrac{63}{256}=\dfrac{193}{256}\approx0.754$ である。$n$ を大きくするとこの確率が $1$ に近づくかどうかは、$u_{2n}$ が $0$ に近づくかどうかで決まる。それを、prop-rwr-recur の積の形を上下から評価して確かめる。

$u_{2n}$ の上下の評価

$n\ge1$ について
$$ \frac1{2\sqrt n}\le u_{2n}\le\frac1{\sqrt{2n+1}} $$
である。

方針:積 $u_{2n}=\dfrac12\cdot\dfrac34\cdots\dfrac{2n-1}{2n}$ の 2 乗を考え、各因子の 2 乗を、掛けると約分で消えていく分数で上下からはさむ。
段 1(上からの評価)。$k\ge1$ について、$(2k-1)(2k+1)=4k^2-1<4k^2$ なので、両辺を正の数 $2k(2k+1)$ で割ると
$$ \frac{2k-1}{2k}<\frac{2k}{2k+1} $$
である。両辺に正の数 $\dfrac{2k-1}{2k}$ を掛けると
$$ \left(\frac{2k-1}{2k}\right)^2<\frac{2k-1}{2k}\cdot\frac{2k}{2k+1}=\frac{2k-1}{2k+1} $$
である。これを $k=1,2,\ldots,n$ について掛け合わせると、右辺の積は
$$ \frac13\cdot\frac35\cdot\frac57\cdots\frac{2n-1}{2n+1}=\frac1{2n+1} $$
と、隣り合う分数の分子と分母が約分で消える。よって $u_{2n}^2<\dfrac1{2n+1}$ で、両辺の正の平方根をとって $u_{2n}<\dfrac1{\sqrt{2n+1}}$ である。
段 2(下からの評価)。$k\ge2$ について、$(2k-1)^2=4k^2-4k+1>4k^2-4k=2k(2k-2)$ なので、両辺を正の数 $2k(2k-1)$ で割ると
$$ \frac{2k-1}{2k}>\frac{2k-2}{2k-1} $$
である。両辺に $\dfrac{2k-1}{2k}$ を掛けると
$$ \left(\frac{2k-1}{2k}\right)^2>\frac{2k-1}{2k}\cdot\frac{2k-2}{2k-1}=\frac{2k-2}{2k}=\frac{k-1}k $$
である。$k=1$ の因子は $\left(\dfrac12\right)^2=\dfrac14$ のままにして、$k=2,\ldots,n$ の不等式を掛け合わせると
$$ u_{2n}^2\ge\frac14\cdot\frac12\cdot\frac23\cdots\frac{n-1}n=\frac14\cdot\frac1n=\frac1{4n} $$
である($n=1$ では右辺の積は $\dfrac14$ だけで、等号が成り立つ)。正の平方根をとって $u_{2n}\ge\dfrac1{2\sqrt n}$ である。$\square$

評価を数で確かめる
  1. $n=1$:$\dfrac1{2\sqrt1}=0.5$、$u_2=0.5$、$\dfrac1{\sqrt3}\approx0.577$。下の評価で等号が成り立つ。
  2. $n=5$:$\dfrac1{2\sqrt5}\approx0.224$、$u_{10}=\dfrac{63}{256}\approx0.246$、$\dfrac1{\sqrt{11}}\approx0.302$。
  3. $n=50$:$\dfrac1{2\sqrt{50}}\approx0.0707$、$u_{100}\approx0.0796$、$\dfrac1{\sqrt{101}}\approx0.0995$。
いつかは原点に戻る

$n\ge1$ について、$2n$ 回目までに少なくとも 1 回原点に戻る確率は $1-u_{2n}$ であり、
$$ \lim_{n\to\infty}(1-u_{2n})=1 $$
である。

方針:前半は thm-rwr-main の (2) の余事象である。後半は lem-rwr-bounds の上の評価とはさみうちの原理による。
段 1(前半)。「$2n$ 回目までに少なくとも 1 回原点に戻る」の余事象は「$S_1,\ldots,S_{2n}$ がどれも $0$ でない」である。thm-rwr-main の (2) より余事象の確率は $u_{2n}$ なので、求める確率は $1-u_{2n}$ である。
段 2(後半)。lem-rwr-bounds より $0< u_{2n}\le\dfrac1{\sqrt{2n+1}}$ である。$n\to\infty$ で $\dfrac1{\sqrt{2n+1}}\to0$ なので、はさみうちの原理により $u_{2n}\to0$ であり、$1-u_{2n}\to1$ である。$\square$

「いつかは原点に戻る確率」は、無限に続く歩みについての確率で、長さが有限の歩みを数えるこの記事の方法では直接は扱えない。そこで、GS06 の Example 12.1 と同じく、「$2n$ 回目までに戻る確率」の $n\to\infty$ での極限を「いつかは戻る確率」と定める。thm-rwr-eventual は、この意味で 原点に戻る確率は $1$ であることを示している。
$1-u_{2n}$ は $1$ に近づくが、その近づき方は遅い(図 3)。ex-rwr-bounds の $n=50$ でも、$100$ 回投げて一度も原点に戻らない確率が約 $8\%$ ある。
2n 回目までに原点に戻る確率 1 − u_2n の n による変化。1 に近づくが遅い。帯は上下の評価からわかる範囲 2n 回目までに原点に戻る確率 1 − u_2n の n による変化。1 に近づくが遅い。帯は上下の評価からわかる範囲
$u_{2n}$ の大きさは、もっと正確に言える。Wallisの公式(高校数学) の系「Wallis の積と中央の二項係数」は、$\dbinom{2n}n\sim\dfrac{4^n}{\sqrt{\pi n}}$、すなわち
$$ \lim_{n\to\infty}\frac{\binom{2n}n}{4^n}\sqrt{\pi n}=1 $$
を示している。これを $u_{2n}$ の言葉で書くと
$$ \lim_{n\to\infty}u_{2n}\sqrt{\pi n}=1 $$
である。これは「$u_{2n}$ と $\dfrac1{\sqrt{\pi n}}$ の 比 が $1$ に近づく」という主張である(2 つの差も $0$ に近づくが、どちらも $0$ に近づくので、差が $0$ に近づくことだけでは大きさの見積もりにならない)。比 $u_{2n}\sqrt{\pi n}$ は $n=5$ で約 $0.975$、$n=50$ で約 $0.9975$ である(有限の $n$ の値は目安で、極限の証明は上の系による)。GS06 の第 12.1 節の練習問題 4 は、これを Stirling の公式から導く問題である。lem-rwr-bounds を書きかえると $\dfrac{\sqrt\pi}2\le u_{2n}\sqrt{\pi n}\le\sqrt{\dfrac{\pi n}{2n+1}}<\sqrt{\dfrac\pi2}$ で、比を $0.886$ 以上 $1.254$ 未満にはさむ評価になっている。極限の値 $1$ はこの間にある。

系:原点に戻る回数の期待値

原点に戻る確率が $1$ なら、1 回戻った後も同じことがくり返され、何回も戻ってきそうである。$2n$ 回目までに原点に戻る回数の期待値を求めて、これを確かめる。

原点に戻る回数の期待値

$2n$ 回目までに原点に戻る回数($S_k=0$ となる $k$、$1\le k\le2n$ の個数)を $N_{2n}$ とすると
$$ E[N_{2n}]=u_2+u_4+\cdots+u_{2n}=(2n+1)u_{2n}-1 $$
である。さらに $E[N_{2n}]>\sqrt n-1$ で、$n\to\infty$ のとき $E[N_{2n}]$ は限りなく大きくなる。

方針:$N_{2n}$ を「$k$ 回目に原点にいれば $1$、いなければ $0$」という確率変数の和に分け、期待値の線形性を使う。和の公式は $n$ についての帰納法で示す。
段 1(和に分ける)。$k=1,\ldots,2n$ について、$S_k=0$ なら $1$、そうでなければ $0$ をとる確率変数を $I_k$ とする。$N_{2n}=I_1+I_2+\cdots+I_{2n}$ である。$I_k$ の期待値は $1\cdot P(S_k=0)+0\cdot P(S_k\ne0)=P(S_k=0)$ である。確率変数の期待値と分散 の命題「期待値の線形性」により
$$ E[N_{2n}]=\sum_{k=1}^{2n}P(S_k=0) $$
である。prop-rwr-parity より奇数の $k$ の項は $0$ で、偶数の $k=2j$ の項は $u_{2j}$ なので、$E[N_{2n}]=u_2+u_4+\cdots+u_{2n}$ である。
段 2(和の公式)。$u_2+\cdots+u_{2n}=(2n+1)u_{2n}-1$ を $n$ についての帰納法で示す。$n=1$ では、左辺 $u_2=\dfrac12$、右辺 $3\cdot\dfrac12-1=\dfrac12$ で成り立つ。$n$ で成り立つとすると、$n+1$ では
$$ u_2+\cdots+u_{2n}+u_{2n+2}=(2n+1)u_{2n}-1+u_{2n+2} $$
である。prop-rwr-recur より $u_{2n+2}=\dfrac{2n+1}{2n+2}u_{2n}$、つまり $(2n+1)u_{2n}=(2n+2)u_{2n+2}$ なので、右辺は
$$ (2n+2)u_{2n+2}-1+u_{2n+2}=(2n+3)u_{2n+2}-1=\bigl(2(n+1)+1\bigr)u_{2(n+1)}-1 $$
となり、$n+1$ でも成り立つ。
段 3(限りなく大きくなる)。lem-rwr-bounds の下の評価 $u_{2n}\ge\dfrac1{2\sqrt n}$ から
$$ E[N_{2n}]=(2n+1)u_{2n}-1\ge\frac{2n+1}{2\sqrt n}-1>\frac{2n}{2\sqrt n}-1=\sqrt n-1 $$
である。$n\to\infty$ で $\sqrt n-1$ は限りなく大きくなるので、$E[N_{2n}]$ も限りなく大きくなる。$\square$

公式に数を入れる
  1. $n=2$:長さ $4$ の歩み $16$ 通りについて、原点に戻る回数を全部足す。2 回目に原点にいる歩みは、最初の 2 回が $+-$ か $-+$ で、残り 2 回は自由なので $2\cdot4=8$ 通り、4 回目に原点にいる歩みは ex-rwr-start の $6$ 通りである。戻る回数の合計は $8+6=14$ で、期待値は $\dfrac{14}{16}=\dfrac78$ である。公式でも $u_2+u_4=\dfrac12+\dfrac38=\dfrac78$、$(2\cdot2+1)u_4-1=\dfrac{15}8-1=\dfrac78$ で一致する。
  2. 公式で計算すると次の表になる。
    $n$$u_{2n}$$1-u_{2n}$$E[N_{2n}]$$\sqrt n-1$
    $5$$0.246$$0.754$$1.707$$1.236$
    $10$$0.176$$0.824$$2.700$$2.162$
    $50$$0.0796$$0.920$$7.039$$6.071$
    $100$$0.0564$$0.944$$10.326$$9.000$
    $100$ 回投げると平均で約 $7$ 回、$200$ 回投げると平均で約 $10$ 回原点に戻る。投げる回数を $2$ 倍にしても、戻る回数の期待値は約 $1.5$ 倍にしかならない。

Wallisの公式(高校数学) の系をもう一度使うと、$\dfrac{E[N_{2n}]+1}{2\sqrt{n/\pi}}=\dfrac{2n+1}{2n}\cdot u_{2n}\sqrt{\pi n}$ は $1$ に近づく($\dfrac{2n+1}{2n}\to1$ と、上の極限による)。つまり $E[N_{2n}]+1$ と $2\sqrt{n/\pi}$ の比が $1$ に近づき、戻る回数の期待値は $\sqrt n$ の程度で増える。GS06 の Example 12.3 は、同じ期待値を母関数を使って求めている。

Banach のマッチ箱問題との一致

Banachのマッチ箱問題 の定理「残りの本数の期待値」は、マッチが $N$ 本ずつ入った箱で、空に気づいたときにもう一方に残っている本数 $R$ の期待値を $E[R]=\dfrac{(2N+1)\binom{2N}N}{4^N}-1$ と求めている。これは cor-rwr-count の $(2n+1)u_{2n}-1$ で $n=N$ としたものと同じ式である。$N=n=1,\ldots,8$ について全部の歩みを数えて比べると、期待値だけでなく、「$2n$ 回目までに原点に戻る回数が $r$ である確率」と「残りの本数が $r$ である確率」$\dbinom{2n-r}n\left(\dfrac12\right)^{2n-r}$ が、すべての $r$ で一致した。一般の $n$ での分布の一致は、この記事では扱わない(証明しない)。

例と反例:条件を変えると何が変わるか

thm-rwr-main の (2) の等式「$2n$ 回目まで一度も原点に戻らない確率 $=$ $2n$ 回目に原点にいる確率」は、def-rwr-walk の条件と、回数が偶数であることに支えられている。

外す条件反例成り立たなくなること
条件 (i) の確率 $\frac12$(左右の確率を偏らせる。独立性は保つ)ex-rwr-cx-biased:表の確率 $p=\frac23$、$n=1$戻らない確率 $\frac59$ と原点にいる確率 $\frac49$ が等しくない。戻る回数の期待値も $8$ を超えない
条件 (i) の「$+1$ か $-1$」(その場にとどまる動きを許す。各動きの確率も $\frac13$ に変わる)ex-rwr-cx-lazy:$-1,0,+1$ を $\frac13$ ずつ、$n=1$戻らない確率 $\frac49$ と原点にいる確率 $\frac13$ が等しくない
回数が偶数 $2n$ であることex-rwr-cx-odd:回数 $3$3 回目まで戻らない確率は $\frac12$、3 回目に原点にいる確率は $0$
反例:表の確率が $\frac12$ でない

条件 (i) を変えて、$X_k=+1$ となる確率を $p=\dfrac23$、$X_k=-1$ となる確率を $q=\dfrac13$ とする(独立性はそのまま)。
(1) $n=1$(2 回)で比べる。2 回目まで一度も原点に戻らないのは $++$ と $--$ のときで、確率は $p^2+q^2=\dfrac49+\dfrac19=\dfrac59$ である。2 回目に原点にいるのは $+-$ と $-+$ のときで、確率は $2pq=\dfrac49$ である。$\dfrac59\ne\dfrac49$ なので、thm-rwr-main の (2) の等式は成り立たない。証明の段 3 では「どの歩みも同じ確率」を使って、確率を歩みの数に置きかえていた。$p\ne q$ では $+1$ の多い歩みほど起こりやすく、数えることと確率を求めることが一致しない。
(2) $2n$ 回目に原点にいる確率は、$+1$ と $-1$ が $n$ 回ずつの歩みがどれも確率 $p^nq^n$ で起こるので、$\dbinom{2n}n(pq)^n$ である。$\dbinom{2n}n$ は二項定理の $(1+1)^{2n}=4^n$ を展開した項の 1 つなので $4^n$ 以下で、
$$ \binom{2n}n(pq)^n\le(4pq)^n=\left(\frac89\right)^n $$
である。cor-rwr-count の段 1 と同じ計算で、$2n$ 回目までに原点に戻る回数の期待値は
$$ \sum_{k=1}^n\binom{2k}k(pq)^k\le\sum_{k=1}^n\left(\frac89\right)^k<\frac{8/9}{1-8/9}=8 $$
で、$n$ をどれほど大きくしても $8$ を超えない。$p=\dfrac12$ のときに限りなく大きくなったこと(cor-rwr-count)と対照的である。$p\ne\dfrac12$ では、いつか原点に戻る確率も $1$ より小さくなることが知られているが、この記事では証明しない。

反例:その場にとどまる動きを許す

条件 (i) を変えて、1 回の動きが $-1$、$0$、$+1$ のどれかで、それぞれ確率 $\dfrac13$ ずつとする(独立性はそのまま)。この点は奇数回目にも原点にいることがある(1 回目に $0$ が出れば $S_1=0$)ので、prop-rwr-parity も成り立たない。
2 回目までで比べる。2 回目に原点にいるのは、動きが $(0,0)$、$(+1,-1)$、$(-1,+1)$ の 3 通りで、確率は $\dfrac39=\dfrac13$ である。2 回目まで一度も原点に戻らないのは、1 回目が $\pm1$(2 通り)で、2 回目に原点に戻らない動き(位置 $1$ からなら $0$ か $+1$、位置 $-1$ からなら $0$ か $-1$ の 2 通り)のときで、$2\cdot2=4$ 通り、確率は $\dfrac49$ である。$\dfrac49\ne\dfrac13$ なので、等式は成り立たない。

反例:回数が奇数

def-rwr-walk のランダムウォークのまま、回数を奇数 $3$ にする。prop-rwr-parity より $S_3\ne0$ がつねに成り立つので、「3 回目まで一度も原点に戻らない」は「2 回目まで一度も原点に戻らない」と同じ出来事で、確率は thm-rwr-main の (2) により $u_2=\dfrac12$ である。一方、3 回目に原点にいる確率は $0$ である。等式を奇数の回数 $2n+1$ で使うときは、「$2n+1$ 回目まで一度も原点に戻らない確率は $u_{2n}$」と読みかえる。

大学数学で見る:母関数による証明と、高い次元

GS06 の第 12.1 節は、$2m$ 回目に初めて原点に戻る確率 $f_{2m}$ と $u_{2m}$ の関係(Theorem 12.2)を母関数の式に直して、$f_{2m}$ を求め(Theorem 12.3)、その母関数に $x=1$ を代入して、原点に戻る確率が $1$ であることを示している(Example 12.1)。この記事は母関数を使わず、鏡像原理(投票定理)と、積の上下の評価だけで同じ結論を得た。$f_{2m}$ は 初めて原点に戻る時刻 で、母関数を使わずに差をとる方法で求める。

平面や空間のランダムウォークについて開く

平面の格子の上を、上下左右に確率 $\frac14$ ずつ動く点や、空間の格子の上を 6 方向に確率 $\frac16$ ずつ動く点についても、同じ問いを考えられる。GS06 の Example 12.2 は、平面では原点に戻る確率が $1$ であり、空間では $1$ より小さいことを示している(GS06 はこれらの結果を Pólya によるものとしている)。この記事では証明しない。数直線と平面では「いつか必ず戻る」のに、空間では戻らないことがある、という違いは、ランダムウォーク で扱う話題である。

演習

$n=4$ で数える

長さ $8$ の歩みで、$S_1,\ldots,S_8$ がすべて正のものを、最後の位置 $S_8=2j$($j=1,2,3,4$)ごとに thm-rwr-main の証明の段 3 の式で数え、原点に戻らない歩みの数が $\dbinom84$ になることを確かめよ。

解答を開く

段 3 の式 $\frac j4\binom8{4+j}$ で、$j=1$:$\frac14\cdot56=14$、$j=2$:$\frac24\cdot28=14$、$j=3$:$\frac34\cdot8=6$、$j=4$:$\frac44\cdot1=1$。合計 $14+14+6+1=35$ で、ずっと負のものも $35$ 通りなので、原点に戻らない歩みは $70$ 通り。$\binom84=70$ と一致し、確率は $\frac{70}{256}=\frac{35}{128}=u_8$ である。

9 割の確率で原点に戻るには

$2n$ 回目までに原点に戻る確率が $0.9$ 以上になる最小の $n$ を求めたい。(1) lem-rwr-bounds を使って、$n\ge50$ なら確率が $0.9$ 以上であること、$n\le24$ なら $0.9$ 未満であることを示せ。(2) $u_{2n}$ を計算して、最小の $n$ を求めよ。

解答を開く

(1) 確率は $1-u_{2n}$ なので、$u_{2n}\le0.1$ かどうかを調べればよい。$n\ge50$ なら $2n+1\ge101>100$ で、$u_{2n}\le\frac1{\sqrt{2n+1}}<\frac1{10}$。$n\le24$ なら $u_{2n}\ge\frac1{2\sqrt n}\ge\frac1{2\sqrt{24}}>\frac1{2\cdot5}=0.1$($\sqrt{24}<5$)。

(2) prop-rwr-recur の漸化式で順に計算すると、$u_{62}\approx0.1009$($n=31$)、$u_{64}\approx0.0993$($n=32$)。よって最小の $n$ は $32$ で、$64$ 回投げれば、原点に戻る確率がはじめて $0.9$ を超える。

さらに先へ

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する