2
高校数学議論
文献あり

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

22
1
$$\newcommand{a}[0]{\alpha} \newcommand{b}[0]{\beta} \newcommand{dis}[0]{\displaystyle} \newcommand{L}[0]{L_n^{(k)}} \newcommand{l}[0]{\left} \newcommand{r}[0]{\right} \newcommand{s}[1]{\left\lfloor#1\right\rceil} \newcommand{t}[0]{\theta} $$

はじめに

 この記事では 前回の記事 に続いて$k$-リュカ数の四捨五入表示について私なりに考察した結果を紹介します。
 今回の考察では次のような結果が得られました。

 $k$-リュカ数$L_n^{(k)}$の四捨五入表示
$$L_n^{(k)}=\s{\a_k^n}$$
$n_0\leq n$で成り立つような$n_0$のうち最小のものを$n_k$とおくと
$$n_k=O(k^3)$$
が成り立つ。

 ただし$\a_k$は方程式
$$f_k(x)=x^k-\sum^{k-1}_{j=0}x^j=0$$
の唯一の正の実数解としました。
  前回 では$O(k^3\log k)$であったところが$O(k^3)$になったので多少ましになったと思います。

方針

 今回は 前回の記事 での考察を元にまず以下の結果を示します。

 方程式$f_k(x)=0$$\a_k$以外の解に適当に番号を付けて$\b_1,\b_2,\ldots,\b_{k-1}$とおくと
$$|\b_j|<\sqrt[2k+1]{\frac{1}{5-4\cos\frac{2\pi j}{k}}}$$
が成り立つ。

 これがわかれば$k$-リュカ数$\L$の一般項から
$$\l|L^{(k)}_{(2k+1)n}-\a_k^{(2k+1)n}\r| =\left|\sum^{k-1}_{j=1}\b_j^{(2k+1)n}\right|\leq\sum^{k-1}_{j=1}\frac{1}{(5-4\cos\frac{2\pi j}{k})^n}$$
が成り立つのであとはこれが$\frac12$未満になるためには$n$
$$n_0'=O(k^2)$$
以上であれば十分であることを示して命題1を得ます。

$\arg\b_j$の評価

 まず$\arg\b_j$の取りうる範囲について考えてみます。$\b$を方程式
$$f_k(x)=\frac{x^{k+1}-2x^k+1}{x-1}=0$$
の解とすると
$$\b^k=\frac{1}{2-\b}$$
が成り立つのでした。
 ここで$\b=r(\cos\t+i\sin\t)$と極座標表示すると上式は
$$r^k(\cos k\t+i\sin k\t)=\frac{(2-r\cos\t)+ir\sin\t}{(2-r\cos\t)^2+r^2\cos^2\t}$$
と表せ、この両辺の$(虚部)/(実部)$を取ることで
$$\tan k\t=\frac{r\sin\t}{2-r\cos\t}$$
つまり極方程式
$$r=\frac{2\tan k\t}{\sin\t+\cos\t\tan k\t}=\frac{2\sin k\t}{\sin(k+1)\t}$$
が得られます。この方程式から$\t$の取りうる範囲を考えてみましょう。
 いま$r$$\t$についての周期$2\pi$および偶関数(グラフが$x$軸について対称)なので$0\leq\t\leq\pi$で考えることにします。

$$g(\t)=\frac{2\sin k\t}{\sin(k+1)\t}$$
$0\leq\t\leq\pi$において単調増加である。

 和積の公式から
\begin{align} g'(\t) &=\frac{2(k\cos k\t\sin(k+1)\t-(k+1)\sin k\t\cos(k+1)\t)}{\sin^2(k+1)\t}\\ &=\frac{(2k+1)\sin\t-\sin(2k+1)\t}{\sin^2(k+1)\t} \end{align}
が成り立つことに注意する。
 いま$0\leq\t\leq\frac{\pi}{k+1}$において
\begin{align} ((2k+1)\sin\t-\sin(2k+1)\t)' &=(2k+1)(\cos\t-\cos(2k+1)\t)\\ &=2(2k+1)\sin k\t\sin(k+1)\t\\ &\geq0 \end{align}
より
$$(2k+1)\sin\t-\sin(2k+1)\t\geq(2k+1)\sin0-\sin0=0$$
が成り立ち、また
$$\arcsin\frac1{2k+1}<\frac\pi2\frac1{2k+1}<\frac\pi{k+1}$$
に注意すると$\frac\pi{k+1}\leq\t\leq\frac\pi2$においても
$$(2k+1)\sin\t-\sin(2k+1)\t >(2k+1)\sin\l(\arcsin\frac1{2k+1}\r)-1=0$$
が成り立つ。
 よって$g'(\t)=g'(\pi-\t)$に注意すると$0\leq\t\leq\pi$において$g'(\t)\geq0$が成り立つことがわかる。

 さてこの補題と$r=|\b|>0$であることから$\t$の取りうる範囲は以下のように評価できます。

 方程式$f_k(x)=0$$\a_k$以外の解に適当に番号を付けて$\b_1,\b_2,\ldots,\b_{k-1}$とおくと
$$\frac{2\pi j}{k}<\arg\b_j<\frac{\pi(2j+1)}{k+1}\quad\l(1\leq j<\frac{k}{2}\r)$$
が成り立つ。

 なお$\frac{k}{2}< j< k$においては$\b_j=\overline{\b_{k-j}}$と定めることで
$$-\frac{2\pi j}{k}>\arg\b_j>-\frac{\pi(2j+1)}{k+1}$$
が成り立ちます。

 まず$0<\t<\pi$において$g(\t)>0$なる$\t$の範囲を考える。
 いま$\frac{\pi j}{k+1}$の前後での$\sin k\t$および$\sin(k+1)\t$の符号を考えると
$$\lim_{\t\to\frac{\pi j}{k+1}\pm0}g(\t)=\mp\infty$$
なので$\frac{\pi j}{k+1}<\t<\frac{\pi(j+1)}{k+1}$における$g(\t)$のグラフは
$$\begin{array}{|c|c|c|c|c|c|c|c|} \hline\t&\frac{\pi j}{k+1}&\cdots&\frac{\pi j}{k}&\cdots&\frac{\pi(j+\frac12)}{k+\frac12}&\cdots&\frac{\pi(j+1)}{k+1} \\\hline g(\t)&-\infty&\nearrow&0&\nearrow&2&\nearrow&\infty \\\hline\end{array}$$
となり$\frac{\pi j}{k}<\t<\frac{\pi(j+1)}{k+1}$$g(\t)>0$の条件になる。
 ここで方程式$f_k(x)=0$$\a_k$以外の解
$$\b=r(\cos\t+i\sin\t)\quad(0< r<1,\ 0<\t<\pi)$$
は二つの極方程式
\begin{align} \cos\t&=\frac{1}{4r}\l(r^2+4-\frac{1}{r^{2k}}\r)\\ r&=\frac{2\sin k\t}{\sin(k+1)\t}\quad(=g(\t)) \end{align}
を満たすのであった。
 いまこのそれぞれの極方程式が表すグラフを考えてみると 前回の記事 の補題3や先の増減表から
$$\frac{\pi j}{k}<\t<\frac{\pi(j+1)}{k+1}$$
においてただ一つ交点を持ち($j=0$については$g(\t)>g(0)=\frac{2k}{k+1}>1$より不適)、また 前回の記事 の補題5から
$$\frac{2\pi j}{k}-\frac{\pi}{6k}<\arg\b<\frac{2\pi j}{k}+\frac{\pi}{6k}$$
であったことに注意すると$j$が偶数のときその交点は方程式$f_k(x)=0$の解となることがわかる。

 ちなみに増減表におまけとして示したように$r<1$という条件から
$$\frac{2\pi j}{k}<\arg\b_j<\frac{\pi(2j+\frac12)}{k+\frac12}$$
とも評価できる。

$n_0$の評価

 さて 前回の記事 と同じように
$$R_j=\sqrt[2k+1]{\frac{1}{5-4\cos\frac{2\pi j}{k}}}$$
とおくと
\begin{align} \frac{1}{4|\b_j|}\l(|\b_j|^2+4-\frac{1}{|\b_j|^{2k}}\r) &=\cos\t_j\\ &<\cos\frac{2\pi j}{k}=\frac54-\frac{1}{4R_j^{2k+1}}\\ &<\frac{1}{4R_j}\bigg(R_j^2+4-\frac{1}{R_j^{2k}}\bigg) \end{align}
つまり$|\b_j|< R_j$が成り立ち(これは$\frac{k}2< j< k$においても成り立つ)、$j=\frac{k}2$のときは$\arg\b_j=\pi$なので同様にして$|\b_j|< R_j$がわかる。以上より命題2を得る。
 命題2を得たいま
$$|L^{(k)}_{(2k+1)n}-\a_k^{(2k+1)n}|=\left|\sum^{k-1}_{j=0}\b_j^{(2k+1)n}\right|\leq\sum^{k-1}_{j=1}\frac{1}{(5-4\cos\frac{2\pi j}{k})^n}$$
が成り立つので次はこの右辺を評価していこう。

$$\frac{k^2}{\pi^2(1-\frac{\pi^2}{48})}< n$$
において
$$\sum^{k-1}_{j=1}\frac{1}{(5-4\cos\frac{2\pi j}{k})^n}<\frac12$$
が成り立つ。

 まず$1\leq j<\frac{k}4$において
$$\frac{1}{(5-4\cos\frac{2\pi j}{k})^n}<\frac{1}{4j(j+1)}$$
が成り立つことを示す。
 いま$j$の取り方より
\begin{align} 5-4\cos(\frac{2\pi j}{k}) &>5-4(1-\frac12\left(\frac{2\pi j}{k}\right)^2+\frac{1}{24}\left(\frac{2\pi j}{k}\right)^4)\\ &=1+8\pi^2\frac{j^2}{k^2}\l(1-\frac{\pi^2}{3}\cdot\frac{j^2}{k^2}\r)\\ &>1+8\pi^2\frac{j^2}{k^2}\bigg(1-\frac{\pi^2}{3}\cdot\frac{(\frac{k}4)^2}{k^2}\bigg)\\ &=1+8\pi^2\frac{j^2}{k^2}\l(1-\frac{\pi^2}{48}\r) \end{align}
なので$a=\pi^2(1-\frac{\pi^2}{48})$とおくと$\frac{k^2}{a}< n$より
\begin{align} \l(5-4\cos\frac{2\pi j}{k}\r)^n &>\l(1+8a\frac{j^2}{k^2}\r)^n\\ &>1+8an\frac{j^2}{k^2}\\ &>1+8j^2>4j(j+1) \end{align}
を得る。
 また$\frac{k}4\leq j\leq\frac{k}2$においても
$$\frac{1}{(5-4\cos\frac{2\pi j}{k})^n} \leq\frac{1}{5^{\frac{k^2}{a}}}<\frac{1}{4\frac k2(\frac k2+1)}\leq\frac{1}{4j(j+1)}$$
が成り立っているので$b=\lfloor\frac k2\rfloor$とおくと
\begin{align} \sum^{k-1}_{j=1}\frac{1}{(5-4\cos\frac{2\pi j}{k})^n} &\leq\sum^{b}_{j=1}\frac{2}{(5-4\cos\frac{2\pi j}{k})^n}\\ &<\frac12\sum^b_{j=1}\l(\frac1j-\frac1{j+1}\r)\\ &=\frac12\l(1-\frac1{b+1}\r)<\frac12 \end{align}を得る。

 以上より以下の主張を得る。

 少なくとも
$$\frac{k^2(2k+1)}{\pi^2(1-\frac{\pi^2}{48})}< n$$
において$L_n^{(k)}=\s{\a_k^n}$が成り立つ。

ちなみに$\pi^2(1-\frac{\pi^2}{48})=7.840\ldots$なのでざっくり$\frac{k^3}{3}\leq n$くらいで成り立つことになる。

オーダーの評価

 ここで$R_j$の取り方をもっとよくしたらオーダーが変わらないか考えてみよう。つまるところ
$$\cos\frac{2\pi j}{k}=\frac{1}{4R_j}\bigg(R_j^2+4-\frac1{R_j^{2k}}\bigg)$$
が成り立つように取ってみる。
 このとき
$$|\L-\a_k^n|<\sum^{k-1}_{j=1}R_j^n<\frac12$$
と評価できるとすると、少なくとも$R_1^{n}<\frac12$つまり
$$\frac{\log 2}{-\log R_1}< n$$
でなくてはならないので$R_1=r_k$とおいて$\log r_k$のオーダーを考えてみよう。

$$\log r_k=O\l(\frac1{k^3}\r)$$
が成り立つ。

 いま
\begin{align} r_k^2-\frac1{r_k^{2k}} &=e^{2\log r_k}-e^{-2k\log r_k}\\ &=(1+2\log r_k)-(1-2k\log r_k)+O((k\log r_k)^2)\\ &=2(k+1)\log r_k+O((k\log r_k)^2) \end{align}
が成り立つことに注意すると$r_k$の定義式
$$\cos\frac{2\pi}{k}=\frac{1}{r_k}\bigg(1+\frac14\Big(r_k^2-\frac1{r_k^{2k}}\Big)\bigg)$$
の対数を取ることで
\begin{align} \log\l(\cos\frac{2\pi}{k}\r) &=\log\bigg(1-\l(\frac{2\pi}k\r)^2+O\l(\frac1{k^4}\r)\bigg)\\ &=-\frac{1}{2}\left(\frac{2\pi}{k}\right)^2+O\l(\frac1{k^4}\r)\\ \log\bigg(1+\frac14\Big(r_k^2-\frac1{r_k^{2k}}\Big)\bigg) &=\frac{1}{4}\Big(r_k^2-\frac1{r_k^{2k}}\Big)+O\bigg(\Big(r_k^2-\frac1{r_k^{2k}}\Big)^2\bigg)\\ &=\frac{k+1}2\log r_k+O((k\log r_k)^2) \end{align}
つまり
$$-\frac{1}{2}\left(\frac{2\pi}{k}\right)^2+O\l(\frac1{k^4}\r) =\frac{k-1}2\log r_k+O((k\log r_k)^2)$$
と評価できることから主張を得る。

 したがってこの一連の手法では$n_k=O(k^3)$が最良の評価であることがわかる。

参考文献

投稿日:2021131
更新日:515
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

子葉
子葉
994
231012
主に複素解析、代数学、数論を学んでおります。 私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。 同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中