文献あり

k-リュカ数の四捨五入表示についての考察

はじめに

この記事ではk-リュカ数の四捨五入表示について私なりに考察した結果を紹介します。
かなり荒い評価ですがひとまず次のような結果が得られました。

$k$-リュカ数$L_n^{(k)}$の四捨五入表示$\L=\s{\a_k^n}$が$n\geq n_0$で成り立つような$n_0$のうち最小のものを$n_k$とおくと$n_k=O(k^3\log k)$が成り立つ。(つまり$k$の増加に伴って$n_k$が増大する速さは少なくとも$k^3\log k$以下であるということ。)

ただし$\a_k$は方程式$f_k(x)=x^k-\sum^{k-1}_{j=0}x^j=0$の唯一の正の実数解とした。
具体的に$n_k$がどんな数で抑えられるかは以下で示すとして、まず考察の過程を見ていきましょう。

方針

まず$k$-ナッチ数の四捨五入表示で有効だった確率論的アプローチを考えてみようとすると、これを適応するには少なくとも$k$個の連続する$n$について
$\dis|L^{(k)}_n-\a_k^n|<\frac12$
が成り立つことを示さなければなりませんが、今までわかっている$\a_k$の評価
$\dis 2(1-\frac{1}{2^k})<\a_k<2(1-\frac{1}{2^{k+1}-k})\quad$( 1 の補題6, 2 の命題4参照)
からは$n=k$前後で示したい不等式が成り立たないことがわかったのでひとまずこのアプローチは断念することにしました。
(実際に$\a_k^n$が$L_n^{(k)}=2^n-1\;(1\leq n\leq k)$よりもはみ出ていることをDesmosなどで確かめてみるとわかります。)

なので今度は方程式$f_k(x)=0$の$\a_k$ではない解$\b$の絶対値を評価してみることにしました。
$k$-リュカ数の一般項が
$\dis L^{(k)}_n=\a_k^n+\sum_{\ss{f_k(\b)=0\\\b\neq\a_k}}\b^n$
で表されること$|\b|<1$であったことを考慮するとある$n_0$に
$\dis|\b^{n_0}|<\frac{1}{2(k-1)}$
が成り立ってくれれば$n\geq n_0$において四捨五入表示が成り立つことになります。
そのような$n_0$を考えるにあたって以下の評価を得ることができました。

方程式$f_k(x)=0$の任意の$\a_k$ではない解$\b$について
$\dis|\b|<\sqrt[2k+1]{\frac{1}{5-4\cos(\frac{11\pi}{6k})}}$
が成り立つ。

この不等式を以下で示し、$n_0$をどのようにとればいいか考えていきましょう。

$|\b|$の評価

いま$\b$は方程式$\dis f_k(x)=\frac{x^{k+1}-2x^k+1}{x-1}=0$の解であったのでこれを適当に変形することで
$\dis \b^k=\frac{1}{2-\b}$
が成り立つことがわかります。ここで$\b=r(\cos\t+i\sin\t)$と極座標表示したとき上式の絶対値を取ると
$\dis r^k=\frac{1}{\sqrt{r^2-4r\cos\t+4}}$
が成り立ち、これをまた適当に変形することで方程式
$\dis\cos\t=\frac{1}{4r}(r^2+4-\frac{1}{r^{2k}})$
が得られます。この方程式から$r=|\b|$の大きさを考えてみることにしましょう。

まず$r=|\b|<1$において$r$が$\t$についての関数となることを示しましょう。
(Desmos等で$(x^2+y^2)^k((x-2)^2+y^2)=1$のグラフを見てみると想像が付きやすいと思います。)

$r$についての方程式
$\dis\cos\t=\frac{1}{4r}(r^2+4-\frac{1}{r^{2k}})$
は$0\leq r\leq1$においてただ一つに定まる。

方程式の左辺は$-1\leq\cos\t\leq1$であって
方程式の右辺を$g(r)$とおくと$\dis g(1)=1,\;\lim_{r\to0^+}g(r)=-\infty$および
$\dis g'(r)=\frac{1}{4}-\frac{1}{r^2}+\frac{2k+1}{4r^{2k+2}}=\frac{1}{4}+\frac{1}{r^2}(\frac{2k+1}{4r^{2k}}-1)\geq\frac{1}{4}+(\frac{2k+1}{4}-1)=\frac{k-1}{2}>0$
であることから増減表を考えることで二本のグラフ$y=\cos\t$と$y=g(x)$の交点は$0\leq x\leq1$においてただ一つ存在することがわかる。

次に$r$が$\t$の増減によってどのような挙動をするのか調べます。$r$は$\t$についての偶関数(グラフが$x$軸について対象)なので変数を$\t$の代わりに$\cos\t$を考えることにします。

$\cos\t$が増加するに伴って$r$も増加する。

補題3の証明から$\dis\frac{d\cos\t}{dr}=g'(r)>0$であったので$\dis\frac{dr}{d\cos\t}=\frac{1}{\frac{d\cos\t}{dr}}>0$を得る。

では$\cos\t$ひいては$\t=\arg\b$が取りうる範囲について考えてみましょう。

$\dis|\arg\b|\geq\frac{11\pi}{6k}$が成り立つ。

