k-ナッチ数列の四捨五入表示の確率論的導出
はじめに
前回の記事
では
apu_yokai
氏の記事に触発されk-ナッチ数列の四捨五入表示の確かさをその一般項の公式から直接考察し、その後apu_yokai氏による
フィボナッチ数列を拡張したk-ナッチ数列の一般項についての予想(その4・証明完成)
の記事をもって完全解決に至りましたが、少し引っかかる点がありました。
というのもapu_yokai氏が立てた予想がWikipediaに2008年付の出典付きで既に紹介されていたわけなのです。
私も初めその出典を見たとき書いてある内容が不十分に見えてまーたWikipediaがいい加減な出典拾ってきたのかと半信半疑でしたが、今一度よくよく見てみると確かに四捨五入表示の導出になっていそうだということがわかってきたので私なりに言葉を付け足しつつその内容をまとめようと思います。
ざっくりとその内容を説明するととあるランダムウォークの問題を考えそれに対する期待値$e_n$を考えます。するとその期待値$e_n$が漸化式$e_{n+1}=2e_n-e_{n-k}$を満たすことがわかり、問題の条件を適当に設定することで$e_n$がk-ナッチ数$a_n$と$\a^n/f'(\a)$との誤差とみなせることを示します。そして$e_n$の確率的な解釈から$|e_n|<1/2$であることがわかり四捨五入表示問題の解決となります。
ランダムウォーク問題
最初に$x$軸上の非負整数点$x=n$に人が立っています。この人は$\frac12$の確率で$x$軸正の方向に$1$進み、$\frac12$の確率で負の方向に$k$進むという動作を区間$[0,k-1]$内に至るまで無限に繰り返します。そしてその人が止まった点$x=j\;(0\leq j\leq k-1)$に応じてその人は得点$e(j)$を得るものとします。
このときその人が得る得点の期待値$e_n$、つまりその人が点$x=j$に止まる確率$p_n^{(j)}$について
$e_n=p_n^{(0)}e(0)+p_n^{(1)}e(1)+\cdots+p_n^{(k-1)}e(k-1)$
と定められる値$e_n$が上述のような性質を満たすことを示します。
漸化式$e_{n+1}=2e_n-e_{n-k}$が成り立つ。
各$j$に対して$p^{(j)}_{n+1}=2p^{(j)}_n-p^{(j)}_{n-k}$が成り立つことを言えばよい。
いま$x=n$からスタートして一回目の動作でそれぞれ$\frac12$の確率で$x=n+1$または$x=n-k$の位置に進み、そこから$x=j$に至る確率はそれぞれ$p^{(j)}_{n+1},\;p^{(j)}_{n-k}$なので
$\dis p^{(j)}_n=\frac{p^{(j)}_{n+1}+p^{(j)}_{n-k}}{2}$
つまり$p^{(j)}_{n+1}=2p^{(j)}_n-p^{(j)}_{n-k}$を得る。
さてこの漸化式から方程式$\dis\frac{x^{k+1}-2x^k+1}{x-1}=0$の解を$\b_1,\b_2,\ldots,\b_k$とおくと$e_n$および$p_n^{(j)}$の一般項は$1,\b_1^n,\b_2^n,\ldots,\b_k^n$の線形結合で表されることがわかる。ここで$\b_k=\a$を唯一の正の実数解とするとそれらの線形結合において$\a^n$の係数が$0$であることを示そう。
$e_n$の一般項は$1,\b_1^n,\b_2^n,\ldots,\b_{k-1}^n$の線形結合で表される。
いま$|\b_i|<1\;(1\leq i< k)$および$1<\a$に注意すると、$p_n^{(j)}$における$\a^n$の係数が$0$でないとすると$n\to\infty$において$p_n^{(j)}$は正か負の無限大に発散することになるが$p_n^{(j)}$は確率なので$0\leq p_n^{(j)}\leq1$であることに矛盾。よって$p_n^{(j)}$ひいては$e_n$における$\a^n$の係数は$0$であり主張を得る。
四捨五入表示の導出
いま得点$e(j)$の値については特に定めてなかったので、$a_n$をk-ナッチ数列、$f(x)$を$f(x)=\frac{x^{k+1}-2x^k+1}{x-1}$とし$e(j)$を
$\dis e_n=a_n-\frac{\a^n}{f'(\a)}\;(n=0,1,2,\ldots,k-1)$
となるように定めることにする。具体的には$p_j^{(j)}=1$より$e_j=e(j)$なので$e(j)=a_j-\frac{\a^j}{f'(\a)}$となる。
このとき$e_n$の一般項を
$e_n=c_0+c_1\b_1^n+c_2\b_2^n+\cdots+c_{k-1}\b^n_{k-1}$
と表すと$c_0,c_1,\ldots,c_{k-1}$についての$k$個の方程式
$\dis e_n=a_n-\frac{\a^n}{f'(\a)}=\sum^{k-1}_{j=1}\frac{\b_j^n}{f'(\b_j)}\;(n=0,1,2,\ldots,k-1)$
から$c_0=0,c_j=\frac1{f'(\b_j)}$がわかり、つまり$k\leq n$においても常に
$\dis e_n=\sum^{k-1}_{j=1}\frac{\b^n_j}{f'(\b_j)}=a_n-\frac{\a^n}{f'(\a)}$
が成り立つことになる。
ところで
前回の記事
で示したように$0\leq n< k$において
$\dis|e(n)|=\left|a_n-\frac{\a^n}{f'(\a)}\right|<\frac12$
が成り立つことがわかっているので$k\leq n$においても
$\dis\left|a_n-\frac{\a^n}{f'(\a)}\right|=|e_n|\leq\sum^{k-1}_{j=0}p_n^{(j)}|e(j)|<\frac12\sum^{k-1}_{j=0}p_n^{(j)}=\frac12$
となることがわかる。
以上が確率論的な四捨五入表示の導出となります。私たちとは違うアプローチであるとは言えど既に2008年にはれっきとした証明が付けられていたみたいですね。
この記事を高評価した人
この記事に送られたバッジ
投稿者