いま原点から円$(x-2)^2+y^2=1$に引ける接線は$\dis y=\tan(\pm\frac{\pi}{6})x$の$2$本であることに注意すると$|z|<1$において$\dis|\arg(2-z)|<\frac{\pi}{6}$が成り立ち、
また$\dis\b^k=\frac{1}{2-\b}$より$\dis\arg(\b)=\frac{2\pi j}{k}-\frac{1}{k}\arg(2-\b)$であって$j\neq0$である(？)(2( https://mathlog.info/articles/1395)の命題3参照)ので
$\dis|\arg\b|\geq\frac{2\pi}{k}-\frac{\pi}{6k}=\frac{11\pi}{6k}$
が成り立つ。

いま$\dis\cos(\frac{11\pi}{6k})\leq g(R)$なる$R$を取ってくると
補題5から$\dis g(r)=\cos\t\leq\cos(\frac{11\pi}{6k})\leq g(R)$、補題3の証明から$r=|\b|\leq R$
と$|\b|$を上から評価できるので命題$2$を得るには以下のことを示せばよいことになります。

$\dis R=\sqrt[2k+1]{\frac{1}{5-4\cos(\frac{11\pi}{6k})}}$とおくと$\dis\cos(\frac{11\pi}{6k})\leq g(R)$が成り立つ。

$r<1$において$\dis(\frac{r^2+4}{4r})'=\frac14-\frac{1}{r^2}<\frac14-1<0$であって
また$\dis\frac{1}{5-4\cos(\frac{11\pi}{6k})}<\frac{1}{5-4}=1$より$R<1$なので
$\dis g(R)=\frac{R^2+4}{4R}-\frac{1}{4R^{2k+1}}>\frac{1^2+4}{4\cdot1}-\frac{1}{4R^{2k+1}}=\frac54-\frac{5-4\cos(\frac{11\pi}{6k})}{4}=\cos(\frac{11\pi}{6k})$
を得る。

$n_0$の評価

さて上での議論により命題2
$\dis|\b|<\sqrt[2k+1]{\frac{1}{5-4\cos(\frac{11\pi}{6k})}}$
が成り立つことを示しましたがこの不等式からどのような$n_0$をとれば
$\dis|\b|^{n_0}<\frac{1}{2(k-1)}$
が成り立つのか見ていきましょう。

まず命題2の不等式の右辺がこのままでは扱いづらいので少し変形しましょう。
(流石に$\dis\frac{(2k+1)\log(2(k-1))}{\log(5-4\cos(\frac{11\pi}{6k}))}\leq n_0$の整数解を直接求めるのには無茶があります。)

$\dis|\b|<\sqrt[2k+1]{\frac{1}{1+\frac{24}{k^2}}}$が成り立つ。

いま任意に$\dis0< x\leq\frac{11\pi}{6}$を取ると
$\dis\cos(\frac{11\pi}{6k})\leq\cos(\frac{x}{k})\leq1-\frac{1}{2}(\frac{x}{k})^2+\frac{1}{24}(\frac{x}{k})^4$
が成り立つので$k\geq2$に注意すると
$\dis{5-4\cos(\frac{11\pi}{6k})}\geq{1+\frac{x^2}{k^2}(2-\frac{x^2}{6k^2})}\geq1+\frac{x^2}{k^2}(2-\frac{x^2}{6\cdot2^2})=1+\frac{x^2}{k^2}\cdot\frac{48-x^2}{24}$
となる。そして$\dis0<\sqrt{24}\leq\frac{11\pi}{6}$より$x=\sqrt{24}$とすると
$\dis\sqrt[2k+1]{\frac{1}{5-4\cos(\frac{11\pi}{6k})}}<\sqrt[2k+1]{\frac{1}{1+\frac{24}{k^2}}}$
を得る。

ここで$\dis\frac{1}{(1+\frac{24}{k^2})^{\frac{n_0}{2k+1}}}\leq\frac{1}{2(k-1)}$となるような$n_0$を考えるとこれを適当に変形することで
$\dis\frac{(2k+1)\log(2(k-1))}{\log(1+\frac{24}{k^2})}\leq n_0$
を得るので命題1を得るには以下のことを示せばよい。

$\dis\frac{(2k+1)\log(2(k-1))}{\log(1+\frac{24}{k^2})}=O(k^3\log k)$が成り立つ。

$k\to\infty$において
$2k+1=O(k)$
$\log(2(k-1))=O(\log k)$
$\dis\log(1+\frac{24}{k^2})=\sum^\infty_{j=0}\frac{1}{j}(\frac{24}{k^2})^j=O(\frac{1}{k^2})$
であることから主張を得る。

ちなみに$\cos$を展開しなくても$O(k^3\log k)$が示せるはずなのでこれが命題2から得られる最良の評価だと思います。
では具体的に$n_0$にどのような値が取れるか考えてみましょう。
整数を考えている以上$O(\log k)$をそのまま扱うのは難しいので$O(k)$で考えることにします。

$\dis n_0=\left\lceil\frac{(2k+1)k^3}{12}\right\rceil$とおくと$\dis\frac{1}{(1+\frac{24}{k^2})^{\frac{n_0}{2k+1}}}\leq\frac{1}{2(k-1)}$が成り立つ。

$\dis\frac{(2k+1)k^3}{12}\leq n_0$に注意すると
$\dis\frac{1}{(1+\frac{24}{k^2})^{\frac{n_0}{2k+1}}}\leq\frac{1}{(1+\frac{24}{k^2})^{\frac{k^3}{12}}}\leq\frac{1}{1+\frac{k^3}{12}\cdot\frac{24}{k^2}}=\frac{1}{1+2k}<\frac{1}{2(k-1)}$
を得る。

恐らくここでの議論も精密化することでもっと良い評価を得られると思いますがとりあえず私の考察としては以上にしたいと思います。